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performance  evaluation  for  Rayleigh  fading  channels.  The  expressions  for  the 
probability  of  error  for  these  systems  are  formulated  in  a  way  that  allows  the 
identification  of  the  key  parameters  of  the  communication  system. 

The  performance  of  DPSK  and  FSK  are  evaluated  for  several  signaling  formats. 
The  design  parameters  considered  are  the  shapes  of  the  basic  data-pulse  wave¬ 
forms  for  both  DPSK  and  FSK,  the  modulation  index,  and  the  relative  phase  between 
successive  transmitted  signals  for  FSK.  It  is  shown  that  the  system  error 
probability  is  highly  dependent  on  these  system  parameters. 

The  interaction  between  the  various  elements  of  the  communication  system 
is  examined.  We  show  that  the  average  error  probability  can  be  approximated  in 
terms  of  one  or  two  rms  type  channel  measurements.  A  technique  for  obtaining 
bounds  on  system  performance  in  terms  of  the  key  system  and  channel  parameters 
is  described  and  applied  to  the  evaluation  of  both  DPSK  and  FSK  systems  for 
several  channel  models  and  signaling  formats.  A  method  of  approximating  the 
performance  of  systems  employing  complicated  pulse  shapes  and  channels  which  are 
difficult  to  fully  characterize  is  discussed.  It  is  shown  that  the  techniques 
for  obtaining  bounds  and  approximations  for  Rayleigh  channels  are  easily  applied 
to  the  more  general  frequency-selective  Rician  fading  channels. 

The  applicability  of  adaptive  equalization  techniques  to  digital  communica¬ 
tions  over  WSSUS  frequency-selective  fading  channels  is  discussed.  We  describe 
the  characteristics  of  adaptive  equalizers  that  are  commonly  employed  for 
fading-channel  communications.  It  is  found  that  adaptive  equalizers  can  be  used 
to  establish  a  coherent  communications  environment  as  well  as  to  reduce  the 
effects  of  ISI.  The  results  of  cited  simulation  and  experimental  studies  are 
compared  to  the  analytical  results  for  the  error  probabilities  of  the  DPSK  and 
FSK  systems.  We  develop  a  method  of  obtaining  estimates  of  adaptive  equalizer 
performance  for  practical  systems. 
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ABSTRACT 


The  performance  of  digital  communications  over  selective  wide-sense- 
stationary  nncorrelated-scattering  (WSSUS)  fading  channels  is  investigated. 
The  emphasis  is  on  the  evaluation  of  the  performance  of  binary  differential 
phase-shift  keyed  (DPSK)  and  binary  frequency-shift  keyed  (FSK)  communications 
with  primary  emphasis  on  the  effects  of  intr?symbol  interference  produced  by 
the  frequency-selective  character  of  WSSUS  fading  channels. 


The  error  probabilities  of  DPSK  and  FSK  are  evaluated  for  several  models 
of  WSSUS  frequency-selective  fading  channels.  The  analysis  is  focused  on  the 
performance  evaluation  for  Ravleiah  fading  channels.  The  expressions  for  the 
probability  of  error  for  these  systems  are  formulated  in  a  way  that  allows  the 
identification  of  the  key  parameters  of  the  communication  system. 


The  performances  of  DPSK  and  FSK  are  evaluated  for  several  signaling 
formats.  The  design  parameters  considered  are  the  shapes  of  the  basic  data- 
pulse  waveforms  for  both  DPSK  and  FSK,  the  modulation  index,  and  the  relative 
phase  between  successive  transmitted  signals  for  FSK.  It  is  shown  that  the 
system  error  probability  is  highly  dependent  on  these  system  parameters. 

The  interaction  between  the  various  elements  of  the  communication  system 
is  examined.  We  show  that  the  average  error  probability  can  be  approximated 
in  terms  of  one  or  two  rms  type  channel  measurements.  A  technique  for 


obtaining  bounds  on  system  performance  in  terms  of  the  key  system  and  channel 
parameters  is  described  and  applied  to  the  evaluation  of  both  DPSK  and  FSK 
systems  for  several  channel  models  and  signaling  formats.  A  method  of 
approximating  the  performance  of  systems  employing  complicated  pulse  shapes 
and  channels  which  are  difficult  to  fully  characterize  is  discussed.  It  is 
shown  that  the  techniques  for  obtaining  bounds  and  approximations  for  Rayleigh 
channels  are  easily  applied  to  the  more  general  frequency-selective  Rician 
fading  channels. 

The  applicability  of  adaptive  equalization  techniques  to  digital 
communications  over  VSSUS  frequency-selective  fading  channels  is  discussed. 
Ve  describe  the  characteristics  of  adaptive  equalizers  that  are  commonly 
employed  for  fading-channel  coauuunications.  It  is  found  that  adaptive 
equalizers  can  be  used  to  establish  a  coherent  communications  environment  as 
well  as  to  reduce  the  effects  of  ISI.  The  results  of  cited  simulation  and 
experimental  studies  are  compared  to  the  analytical  results  for  the  error 
probabilities  of  the  DPSK  and  FSK  systems.  Ve  develop  a  method  of  obtaining 
estimates  of  adaptive  equalizer  performance  for  practical  systems. 
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CHAPTER  1 


INTRODUCTION 

Several  problems  arise  in  the  consideration  of  slow- frequency- hopped 
(SFH)  spread-spectrum  communications  which  motivate  the  study  of  the 
performance  of  digital  communications  over  selective  fading  channels  [1-4]. 
In  many  applications  of  SFH  systems,  the  channel  (or  transmission  medium) 
cannot  be  adequately  modeled  as  a  non-dispersive  additive  white  Gaussian  noise 
channel.  In  cases  where  the  characteristics  of  the  channel  are  significantly 
different  from  this  ideal,  the  channel  is  commonly  referred  to  as  a  fading 
channel.  Fading  channels  may  exhibit  such  undesirable  properties  as  a  time- 
varying  amplitude  response,  the  spreading  of  transmitted  signals  in  the 
frequency  domain  (time-selectivity),  and  dispersion  in  time  (frequency- 
selectivity)  which  may  produce  significant  intersymbol  interference  (ISI).  In 
cases  where  it  is  impractical  to  obtain  accurate  channel  estimates  and 
incorporate  these  estimates  in  the  detection  process,  the  random  character  of 
the  fading  channel  precludes  the  use  of  coherent  demodulation.  Binary 
differential  phase-shift  heying  (DPSK)  and  frequency- shift  keying  (FSK)  are  of 
particular  interest  for  applications  of  SFH  systems  in  selective  fading 
channels,  since  these  forms  of  digital  communications  do  not  require  the 
receiver  to  establish  phase  coherence  at  the  beginning  of  each  hop  [2,3,5]. 
The  primary  focus  of  this  thesis  is  the  evaluation  of  the  average  probability 
of  error  for  DPSK  and  FSK  communications  in  selective  fading  channels.  The 
analyses  of  these  systems  are  of  considerable  importance  independent  of 
applications  to  spread-spectrum  communications.  Furthermore,  the  results  on 
the  average  error  probabilities  for  DPSK  and  FSK  can  also  be  used  in  the 
performance  evaluation  of  SFH  systems  [2-4], 
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Previous  analyses  of  various  foras  of  digital  communications  over  fading 
channels  which  are  selective  only  in  time  show  that  system  performance 
parameters  such  as  signal-to-noise  ratio  and  average  error  probability  are  not 
significantly  degraded  unless  the  degree  of  time  selectivity  is  quite  large 
[2,3,6-12].  When  coupled  with  the  measurements  of  the  time-selectivity  of 
practical  channels  [7,13],  these  results  indicate  that  the  time-selective 
nature  of  the  channel  is  not  a  limiting  factor  for  the  performance  of  these 
systems.  In  this  thesis,  we  are  primarily  concerned  with  the  effects  of 
intersymbol  interference  produced  by  the  frequency-selective  character  of 
fading  channels.  Thus,  models  of  strictly  frequency-selective  channels  are 
used  in  the  evaluation  of  the  average  probability  of  error. 

Two  common  examples  where  fading  phenomena  are  encountered  in  practice 
are  ionospheric  high-frequency  "skywave"  and  tropospheric  scatter  channels 
[6-8].  A  channel  model  which  accurately  describes  the  characteristics  of 
these  and  other  examples  of  fading  environments  is  the  wide-sense-stationary 
uncorrelated-scattering  (WSSUS)  fading  channel,  which  is  discussed  at  length 
in  [9]  and  [10].  This  model  is  quite  general  and  includes,  for  example,  the 
doubly-selective  Rician  channel  as  a  special  case. 

Experimental  investigations  [8,14,13]  of  the  multipath  characteristics  of 
fading  channels  show  that  no  one  multipath  model  adequately  describes  the 
properties  of  various  frequency-selective  channels  encountered  in  practice. 
These  investigations  also  indicate  that  there  are  many  situations  where  the 
multipath  parameters  of  the  channel  do  not  remain  constant  during  the  time 
required  for  transmission  of  a  long  data  sequence.  In  Chapter  2,  four 
examples  of  multipath  propagation  models  for  WSSOS  frequency-selective 
channels  are  described. 
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In  Chapters  3  and  4,  ve  evaluate  the  average  probabilites  of  error  for 
binary  DPSK  and  FSK  over  WSSUS  frequency-selective  fading  channels.  The 
analysis  is  focused  on  Ravleiah  fading  channels,  since  the  adverse  effects  of 
frequency-selective  fading  are  most  evident  for  such  channels.  The  primary 
goal  of  the  analysis  in  Chapters  3  and  4  is  to  formulate  expressions  for  the 
system  error  probability  that  provide  insight  to  the  nature  of  the  fading 
mechanism.  The  expressions  for  the  probability  of  error  for  these  systems  ar 
developed  in  a  vay  that  allows  the  identification  of  the  key  parameters  of  tl 
communication  system. 

One  of  the  fundamental  design  parameters  for  DPSK  systems  is  the  shape  of 
the  data-pulse  waveform.  Since  the  spectral  characteristics  of  the 
transmitted  signal  are  largely  determined  by  the  properties  of  this  waveform, 
it  is  reasonable  to  expect  that  the  performance  of  DPSK  in  a  frequency- 
selective  fading  environment  is  highly  dependent  on  the  choice  of  the  pulse 
waveform.  For  FSK  communications,  the  character  of  the  transmitted  signal 
depends  on  the  frequency  separation  between  the  two  FSK  tones  and  the  relative 
phases  between  successive  transmitted  signals  as  well  as  the  shape  of  the 
data-pulse  waveform.  Previous  investigations  of  DPSK  and  FSK  over  frequency- 
selective  fading  channels  [2,5,16-18]  indicate  that  the  error  probabilities 
are  strongly  dependent  on  these  parameters. 

Ve  examine  the  interaction  between  the  characteristics  of  the  fading 
channel,  the  transmitted  signal,  and  the  nonlinear  detection  that  arises  in 
both  DPSK  and  FSK  demodulation.  We  show  that  in  many  cases  of  practical 
interest,  the  performance  of  these  systems  can  be  approximated  in  terms  of  one 
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or  two  parameters  which  can  be  obtained  from  rms  type  channel  measurements 
(e.g.,  see  [13]).  It  is  found  that  the  error  probabilities  for  DPSK  and  FSK 
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depend  on  a  number  of  common  factors  and  that,  in  certain  cases,  the  effects 
of  frequency  selectivity  on  both  DPSK  and  FSK  systems  can  be  characterized  by 
the  same  channel  measurement.  A  technique  for  obtaining  bounds  on  system 
performance  is  described  and  applied  to  the  evaluation  of  both  DPSK  and  FSK 
systems  for  several  examples  of  channel  models  and  signaling  formats.  We  then 
present  a  method  for  approximating  the  performance  of  systems  employing 
complicated  pulse  shapes  and  channels  that  are  difficult  to  fully 
characterize.  By  using  the  results  of  [19],  it  is  shown  that  the  techniques 
for  obtaining  bounds  and  approximations  for  Rayleigh  channels  can  easily  be 
applied  to  the  more  general  frequency-selective  Rician  fading  channels. 

In  many  practical  systems,  the  effects  of  intersymbol  interference  (ISI) 
in  a  frequency-selective  fading  channel  can  severely  limit  the  performance  of 
conventional  digital  communications.  For  a  number  of  years,  considerable 
attention  has  been  given  to  adaptive  signal  processing  or  "equalization” 
techniques  for  digital  cosimuni cat ions  over  certain  non-fading  channels  [20- 
21].  For  channels  of  this  type,  such  as  telephone  lines  and  line-of-sight 
(LOS)  microwave  links,  both  linear  and  nonlinear  equalizers  (typically  in  the 
configuration  of  tapped-delay-line  (TDL)  filters)  have  been  used  effectively 
to  reduce  the  effects  of  ISI  [22].  More  recently,  there  has  been  considerable 
interest  in  applications  of  similar  equalization  techniques  to  improve  the 
performance  of  digital  communications  over  WSSUS  frequency-selective  fading 
channels  [23].  While  the  random  character  of  fading  channels  presents  a 
number  of  additional  difficulties  which  must  be  overcome,  both  simulation 
studies  [24]  and  experimental  evidence  [14]  indicate  that  reliable 
communications  can  be  achieved  even  in  fading  environments  that  produce 
unacceptably  high  error  probabilities  for  unequalized  systems. 


In  Chapter  5,  ve  describe  adaptive  equalizers  that  are  commonly  employed 
for  fading-channel  communications.  It  is  found  that  the  main  attributes  of 
adaptive  equalizers  are  the  reduction  of  the  effects  of  ISI  and  the  ability  to 
establish  coherent  communications.  Ve  briefly  discuss  the  results  of  several 
simulation  studies  and  experimental  investigations.  These  results  are 
compared  to  the  results  presented  in  Chapters  3  and  4.  By  examining  the  basic 
properties  of  TDL  equalizers,  ve  develop  a  method  of  obtaining  estimates  of 
adaptive  equalizer  performance  for  practical  systems. 
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CHAPTER  2 

CHANNEL  MODELS 

In  this  chapter,  ve  discuss  the  key  characteristics  of  the  wide-sense¬ 
stationary  oncorrelated-scattering  (WSSUS)  fading  channel  model  used  in  the 
performance  evaluation  of  the  communication  systems  discussed  in  subsequent 
chapters.  This  channel  model  is  described  in  detail  in  [9]  and  is  employed  in 
the  analysis  of  a  variety  of  digital  communication  systems  in 
[2,5,9,11.16,17,24,25]. 

Ve  employ  narrowband  signal  models  (see  [6])  so  that  if  the  input  to  the 
channel  is 


s(t)  -  Re{s(t)exp(j2nf  t)}  , 

c 

then  the  output  is  given  by 


(2.1) 


where 


r(t)-=  Re(r(t)exp(j2nf  t)}  , 

c 


:(t)  =  as(t)  +  J  h(t,5)s(t-5)  d$  +  n(t)  , 


(2.2) 


(2.3) 


and  n(t)  is  the  equivalent  (see  [11])  low-pass  Gaussian  noise  with  (one-sided) 
spectral  intensity  Nq.  For  the  general  case  of  Rician  fading,  the  received 
signal  r(t)  consists  of  three  components:  a  single  specular  component,  a 
diffuse  or  Rayleigh-faded  component,  and  the  channel  noise.  The  fading 
process  is  characterized  by  the  ensemble  autocovariance  of  the  response 
function  h(t,g) , 

E[h(t,x)h*(x.{)]  -  2a2p(t-x,{)6(x-$)  .  (2.4) 
where  6(»)  is  the  Dirac  delta  function  and  2o^  is  the  total  power  in  the 
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fading  process,  so  that  J  p(0,$)  d{  ■  1.  Alternatively,  the  quantity  o2  can 

&GD 

be  interpreted  as  the  average  power  received  when  a  (real)  sinusoid  of  unity 
peak  value  is  transmitted  over  the  channel. 

If  a  >  0  in  (2.1),  there  is  a  specular  (non-faded)  component  present  in 
the  channel  output.  In  this  case  the  channel  is  termed  a  Rician  fading 
channel  (as  in  [12,19]).  If  a  =  0,  there  is  no  specular  component  in  the 
received  signal  and  the  channel  is  a  Rayleigh  fading  channel  as  is  considered 
in  [2,5,11,16,17].  The  autocovariance  function  in  (2.4)  represents  the 

second-order  statistics  of  a  channel  that  is  selective  in  both  time  and 
frequency  (i.e.,  doubly  selective).  Donbly  selective  channels  are  the  most 
general  examples  in  the  class  of  WSSOS  fading  channels  [9]. 

The  effects  of  time-selective  fading  on  DPSK  communications  are 
considered  in  [3]  for  several  data  pulse  shapes  and  a  variety  of  fading 

channel  models.  In  [11],  the  performance  of  binary  FSK  communication  via 
time-selective  fading  channels  is  calculated  for  one  example  of  a  fading 
channel.  In  each  case  the  performance  is  evaluated  as  a  function  of 

parameters  related  to  the  bandwidth  of  the  Doppler  power  density  spectrum 

[11].  While  it  is  probably  true  that  most  fading  channels  are  time-selective 
to  some  degree,  it  is  also  true  that,  for  a  given  fading  channel,  the 
normalized  Doppler  spread  (see  [3])  decreases  with  increasing  channel  data 
rate.  Moreover,  the  results  in  [11]  and  [3],  when  coupled  with  measured 
Doppler  spreads  for  typical  channels  [7], [13],  indicate  that  the  time- 
selective  nature  of  the  channel  is  not  a  limiting  factor  in  determining  the 
performance  of  practical  systems.  In  the  analyses  that  follow,  we  consider 
channels  that  are  selective  only  in  frequency,  i.e.,  that  are  strictly 


frequency-selective , 


For  strictly  frequency-selective  fading  channels,  the  antocovariance 
function  (2.4)  of  the  fading  process  becomes  [9] 

2a2  p(t-r.?)6(r-5)  =  2<r2g($)8 (t-^)  ,  (2.5) 

where  g(£)  is  the  delay  oower-densitv  spectrum  of  the  fading  process.  The 
inverse  Fourier  transform  of  g($)  is  called  the  frequency  correlation  function 
[9]  and  is  given  by 

4D 

G(fl)  -  j  g($)ej2na5  d?  .  (2.6) 

£«D 

There  are  several  ways  in  which  the  degree  of  frequency-selectivity  can 
be  specified.  One  is  to  define  the  selectivity  as  the  "bandwidth"  of  the 
frequency  correlation  function  G(0) .  For  example,  in  [16]  and  [17],  the 
distance  between  the  "1/e"  points  of  G(Q)  is  used  as  a  measure  of  frequency- 
selectivity.  Alternatively,  the  degree  of  selectivity  can  be  defined  in  terms 
of  the  delay  power-density  spectrum  g($). 

In  any  physical  channel  the  transmitted  signal  undergoes  a  propagation 
delay  x d  (say)  which,  for  the  analyses  of  non-selective  or  non-fading 
channels,  is  usually  assumed  to  be  a  known  deterministic  quantity  that  is 
compensated  for  at  the  receiver.  In  the  case  of  frequency-selective  channels, 
the  transmitted  signal  may  experience  a  continuum  of  random  propagation 
delays,  with  mean  value  given  by 

*d  “  J  UU>  d$  .  (2.7) 

• €0 

which  is  referred  to  as  the  mean  path  delay.  A  parameter  that  is  commonly 
used  to  specify  the  selectivity  of  the  channel  in  terms  of  relative  delay  is 
the  "multipath  spread"  [5]  which  is  defined  by 


B.  S. 


2{L  U  '  5d)2g<?)  d?}1/2 


(2.8) 


(see  [13],  where  techniques  for  measuring  this  channel  parameter  are 
discussed) . 

We  are  concerned  with  relative  delays  only  and  we  assume  that  the 
receiver  has  compensated  for  the  mean  path  delay  Equivalently,  we  let 

*  0  in  order  to  establish  a  time  reference.  Parameters  that  are  used  in  the 
sequel  to  specify  the  degree  of  frequency-selectivity  are  the  rms  delay. 

JO 

*  “  if  $2*lU)  d$}1/2  (2.9) 

which  is  equal  to  half  the  spread  parameter  defined  in  [5]  and  the  rms 
multipath  spread,  defined  by  p  **  M/T,  where  T  is  the  data  symbol  duration. 
This  latter  parameter  is  related  to  the  normalized  data  rate  [16,17],  since  it 
is  a  function  of  the  ratio  of  the  transmitted  data  rate  and  the  channel 
correlation  bandwidth. 

Since  the  delay  power-density  spectrum,  g(£)  in  (2.5)  (and  hence  the 
function  p(0,{)  in  (2.4))  can  be  viewed  as  the  Fourier  transform  of  the 
correlation  function  G(Q)  in  (2.6),  it  is  necessarily  a  nonnegative,  real¬ 
valued  function.  We  also  assume  that  g(£)  is  symmetric,  i.e.,  g(£)  =  g(~£)> 
This  assumption  can  be  made  without  loss  of  generality  since  we  consider  only 
symmetric  binary  signaling,  and  the  average  probability  of  error  depends  on 
the  fading  channel  model  only  through  an  ensemble  average  (2.4)  of  the  channel 
statistics. 

In  the  analysis  that  follows,  the  evaluation  of  the  average  probability 
of  error  for  DPSK  and  FSK  is  considered  for  four  examples  of  delay  power- 
density  spectra:  the  Gaussian,  exponential,  triangular  and  rectangular  delay 


10 


The  parameter  N  in  (2.10)  end  (2.11)  is  the  ras  delay  given  by  (2.9).  The 
exponential  delay  power-density  spectrum  is  given  by 

g(?)  *  *  ,  (2.12) 

and  the  corresponding  frequency  correlation  function  is 

G(fi)  -  (1  +  2(nMfl)2)-1.  (2.13) 


The  triangular  delay  spectrum,  which  was  employed  in  [4]  and  [18]  for  the 
analysis  of  a  slow-frequency-hopped  FSE  multiple-access  system,  is  given  by 


Ul  <  T0. 

otherwise  . 


(2.14) 


The  frequency  correlation  function  for  the  triangular  spectrum  is 

G(Q)  -  sinc2(0T0)  ,  (2.15) 

where  sinc(x)  ■  sin(nx)/(nx) ;  and  the  ras  delay  is  given  by  M  *  TQ/  'T&.  The 
rectangular  delay  power-density  spectrum  is  given  by 


Li 
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%  s 

p 

?:  " 


z  3 
&  v> 


"ft  d 


V 

I  S 

1  § 


ft  M 


<2T0)_1  ;  |?l  <  Tq 


(2.16) 


;  otherwise  , 


with  corresponding  frequency  correlation  function 


G(fl)  -  sinc(2flT0)  , 


(2.17) 


where  M  -  T 0/VT. 


For  a  general  frequency-selective  channel,  the  detection  of  a  given 
information  bit  may  depend  on  a  number  of  consecutive  data  pulses,  because  of 
intersymbol  interference  (ISI).  Although  the  analysis  presented  below  is 
easily  extended  to  include  more  severe  intersymbol  interference,  we  assume  as 
in  [S,16.17]  that  the  degree  of  frequency-selective  fading  is  small  enough 
that  the  intersymbol  interference  affects  only  adjacent  data  pulses;  this  is 
if  p(t-x,£)  •»  0  for  |$|  >  T  (cf.  (2.4)).  In  this  case,  the  channel 
is  referred  to  as  an  adiaccnt-nulso-lfnited  ISI  channel.  This  assumption  is 
made  for  two  reasons.  If  the  intersymbol  interference  is  not  limited  to 
adjacent  data  pulses,  it  is  much  more  difficult  to  formulate  tractable 
expressions  for  the  system  error  probability.  Second,  for  all  models  of  the 
delay  power-density  spectra  considered  here,  the  complementary  assumption 
implies  that  the  resulting  average  probability  of  error  is  unacceptably  large. 

Our  assumption  of  adjacent— pulse-limited  ISI  places  a  restriction  on  the 
maximum  value  of  the  rms  delay.  In  particular,  if  we  require  that  at  least 
90%  of  the  total  energy  of  the  delay  spectrum  lies  within  the  range 
t~T  j  1  T]  for  the  Gaussian  delay  density,  (i.e.,  if  we  require  that  in  the 
absence  of  additive  noise,  less  than  10%  of  the  energy  received  in  the 
interval  [(i-l)T, iT]  is  due  to  pulses  transmitted  outside  the  interval 


[(i-2)T,  (i+l)T]) ,  then  ire  must  have 


*  *  m  f  <  0.43  . 

where  |i  is  the  rms  multipath  spread.  For  the  exponential  delay  power-density 
spectrum  the  requirement  that  90%  of  the  total  energy  of  the  delay  spectrum  is 
within  the  range  [-T  i.  (  i.  X]  implies  that  the  rms  multipath  spread  satisfies 


4 


ln(.l) 


0.62  . 


In  the  case  of  a  triangular  delay  spectrum,  intersymbol  interference  is 
completely  limited  to  adjacent  pulses  if 

2  Zo 

H  “  gf  i  1/6. 


Finally,  it  ia  easy  to  see  that  intersymbol  interference  for  the  rectangular 
spectrum  is  limited  to  adjacent  pulses  if 

2  Ifi. 

*  *  3T  1/3  • 

In  the  analysis  that  follows,  we  show  that  the  "shape"  of  the  delay 
power-density  spectrum  (as  well  as  the  rms  multipath  spread)  can  have 
considerable  influence  on  the  error  probability  of  both  differentially 
coherent  and  noncoherent  communi cat ions. 


CHAPTER  3 


PERFORMANCE  OF  BINARY  DPSK  COMMUNICATIONS 


OVER  FREQUENCY- SELECTIVE  FADING  CHANNELS 


I  * 

V 

y  ? 

■3  > 

'•  v 


In  this  chapter,  we  consider  the  evaluation  of  the  average  probability  of 
error  for  binary  DPSK  communi cat ions  over  wide-sense-stationary  uncorrelated- 
scattering  frequency-selective  fading  channels.  In  what  follows,  we  primarily 
consider  the  performance  of  DPSK  over  WSSUS  frequency-selective  Ravleiah 
fading  channels.  It  is  in  this  case  that  the  adverse  effects  of  frequency- 
selective  fading  and  the  dependence  of  system  performance  on  the  channel  delay 
power-denaity  spectrum  is  most  evident. 


The  dependence  of  the  error  probability  for  Rayleigh  fading  on  the  shape 
of  the  data-pulse  waveform  used  as  the  DPSK  signal  has  been  considered  in 
[3, Id, 17].  In  particular,  the  rectangular  pulse  as  well  as  the  sine  pulse, 
which  is  the  basic  pulse  shape  for  min imua- shift-keying  (MSK)  modulation,  is 
considered  in  [3]  where  the  average  probability  of  error  is  obtained  for  the 
Gaussian,  exponential,  triangular,  and  rectangular  delay  power-density 


spectrum  models.  The  result  for  DPSK  using  the  rectangular  pulse  for  a 
Gaussian  delay  spectrum  has  been  previously  obtained  by  Bello  and  Nelin  [17] 
and  by  Bailey  and  Lindenlaub  [16].  In  [16],  the  authors  also  consider  the 
rectangular  pulse  and  the  raised-cosine-spectrum  pulse  for  a  rectangular  delay 
power-density  spectrum. 

A  number  of  conclusions  may  be  drawn  from  these  results: 
i)  For  all  models  of  delay  power-density  spectra  considered,  significant 
gains  in  performance  can  be  achieved  by  proper  choice  of  pulse  shape. 
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ii)  Closed-fora  expressions  for  the  average  probability  of  error  can  be 
unwieldy  even  for  the  simplest  pulse  shapes  because  of  the  quadratic 
nature  of  interference  inherent  in  differential  detection  in  a 

frequency-selective  environment. 

iii)  Depending  on  the  data  pulse  and  the  choice  of  normalization  for  multipath 
spread,  the  error  probability  can  be  very  sensitive  to  the  parameters  of 
the  fading  channel,  i.e.,  delay  power-density  spectrum. 

The  first  observation  indicates  that  the  error  probability  of  DPSK  in 
frequency-selective  fading  channel  can  be  substantially  reduced  by  a  judicious 
choice  of  pulse  waveform.  Unfortunately,  the  latter  observations  suggest  that 
identification  of  a  "good"  pulse  shape  is,  at  best,  analytically  cumbersome 
and  dependent  on  a  number  of  channel  and  system  parameters  (as  indicated  by 
the  complicated  form  of  the  results  in  [5,16,17]).  Moreover,  it  is  not 
necessarily  true  that  analyses  of  this  type  provide  a  reasonable  indication  of 
the  performance  for  a  physical  channel  unless  an  exact  mathematical 
description  of  the  channel  is  available.  This  is  especially  unfortunate  since 
a  complete  characterization  may  not  be  possible  [13],  and  for  many  practical 
channels,  the  characteristics  are  not  likely  to  remain  constant  during  the 
time  required  for  transmission  of  a  long  data  sequence  [9].  Hence,  it  is 
possible  that  several  statistical  models  could  be  used  to  describe  the  same 
frequency-selective  Rayleigh  channel  for  repeated  transmissions. 

In  this  chapter,  we  show  that  the  performance  of  DPSK  communications  over 
frequency-selective  Rayleigh  fading  channels  can  be  closely  approximated  in 
terms  of  one  or  two  parameters  which  can  be  obtained  from  rms-type  channel 
measurements.  A  technique  for  obtaining  bounds  on  system  performance  in  terms 
of  the  key  channel  parameters  is  described.  We  present  a  method  for 


approximating  the  performance  of  systems  employing  complicated  pulse  shapes 
and  channels  that  are  difficult  to  characterize  fully.  Finally,  it  is  shown 
using  the  results  of  [19],  that  our  results  on  bounds  and  approximations  of 
the  probability  of  error  for  Rayleigh  channels  can  be  applied  to  the  more 
general  case  of  frequency-selective  Rician  fading  channels. 

3.1  Svatem  Model 

The  system  consists  of  an  information  source,  a  DPSK  transmitter,  and  a 
differentially  coherent  receiver.  The  information  to  be  transmitted  is 
modeled  as  a  sequence  (b^)  of  mutually  independent  random  variables,  each 
taking  values  0  or  1  with  equal  probability.  The  binary  data  sequence  (b^), 
with  elements  in  {-1.1},  is  formed  by  differentially  encoding  the  information 
sequence  (b^).  Thus,  (bj)  is  a  sequence  of  mutually  independent  random 
variables,  each  taking  on  the  values  -1  and  +1  with  equal  probability.  The 
data  signal  is  given  by 


b(t)  =  b.v(t-iT) 

for  each  integer  i.  The  data  pulse-waveform  v(t)  is  assumed  to  be  time 
limited  to  the  interval  [0,T]  such  that 


=  1  ; 


thus  for  each  integer  i,  the  data  signal  is  a  positive  or  negative  version  of 
the  basic  pulse  shape  v(t).  The  transmitted  signal  s(t)  is  given  by  (2.1) 
with  s(t)  defined  by 


r  4  -T  »  a  w  £  ^  4  ■ 
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s(t) 


-  [JK 
\J  T 


b(t)  . 


(3.1) 


where  E  is  the  energy  per  data  bit  of  the  transmitted  signal. 


In  the  present  analysis,  three  examples  of  the  pulse  waveform  v(t)  are 
considered.  These  are  the  rectangular  pulse 


v(t)  -  pT(t) 


1.  0  <  t  <  T 

0,  otherwise  , 


(3.2) 


the  sine  pulse. 


v(t)  -  \/T  sin(nt/T)pT(t)  , 
and  the  rectangular  phase-coded  pulse  [26]  given  by 


(3.3) 


v(t)  *  2  ‘i  *T  <t"iTc) 
i-0  c 


(3.4) 


for  0  ^  t  i  T.  The  sequence  (a^),  referred  to  as  the  signature  sequence,  is  a 
sequence  of  elements  of  (+1,-1).  The  chip  duration  Tc  is  related  to  T  by 
T  "  NTc  where  N  is  the  integer  number  of  chips  per  data  pulse.  The  sine  pulse 
in  (3.3)  is  the  basic  waveform  used  for  minimum-shift-keying  systems  while  the 
phase-coded  pulse  is  referred  to  as  the  spectral-spreading  signal  in  direct- 
sequence  (DS)  or  code-division  spread-spectrum  multiple-access  (SSMA) 
communications  [26-29]  and  hybrid  SFH/DS  SSMA  [4] . 


The  DPSK  receiver  is  the  differentially  coherent  matched  filter  receiver 
shown  in  Fig.  3.1.  At  the  end  of  the  i-th  data-pulse  interval  the  receiver 
forms  the  decision  statistic  represented  (in  terms  of  narrowband  models)  by 


r  f(i+l)T  fiT  *  -| 

2Re  I  4  r(t)v*(t)dt  f  r"(t)v(t)dt  1 

L  JiT  J(i-1)T  J 


UV*  +  0*V 


(3.5) 


m 


m 

£ 


* 


>,  A 
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where 


D  =  2 


V  =  2 


Not  i 


.( i+l)T 

* 

iT 

r( t) v  (t)  dt 

(3.6a) 

■iT 

(i-l)T 

r(t)v  (t)  dt  . 

(3.6b) 

tic  Z . 

l 

is  expressed  as  a  quadratic 

form  of 

U  and 

V. 

Under  the  assumption  that  p(t-x,?)  *  0;  |£|  >  T,  i.e.,  adjacent-pulse- 
limited  ISI,  the  output  statistic  depends  on  at  most  four  consecutive  data 
bits  represented  by  =  ^i-2'  ^i-1*  ^i'  ^i+1^  an<*  t^ie  probability  of  error 
can  be  written  in  terms  of  probabilities  conditioned  on  the  event  that  certain 
sequences  of  data  bits  were  transmitted,  i.e.. 


p.  -  u  ?  p<ii>  ■ 


ii 


(3.7) 


where  P(b^.)  =  Pr  (error  occurslb^  transmitted}.  In  [11,17]  Bello  and  Nelin 
show  that  these  conditional  error  probabilities  are  given  by 


-  (2  +  y^))'1  . 


(3.8) 


where  y(*)  is  the  "equivalent"  signal-to-noise  ratio,  which  may  be  written  in 
terms  of  the  conditional  moments  ^(bj)  =  E(XY*lb.}  of  the  random  variables  U 
and  V  as 


2mUV(fej) 


r<ii>  - 


*  mUU^i^”VV^i^  “  “UV^i^ 
_ _2mUV^i^ 


;  bi-i  =  bi 


(3.9) 


'Z  “uu^i” 


;  bi  *  b^-, 


®VV^i^  +  “UV^i^ 


.v 


- 


v"v 

$ 

■  c 


Using  (1.5)  and  (3.6),  the  necessary  moments  are  written  as 


,r(i+l)T 

f(  i+l)T  -T 

8o 

J  J  P(t 

J  iT 

J  iT  *T 

+  8Nqt 

(3.10a) 


( i+l)T  iT  T 

“^(b . )  =  8o2  f  f  f  p(t-*.{)s(t-5)v*(t)s*(x-§)v(x)  d(dxdt.  (3.10b) 

UV  1  JiT  J ( i-l)T  -T 


The  moment  myy(b^)  is  defined  similarly  to  ®pxj(b£)  with  the  region  of 
integration  of  the  outer  two  integrals  given  by  [(i-l)T,  iT] ,  and  uy^t^)  = 

“WV  * 

By  substituting  g(>)  for  p(<.<)  in  (3.10)  and  using  the  assumption  that  g 
is  symmetric,  the  moments  needed  to  evaluate  r(b^)  are  found  to  be 

2 

“^(b^  =  Hll  lny  +  T,y  +  ibi_1bi  +  bib^)^]  +  8N0T,  (3.11a) 

2 

“w^j)  -  +  %  +  ^i-2bi-i  +  bi-lbi)T»v]  +  8N0T-  (3  .lib) 
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mUV(^i)  =  — ~  [  2bibi-lT'v  +  (bi-2bi-l  +  bibi+l},lv 


where 


*  <2  *  bi-2l>i  *  »i-ib1+i)<J  . 


\  =  J  g<5)5j<?) 

x 

%  =  J  gdt)Rv<5)  as. 


K  =  J0  g(5)Rv<5)  RVU>  , 


(3.11c) 

(3.12a) 

(3.12b) 

(3.12c) 


where  the  parameters  Ry($)  and  ^(5)  are  defined  in  (3.14).  By  examining 

(3.9)  and  (3.11),  it  is  easy  to  see  that  y(iA)  exhibits  certain  symmetry 

* 

relations  in  terms  of  the  data  pulse  sequence  b^.  Making  use  of  these 
properties,  P  in  (3.7)  can  be  reduced  to 

v 

Pe  =  8  [P(-l.+l.-l.+D  +  P(+l, +1.-1.-1)  +  2P(-1, +1.-1.-1) 

+  P(-l,+l,+l,-l)  +  P(+l,+l,+l,+l)  +  2P(-1 ,+l ,+l ,+l) 1 .  (3.13) 


The  functions  R  (•)  and  R  (•)  in  (3.12),  which  are  the  usual  continuous- 
v  v 

time  aperiodic  autocorrelation  functions  [26,27]  for  time-limited  signals, 
depend  only  on  the  pulse  waveform  v(t).  They  are  defined  by 


and 


f  v(t)  v(t-£)  dt 


*VU)  ■  f  v<*>  v(t+T-«)  dt 
J0 


(3.14a) 


(3.14b) 


for  0  ^  <  T.  These  functions  are  easily  evaluated  for  the  first  two  pulse 

shapes.  They  are  found  [27]  to  be 


RyU)  =  T  -  5  (3.15a) 

and 

Rv<$)  =  (3.15b) 

for  the  rectangular  pulse,  and  they  are  given  by 

Ry(0  =  (T-5)cos(n5/T)  +  (T/n)sin(n5/T)  (3.16a) 

and 


Ry($)  *  -5  cos(rt5/T)  +  (T/n)sin(n£/T) 


(3.16b) 
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for  the  sine  pulse.  The  autocorrelation  functions  R^f.)  and  Rv(«) 
phase-coded  pulse  (3.4)  are  given  for  0  <.  kT  <  £  <  (k+l)T  <  T  by 

Rv(£)  =  Ca(k)((k+l)Tc  -  £)  +  ca(k+l)(£  -  kTc) 

and 

=  Ca(N-k)((k+l)Tc  -  £)  +  Ca(N-k-l)(£  -  kTc>  . 

where  is  the  discrete  aperiodic  autocorrelation  function  for  the 
(a^)  of  length  N  defined  in  [26]  by 

N-l-k 

(  1  aiai+k'  0  <  k  <  N-l 

i=0  1  * 

J  N-l+k 

C  (k)  -  \  I  a j_ka • ,  1-N  <  k  <  0 

i-0 

v  0,  Ikl  >.  N  . 

This  function  iS  considered  in  [26—29]  and  many  of  its  properties  are  given 
in  [30].  In  what  follows,  it  is  shown  that  the  performance  of  the  phase-coded 
pulse  depends  on  the  correlation  properties  of  the  signature  sequence  (a.)  as 
well  as  other  system  characteristics  discussed  earlier.  Notice  that  the 
simple  rectangular  pulse  (3.2)  is  obtained  as  a  special  case  of  (3.4)  for  any 
N  2  1  by  letting  a.  =  +l  for  all  i  €  {0,  ....  N-l).  We  would  also  point  out 
that  the  approximation  [2,3]  of  the  average  error  probability  for  SFH/SSMA  in 
terms  of  the  hopping  pattern  parameters  and  the  probability  of  error  given 
there  are  no  hits  does  not  generally  apply  to  phase-coded  pulse  waveforms. 
Indeed,  the  principal  reasons  for  considering  this  class  of  pulse  shapes  for 
coherent  PSK  systems  include  increased  multiple-access  capability  [4,26-29] 
and  immunity  to  the  interference  effects  caused  by  multiple  and/or  diffuse 


(3.17b) 

sequence 


for  the 


(3.17a) 


propagation  paths  [12,31].  Since,  at  present,  there  is  no  evidence  to  the 
contrary,  one  would  expect  similar  gains  to  be  realized  in  differentially 
coherent  systems.  Hence,  a  proper  analysis  of  differentially  coherent  hybrid 
slow-frequency-hopped  direct-sequence  SSMA  would  be  influenced  by  these 
characteristics. 

Notice  that  the  moments  necessary  to  specify  the  average  probability  of 
error  depend  on  the  pulse  shape  and  the  delay  power-density  spectrum  only 
through  the  integral  expressions  (3.12).  Thus,  the  average  probability  of 
error  is  easily  obtained  once  these  three  basic  integrals  are  evaluated.  In 
[3]  closed-form  expressions  for  the  integrals  in  (3.12)  are  obtained  for  the 
rectangular  and  sine  pulses,  and  for  the  four  examples  of  delay  spectra 
discussed  in  Chapter  2. 

From  (3.9)  and  (3.11)  it  is  clesr  that  the  probability  of  error  P  for  a 

c 

WSSUS  frequency-selective  fading  channel  depends  on  a  number  of  system  and 
channel  parameters.  In  order  to  evaluate  the  effect  of  the  fading  process,  it 
is  instructive  to  first  consider  a  limiting  case  in  which  the  additive  channel 
noise  is  neglected.  In  particular.  Fig.  3.2  shows  Pe  as  a  function  of  the 
signal-to-noise  ratio  S  =*  2<t^E/Nq,  which  is  the  ratio  of  the  average  energy 
per  data  pulse  of  the  received  signal  to  the  real  noise  power  spectral  density 

Nq/2.  The  rectangular  pulse  waveform  and  the  Gaussian  model  for  the  delay 
power-density  spectrum  with  rms  multipath  spread  p  =  0.05  is  assumed  for  this 
example.  Also  shown  in  Fig.  3.2  is  Pn8,  the  error  probability  for  DPSK  in  a 
nonselective  Rayleigh  fading  channel,  which  is  [6,  Eq.  9-5-25], 

?n#  -  (2  +  4o2E/N0)_1  .  (3.18) 


Notice  that  when  the  signal-to-noise  ratio  S  is  small,  the  average  probability 


of  error  is  well  spproziaated  by  PQ#.  This  is  because  at  low  signal-to-noise 
ratios,  the  decision  errors  are  largely  dne  to  additive  Ganssian  noise  rather 
than  intersymbol  interference. 

For  large  signal-to-noise  ratios,  the  errors  caused  by  the  effects  of 

intersyabol  interference  doainate;  the  average  error  probability  approaches  a 

liaiting  error  probability  Pj  as  s  ->  «.  This  liaiting  error  probability  is 

known  as  the  irrednc ible  error  probability  for  the  frequency-selective  WSSUS 

Rayleigh  fading  channel  [5,16,17].  For  large  signal-to-noise  ratios,  the 

irreducible  error  probability  is  a  good  (although  pessimistic)  estimate  of  the 

average  probability  of  error.  If  Pj  and  Pns  are  known,  we  can  obtain  a  good 

approximation  to  the  average  probability  of  error  for  both  large  and  small 

signal-to-noise  ratios.  Moreover.  Pj  alone  is  a  useful  indication  of 

performance  for  practical  systems  since  it  is  a  lower  bound  on  P  for  all 

© 

signal-to-noise  ratios  and  a  good  estimate  of  P  for  high  signal-to-noise 
ratios. 

In  [5]  the  results  of  numerical  evaluations  of  Pj  =  Pe|g  _y  m  are  given 
for  the  rectangular  and  sine  pulses  and  the  four  delay  power-density  spectra 
discussed  above.  The  parameter  d,  given  by 

d  =  |  {  J  52gU>  d?}1/2  .  (3.19) 

is  used  as  a  basis  for  comparing  the  system  performance  for  the  various 
channel  models.  This  parameter  is  just  twice  the  rms  multipath  spread  p  =  M/T 
where  M  is  the  rms  delay  given  by  (2.9). 

In  Fig.  3.3,  the  irreducible  error  probabilities  Pj.  are  shown  as  a 
function  of  the  rms  multipath  spread  p  for  the  various  combinations  of  basic 
pulse  shapes  and  delay  power-density  spectra  mentioned  above.  The  values  of 


Pj  for  the  phase-coded  pulse  are  calculated  by  standard  numerical  integration 
techniques.  The  signature  sequence  used  for  this  waveform  is  an  m-sequence  of 
length  N  -  31  in  its  characteristic  phase  (see  definition  and  Fig.  7  in  [32]). 
Unless  stated  otherwise,  comparisons  between  irreducible  error  probabilities 
are  always  made  with  respect  to  a  fixed  value  of  rms  multipath  spread. 

Notice  that  while  the  irreducible  error  probabilities  for  the  rectangular 
and  sine  pulse  are  "well-behaved"  functions  of  the  rms  multipath  spread,  the 
limiting  error  probability  for  the  phase-coded  pulse  varies  in  a  more  erratic 
way  making  a  relative  evaluation  of  the  effects  of  intersymbol  interference 
difficult.  In  the  next  section,  we  show  that  such  erratic  behavior  is 
characteristic  of  the  irreducible  error  probability  for  phase-coded  pulses 
when  m-sequence s  (or  any  binary  sequences  that  produce  large  time-bandwidth 
product  signals)  are  used  as  signature  sequences.  These  results  do,  however, 
indicate  that  rectangular  pulse  DPSK  exhibits  small  variations  in  the 
irreducible  error  probability  with  respect  to  changes  in  the  model  for  the 
delay  power-density  spectrum.  In  contrast,  the  irreducible  error 
probabilities  for  systems  employing  the  sine  pulse  and  the  phase-coded  pulse 
vary  by  more  than  a  factor  of  10.  In  the  next  section,  it  is  shown  that  this 
apparent  insensitivity  of  the  rectangular  pulse  to  changes  in  the  shape  of  the 
delay-spectrum  is  due  primarily  to  the  choice  of  delay-spread  normalization 
rather  than  any  inherent  properties  of  the  pulse  itself.  For  example,  if  the 
normalization  is  defined  in  terms  of  the  "1/e  points"  of  the  corresponding 
frequency  correlation  function  6(0) ,  as  in  [1(5,17],  the  irreducible  error 
probabilities  for  the  rectangular  pulse  vary  by  a  factor  of  four.  In  the  next 
section,  we  examine  the  implications  of  the  choice  of  normalization  on  the 
evaluation  of  the  irreducible  and  average  error  probabilities. 


In  order  to  gain  a  better  understanding  of  the  relationships  between  data 
pnlse  shapes,  models  for  the  delay  power-density  spectrum,  the  signal-to-noise 
ratio,  and  other  system  design  parameters,  it  is  useful  to  examine  the 
frequency-selective  fading  mechanism  from  a  "cause-and-ef feet"  point  of  view. 
In  this  section,  we  investigate  the  influence  of  the  fading  process  on  the 
average  error  probability  by  separately  considering  the  effects  of  WSSUS 
frequency— selective  fading  in  the  time  domain  (or  delay  domain)  and  in  the 
Fourier  transform  or  frequency  domain.  It  is  true,  of  course,  that  these 
representations  each  merely  describe  the  same  fading  mechanism  from  different 
points  of  view.  However,  the  separate  consideration  of  these  two  effects 
provides  insight  to  the  causes  of  degradation  of  system  performance  inherent 
in  a  frequency-selective  environment.  Moreover,  investigating  the  fading 
characteristics  in  this  way  aids  in  the  identification  of  the  nature  of  the 
fundamental  trade-offs  between  the  key  system  parameters.  In  the  remainder  of 
this  section,  we  first  consider  the  characterization  of  frequency-selective 
fading  in  the  delay  domain  and  identify  the  key  parameters  for  determining  the 
effects  of  intersymbol  interference.  (In  all  that  follows,  we  assume  that  the 
delay  density  g(()  is  symmetric.)  We  then  focus  on  the  frequency-domain 
characterization  and  identify  the  system  parameters  that  provide  an  indication 
of  performance  degradation  due  to  additive  noise. 

3*2.1  Effects  £f  Intersymbol  Interference :  Characterization 

Consider  the  detection  of  an  information  bit  corresponding  to  the  data- 
bit  pair  ^Q.bj).  If  we  ignore  the  contributions  of  the  desired  signal  as 


well  as  the  additive  noise  and  only  consider  the  effects  of  intersymbol 


interference  resulting  from  positive  (relative)  propagation  delays,  the 
decision  statistic  Zj  =  uv*  +  U*V  with  the  random  variables  D  and  V  given  by 


and 


0  =  2  J^.v(t~T)  j  Mt-T.t)  /^b0v(t-5)d5dt 


v  "  2  f  v(t)  f  h(t.5>  y?b-iv(t+T-5)d5dt  . 

"o  *o  * 


(3.20a) 


(3.20b) 


Using  (2.4),  (2.5)  and  (3.20),  the  moments  necessary  for  the  evaluation 
of  the  average  probability  of  error  are  found  as 


m  16o2e 

“Vu  =  myy  =  — -  -  nv 

T 


(3.21a) 


and 


“uv 


16o2E 

T 


b_ibQ  nv  • 


(3.21b) 


Under  the  assumption  that  the  delay  power-density  spectrum  g(£)  is  symmetric, 
it  is  easy  to  see  that  identical  expressions  result  from  the  consideration  of 
the  effects  of  intersymbol  interference  produced  by  negative  values  of 
relative  delay  (with  a  change  of  data-bit  indices  in  (3.21b)).  Alternatively, 
if  we  consider  the  matched  filter  output  for  both  positive  and  negative 
relative  delays  when  b^  =  (+1 ,+l ,+l ,+l) ,  the  moments  in  (3.11)  are  found  to  be 

2  2 

®UU  *  ®w  ”  “uv  =  — ~ —  2tnv  +  ^v  +  2l*v^  =  — Y  E,Hv  (3.22) 

so  that  in  the  absence  of  additive  noise,  (3.22)  represents  the  total  "power" 
at  the  filter  output  due  to  the  transmitted  signal  with  no  data  modulation. 

Thus,  2i)y/2(ny  +  +  2t|^)  ■  2ny/Hy  represents  the  portion  of  filter  output 

due  to  interfering  signals  relative  to  the  total  output.  Hence,  the  parameter 
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^/Ily  provides  a  characterization  of  the  effects  of  intersymbol  interference 
on  DPSK  communications  in  frequency-selective  fading. 

From  the  discussion  in  Section  3.1  it  is  apparent  that  the  parameter  most 
indicative  of  the  effects  of  intersymbol  interference  is  the  irreducible  error 
probability  Pj.  Hence  we  would  expect  that  if  the  delay-spread  is  normalized 
with  respect  to  2qv/Hv,  the  limiting  error  probability  Pj  should  exhibit 
minimal  sensitivity  to  the  actual  shape  of  the  delay  power-density  spectrum 
g(£).  In  particular,  notice  that,  under  the  assumption  of  adjacent-pulse- 
limited  ISI,  the  rms  mnltipath  spread  can  be  written  as 

T 

M  =  Pr  =  {  I  f  *<«>*!<«>  d«>1/2 

r  J0  r 

“  T  {  2  ^r}1/2  ’  (3.23) 

since  =  j2  for  the  rectangular  pulse,  where  Rr(£)  =  R^)  is  given  by 

(3.15).  (The  subscript  "r"  denotes  the  rectangular  pulse.)  Thus,  it  is  not 
'surprising  that  the  normalized  rms  mnltipath  spread  provides  an  effective 
normalization  of  the  various  delay  spectra  for  the  rectangular  pulse  as 
suggested  by  the  results  in  Fig.  3.3. 

Similarly,  we  define  the  normalized  rms  multipath  spread  ps  for  the  sine 
pulse  by 

^S  =  l  i  JT  «u>r2u>  w1/2 

s  J0  s 

*  1  I  -V1/2  '  (3.24) 

s 

with  R#($)  *  Ry($)  given  by  (3.16b).  In  Table  3.1,  numerical  values  of  Hg/T2 
are  listed  for  the  four  examples  of  delay  power-density  spectra  discussed  in 
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Chapter  2  and  various  values  of  rms  delay  M.  It  is  easy  to  see  that  I^/t2  2 

0.6  for  all  values  of  rms  delay,  and  that  is  not  very  sensitive  to  the 

shape  of  the  delay  spectra.  Moreover,  H  =  for  all  practical  values  of  rms 

s 

delay  so  that  qg  alone  provides  a  good  indication  of  the  effects  of 
intersymbol  interference  for  the  sine-pulse. 

Finally,  the  normalized  rms  multipath  spread  for  the  phase-coded  pulse 
is  given  by 


T 

*T  ~  <  1  l„  «‘5>8pU>  IS) 


1/2 


P  '0 

- { I -v1/2 

HP  P 


(3.25) 


*ith  RpU)  *  RyU)  given  by  (3.17b).  The  values  of  Hp/T2  for  the  m-sequence 
of  length  N  «  31  are  listed  in  Table  3.2  as  a  function  of  the  rms  delay  M. 


In  contrast  to  the  results  for  the  rectangular  and  sine  pulse,  the 

parameter  is  very  sensitive  to  the  shape  of  the  delay  spectrum  and  can  be 

much  smaller  than  T^  for  relatively  small  values  of  rms  delay.  This  is  true 

for  two  reasons.  The  parameter  q  can  be  substantially  smaller  than  the 

corresponding  parameters  for  either  the  rectangular  or  sine  pulse;  this  fact 

has  profound  implications  on  the  performance  of  phase— coded  pulse  waveforms  in 

the  presence  of  additive  noise.  Using  the  properties  of  the  discrete  periodic 

autocorrelation  function  for  m-sequences  [30],  it  is  easy  to  show  that  q^ 

given  by  (3.12c)  is  never  positive.  In  fact,  depending  on  the  shape  of  the 

delay  density,  the  negative  contribution  of  2q^  can  be  nearly  as  large  as  the 

sum  of  q^  and  This  is  in  sharp  contrast  to  the  analogous  situation  for 

the  rectangular  and  sine  pulse  where  q^  is  non— negative  and  generally  much 

smaller  than  q  , 
v 


S'. 

upl 


The  negative  contributions  of  in  the  absence  of  data  modulation  can  be 
vieved  as  intrasvmbol  interference  which  is  characteristic  of  phase-coded 
pulses  when  m-sequences  are  used  as  signature  sequences.  While  this 
characteristic  is  not  exclusive  to  m-sequences,  it  is  not  generally  true  for 
other  classes  of  binary  sequences.  Because  of  this  intrasymbol  interference, 
phase-coded  pulses  employing  m-sequences  are  a  particularly  poor  choice  of 
pulse  waveform. 

In  Figs.  3.4  and  3.5,  the  results  of  numerical  evaluations  of  Pj.  for  the 

sine  pulse  (as  a  function  of  ps)  and  the  phase-coded  pulse  for  N  =  31  (as  a 

function  of  are  shown  for  the  four  models  for  the  delay  power-density 

spectra.  For  comparison,  the  irreducible  error  probabilities  obtained  by 

approximating  the  parameter  as  are  also  shown.  By  examining  the  results 

in  Figs.  3.4  and  3.5  for  the  cases  when  the  true  value  of  Hy  is  included  in 

the  computation,  and  when  Hy  is  approximated  as  T2  (in  which  case  |i  = 

(2qv)  /T),  it  is  clear  that  this  parameter  is  a  significant  factor  only  for 

phase— coded  pulses.  In  fact,  the  error  introduced  by  the  approximation  for 

the  evaluation  of  the  irreducible  error  probability  for  the  sine— puls?  is 

quite  small.  However,  assuming  H  =  T2  for  the  evaluation  of  the  irreducible 

P 

error  probability  for  a  phase-coded  pulse  with  N  =  31  can  produce  a  result  in 
error  by  more  than  two  orders  of  magnitude! 

Notice  that  some  of  the  irreducible  error  probabilities  in  Figs.  3.4  and 
3.5  are  shown  for  different  ranges  of  normalized  rms  multipath  spread  for  the 
cases  where  is  approximated  as  T^.  This  was  necessary  to  ensure  that  the 
assumption  of  adjacent-pulse-limited  ISI  holds  but  is  also  indicative  of  the 
relative  values  of  qy  for  different  data-pulse  shapes.  That  is,  for  a  fixed 
channel  model  and  fixed  rms  delay  M,  there  is  substantial  variation  between 


the  values  of  qv  for  different  pulse  shapes.  The  roles  of  Hy  and  qy  in 
determining  the  normalized  rms  multipath  spread  is  made  apparent  by  examining 
the  relationship  between  these  parameters  for  different  pulse  shapes  with  a 
fixed  channel  model.  The  values  of  (2tiv)*^/T  for  the  Gaussian  delay  spectrum 
are  given  in  Table  3.3  for  the  rectangular,  sine,  and  phase-coded  pulse  shapes 
using  characteristic  m-sequences  of  lengths  N  =  7,  IS,  31,  and  63  given  in 
[32,  Fig.  7].  The  values  of  the  normalized  rms  multipath  spread  pv  are  listed 
in  parentheses  in  Table  3.3.  Notice  that  for  large  rms  delays,  the  parameter 
Up  for  the  phase-coded  pulse  can  be  much  smaller  than  the  corresponding 
parameter  for  the  rectangular  pulse.  However,  the  values  of  normalized  rms 
multipath  spread  for  the  phase-coded  pulse  can  be  significantly  larger  than 
the  corresponding  values  for  either  the  rectangular  or  sine  pulse.  Also, 
notice  that  both  and  (i seem  to  be  independent  of  the  sequence  length  N  for 
any  fixed  value  of  rms  delay.  This  is  a  bit  surprising  since  it  might  be 
expected  that  the  effects  of  intersymbol  interference  for  phase-coded  pulses 
employing  m-sequences  could  be  substantially  reduced,  for  any  rms  delay,  by 
simply  using  an  m-sequence  of  greater  length.  Unfortunately,  the  results  in 
Table  3.3  indicate  that  unless  the  sequences  are  carefully  selected, 
increasing  the  sequence  length  can  increase  the  sensitivity  to  delay  spread, 
and  that  the  phase-coded  pulse  is  at  least  as  vulnerable  to  the  effects  of 
intersymbol  interference  as  the  rectangular  pulse. 

The  above  results  for  the  phase-coded  pulse  are  a  consequence  of  both  the 
channel  model  we  employ,  (i.e.,  that  of  a  continuous  delay  power-density 
spectrum),  and  the  interaction  between  the  data  pulse  correlation  function 
&v({)  and  the  delay  spectrum  g({).  In  particular,  for  very  small  values  of 
rms  delay,  the  power  spectrum  of  the  delay  is  concentrated  close  to  the  origin 


Table  3.3.  (2ny)  /T  and  (normalized  rms  multipath  spread  |iy)  as  a  function 

of  the  ras  delay  H  for  the  rectangular  and  sine  pulses  and  for  the 
phase-coded  pulse  with  sequence  lengths  N  =  7.15,31,63 


M/T  rect 


N  -  7 


N  =  15  N  =  31 


N  =  63 


(0.01) 

(0.00001) 

(0.01067) 

(0.01139) 

(0.01297) 

(0.01639) 

0.05 

(0.05) 

0.0016 

(0.00158) 

0.0496 

(0.06874) 

0.0499 

(0.08980) 

0.0191 

(0.04728) 

0.0423 

(0.14683) 

0.1 

(0.1) 

0.0119 

(0.01235) 

0.0782 

(0.13992) 

0.0929 

(0.22867) 

0.0278 

(0.09631) 

0.0495 

(0.24218) 
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so  that  the  parameter  jip  depends  most  heavily  on  the  portion  of  Rp(£)  that 
corresponds  to  very  small  values  of  £.  From  (3.17b)  ve  see  that,  for  values 
of  $  in  the  range  [O.TJ.  |Rp(§)  |  =  lRr($)  |  =  $.  Unless  the  signature 
sequences  are  carefully  chosen,  this  relationship  can  hold  for  several 
consecutive  chip  intervals.  Thus,  it  is  clear  that  for  small  values  of  rms 
delay,  the  normalized  rms  multipath  spreads,  and  hence  the  irreducible  error 
probabilities  are  nearly  the  same  for  the  rectangular  and  phase-coded  pulse 
shapes. 

Finally,  the  results  in  Figs.  3.3-3.S  indicate  that  the  irreducible  error 
probability  strongly  depends  on  the  normalized  rms  multipath  spread  From 
these  results  and  the  results  in  Table  3.3,  one  might  erroneously  conclude 
that  the  performance  of  the  phase-coded  pulse  is  approximated  by  the 
performance  of  the  rectangular  pulse.  In  what  follows,  we  demonstrate  that 
pulse  shapes  that  produce  similar  irreduc ible  error  probabilities,  do  not 
necessarily  exhibit  similar  performance  for  practical  signal-to-noise  ratios. 
These  results  do  however  show  that  the  irreducible  error  probability,  as  a 
function  of  the  normalized  rms  multipath  spread,  is  insensitive  to  variations 
in  the  shape  of  the  delay  power-density  spectrum.  In  fact,  the  maximum 
variation  of  Pj  for  the  pulse  shapes  evaluated  in  Figs.  3. 3-3. 5  is  less  than  a 
factor  of  1.2  for  fixed  normalized  rms  multipath  spread.  This  latter 
observation  indicates  that  the  limiting  error  probability  for  a  particular 
frequency-selective  Rayleigh  fading  channel  can  be  minimized  by  simply 
choosing  a  data  pulse  waveform  which  minimizes  the  normalized  rms  multipath 
spread  (defined  as  in  (3 . 23 ) — ( 3 .25) )  for  the  delay  power-density  spectrum  of 
the  fading  channel. 
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3.2.2  Effects  of 


Interference :  Bounds 


The  relative  insensitivity  of  the  irreducible  error  probability  to  the 
shape  of  the  delay  spectrum  for  a  fixed  normalized  rms  multipath  spread 
suggests  the  possibility  of  obtaining  meaningful  bounds  on  the  probability  of 


error.  In  this  section  we  describe  a  method  of  obtaining  bounds  on  P^  for 
frequency-selective  Rayleigh  fading  channels. 

3. 2. 2.1  Rectanaular  pulse 

Under  the  assumption  of  adjacent-pulse-limited  ISI>  the  functions  defined 
by  (3.12)  can  be  written  (for  the  rectangular  pulse)  as 


ir  =  I*/2  -  2Tf (g)  +  p*/ 2  . 


(3.26a) 


£ 


nr  *  M2' 


n;  -  Tf(g)  -  n;/2 


(3.26b) 


(3.26c) 


with  the  normalized  rms  multipath  spread  4^  given  by  (3.23)  and  where  the 


functional  f(g)  is  given  by 


fT 

f(*>  *  J  $  *U>  d?  . 

w 


(3.27) 


Bence,  for  fixed  4^,  the  irreducible  error  probability  is  completely  specified 
up  to  the  determination  of  f(g).  Equations  (3.13)  and  (3.26)  imply  that  the 
irreducible  error  probability  is  a  continuous  function  of  f(g)  so  that  it  is 
possible  to  characterize  Pj  over  the  range  of  values  taken  on  by  the 
functional  f(g).  Thus,  upper  and  lower  bounds  on  Pj  can  be  obtained  if  the 
minimum  and  maximum  values  of  f(g)  can  be  found  for  a  fixed  normalized  rms 
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multipath  spread.  The  characterization  of  the  variation  of  f(g)  is 
accomplished  through  an  application  of  a  result  due  to  Dubins  [33]. 

Let  M  be  the  set  of  all  noir-negative  measures  m  on  the  Borel  sets  of 
[0,T]  such  that  m([0,T])  =  1/2.  It  can  be  shown  [33]  that  the  extreme  points 
of  N  (denoted  by  ex  M)  are  the  point  masses  on  [0,T]  and  that  M  is  compact  on 
the  weak*  topology  [341.  If  m  €  M,  the  hyperplane  H  defined  by 

f  x(t)  m(dt)  *  c 
J0 

for  constant  c  is  closed  and  bounded  for  x  €  C[0,T]  [35].  Dubins'  Theorem 
guarantees  that  every  extreme  point  of  M’  =  {mDH}  is  a  convex  combination  of 
at  most  two  extreme  points  of  M. 

In  terms  of  the  present  application,  upper  and  lower  bounds  on  Pj.  can  be 
obtained  for  each  value  of  normalized  rms  multipath  spread  by  finding  the 
minimum  and  maximum  values  of 


f(m)  =  f  5  B(d?)  , 
J0 


(3.28) 


for  m  €  M  subject  to  the  constraint 


r2-f  52  m(d§)  =  h2/2, 

J0 


(3.29) 
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where  m(d()  *  g(()d£.  Since  the  functional  f  in  (3.27)  is  a  continuous  linear 
functional  on  m.  it  obtains  its  minimum  value  on  an  extreme  point  of 
N*  =  HOH.  Hence,  by  Dubins'  Theorem,  the  minimum  of  f  is  found  by  searching 
over  all  convex  combinations  of  two  point  masses  (impulse  functions)  in  M'  for 
the  smallest  value  of  f(m);  i.e., 

fmin  "  minf  :  b  €  ex  M'  }  . 
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This  search  in  [0,T]2  is  easily  implemented.  In  fact,  for  the 

rectangular  pulse  it  can  be  shown  that  the  pair  (f  .  ,f  )  i*  eiven  bv 

min  m&x  °  J 

2 

(Rj /2,|ir/2}  for  pr  <  1/2.  The  resulting  bounds  for  the  irreducible  error 
probability  are  then  found  by  searching  over  the  interval  tfm^n«fma  ]  for  the 
minimum  and  maximum  values  of  Pj,  it  is  necessary  to  carry  out  this  final 
search  since  Pj  is  not  necessarily  a  monotonic  function  of  f(g).  In  Fig.  3.6, 
the  minimum  and  maximum  irreducible  error  probabilities  for  the  rectangular 
pulse  are  shown  as  a  function  of  pr.  Notice  that  for  a  fixed  value  of  pr,  the 
total  variation  from  the  minimum  to  the  maximum  value  of  Pj.  is  very  small;  the 
two  bounds  differ  by  a  factor  of  less  than  1.5. 

3. 2. 2. 2  Sine  pulse 

Using  the  fact  that  Hs  =  T2 ,  and  assuming  adjacent-pulse-limited  ISI,  we 
can  write  the  functions  in  (3.12)  for  the  sine  pulse  as 


A  =  T^f '  (g)  -  2Tf"(g)  +  h2/2  , 


and 


ns  =  *J2’ 


Ti;  =  Tf"(g)  -  n“/2 


(3.30a) 

(3.30b) 


(3.30c) 


with  the  normalized  rms  multipath  spread  ns  given  by  (3.24)  and  where  the 
functionals  f'(g)  and  f"(g)  are  given  by 

0 

and 


T 

f'(g)  «  f  cos2(n$/T)  g(5)  d? 

»  A 


(3.31a) 


*”<«)  “  f  [-?cos2(n«/T)  +  T/n  sinU5/T)cos(n$/T)]  g($)  d$  .  (3.31b) 

0 
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Tims,  we  may  proceed  as  in  the  previous  section  where  the  constraint  &  is  now 
given  by 

T 

t‘M  [-$cos(n$/T)  +  T/n  sin(n£/T)  ]2  &<$)  d$  =  n*/2.  (3.32) 

J0  s 

except  that  in  this  case  the  set  {f ’ (m) ,f”(m) ;  m  €  ex  M'J  is  a  region  in  . 

Hence,  the  minimum  and  maximum  obtainable  error  probabilities  are  found  by 

searching  over  the  convex  region 


Q  =  co  { (f * (m) , f"(m) )  :  m  €  ex  M'}  €  R2  , 


where  co  (S)  denotes  the  convex  hull  of  S.  However,  in  order  to  avoid  the 
task  of  completely  characterizing  Q,  we  perform  the  two-dimensional  search  for 
the  minimum  and  maximum  Pj  over  the  rectangular  region 


0 


1 


X  ff".  f"  11 
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which  contains  Q.  The  resulting  bounds  for  the  irreducible  error  probability 
for  the  sine-pulse  are  shown  in  Table  3.4  as  a  function  of  the  normalized  rms 
multipath  spread  ns.  Notice  that  even  though  we  have  "weakened"  the  bounds  on 
Pj  by  performing  the  search  over  Q  rather  than  Q  and  assuming  that  Hs  =  T^ , 
the  total  variation  between  upper  and  lower  bounds  is  still  no  larger  than  a 
factor  of  4  (for  fixed  p$). 

3. 2. 2. 3  Phase-coded  pulse 

Finally,  we  address  the  problem  of  obtaining  bounds  on  Pj  for  the  phase- 
coded  pulse,  which  is  the  most  complicated  of  the  three  pulse  shapes 
considered.  Under  the  assumption  of  adjacent-pulse-limited  ISI,  the  function 
ilp  can  be  written  as 


Table  3.4.  Irreducible 
respect  to 

error  probability  bounds  for 
tbe  normalized  rms  multipath 

sine  pulse  DPSK  with 
spread 

- 

a 

Lower  Bound 

Upper  Bound 

0.0005 

0.8361 *10~7 

1.2273.10-7 

a’ 

0.001 

3.7526*10 ~7 

5.2812.10-7 

..  a 

v  i* 

0.005 

1.0311*10~5 

1.4204-10-5 

0.01 

4.1277«10~5 

6.1089.10-5 

0.05 

1.0262.10-3 

2.4748-10"3 

0.1 

0.4104*10~2 

1.6259-10-2 

* 

w 
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Unfortunately,  there  is  no  obvious  way  to  characterize  the  functions  n  and  n' 

P  'P 


in  terns  of  the  normalized  rms  multipath  spread  and  a  few  linear 


functionals  of  g.  Moreover,  since  strongly  depends  on  the  shape  of  the 


delay  power-density  spectrum  and  can,  in  fact,  equal  zero  for  certain  examples 


of  multipath  delay  spectra,  there  is  not  a  clear  relationship  between  and 


(as  evidenced  by  the  data  in  Table  3.3).  Thus,  a  given  value  of  normalized 


rms  multipath  spread  does  not  necessarily  correspond  to  a  linear  constraint  on 
the  set  of  delay  power-density  spectra.  Alternatively,  if  we  assume  that 


*  1^  and  proceed  as  before,  the  resulting  upper  bound  for  the  irreducible 


error  probability  is  approximately  1/2. 


However,  if  we  make  two  additional  assumptions  concerning  the 
characterization  of  the  channel,  both  upper  and  lower  bounds  on  the 

performance  of  phase-coded  pulse  DPSK  can  be  obtained.  We  assume  that  a 

fraction  of  at  least  1/2N  of  the  total  power  of  the  delay  spectrum  lies  in  the 
range  I**Tc,T  ]  where  N  is  the  length  of  the  signature  sequence.  Second,  we 
assume  that  the  frequency-selective  channel  is  represented  by  a  specular 
multipath  channel  with  paths  corresponding  to  delays  in  the  set 

(0,+Tc,+2Tc, • • • ,+(N-l)Tc}  so  that  the  delay  power-density  spectrum  g(lj)  is 

represented  by  a  discrete  density  spectrum  g(kTc).  This  model  may,  in  a 
sense,  be  more  appropriate  for  phase-coded  pulse  since  systems  employing 
phase-coded  waveforms  can  resolve  the  multipath  components  of  the  delay 
spectrum  for  many  channels  of  interest.  We  proceed  by  finding  the  minimum  and 
maximum  values  of  and  under  the  constraint, 

m~2  n  _  . 

H  >  C,(N-k)*g(kT  )  =  yil/2. 


(3.34) 


and  performing  the  search  for  the  minimum  and  maximum  Pj  over  the  region 


Q  =  Ui i 


p.min’^p.max^  ^  ^p.min'^p.max^  ^  ^p.min'^p.max 


]} 


where  q  =  2c^T^/N.  The  resulting  bounds  for  the  irreducible  error 

pi  min 

probability  for  the  phase-coded  pulse  are  given  in  Tables  3.5  and  3.6, 
respectively,  for  m-sequences  of  length  N  =  7  and  N  =  31  discussed  above. 


Notice  that  there  is  almost  no  variation  between  upper  and  lower  bounds, 

and  that  the  bounds  for  the  phase-coded  pulses  are  very  similar  to  those 

obtained  for  the  rectangular  and  sine  pulses.  The  small  variation  between  the 

bounds  in  Tables  3.5  and  3.6  is  partially  a  result  of  the  first  assumption; 

there  is  a  path  with  non-zero  energy  corresponding  to  zero  delay.  However, 

the  "closeness"  of  these  bounds  is  primarily  due  to  the  fact  that  modeling  the 

delay  density  as  a  discrete  spectrum  allows  the  effects  of  the  parameter  to 

be  incorporated  in  the  bounding  procedure.  (Recall  that  this  parameter  is  a 

significant  factor  in  determining  the  value  of  the  normalized  ms  multipath 

spread.)  In  fact,  it  is  easy  to  demonstrate  that  if  we  assume  H  =  Tr  ,  the 

P 

resulting  upper  and  lower  bounds  differ  by  a  factor  of  N  (the  sequence  length) 
for  normalized  rms  multipath  spreads  of  about  0.01. 


Since  the  bounds  in  Tables  3.5  and  3.6  are  obtained  under  assumptions 
that  are  more  restrictive  than  those  used  to  obtain  the  bounds  for  the  other 
pulse  shapes,  it  might  be  argued  that  the  results  for  phase-coded  pulses  are 
of  little  consequence.  However,  it  seems  reasonable  to  assume  that  some 
portion  of  the  received  signal  energy  corresponds  to  zero  (or  near-zero) 
propagation  delays,  and  that  this  portion  is  at  least  as  large  as  the  average 
of  the  energy  received  in  any  interval  of  width  T  Since  the  upper  bound  on 
Pj  for  phase-coded  pulses  can  approach  1/2  if  there  is  no  signal  component  in 


B 

%  Table  3.5.  Irreducible  error  probability  bounds  for  N  =  7  phase-coded  pulse 

>  with  respect  to  the  normalized  rms  multipath  spread  p 


".A 


D 


"I 


v. 


V 


cu 


MP 

Lower  Bound 

Upper  Bound 

0.0005 

1.1176.10-7 

1.303 9. 10-7 

0.001 

4.9919.10-7 

5.1781.10-7 

0.005 

1.2543.10-5 

1.2565.10-5 

0.01 

4.9754.10-5 

5.0252.10"5 

0.05 

1 .1186  *10-3 

1.4069.10-3 

0.1 

3.3669.10~3 

7.3647-10"3 

0.5 

3.7202.10-2 

2.5337  «10-1 

Table  3.6.  Irreducible  error  probability  bounds  for  N  =  31  phase-coded  pulse 
with  respect  to  the  normalized  rms  multipath  spread  p 


P 


P 


Lower  Bound 


Upper  Bound 


0.0005 

0.001 

0.005 

0.01 

0.05 

0.1 


1.1176.10- 7 

4.9919.10- 7 
1.2267.10"5 

4.5614.10- 5 
0.3679.10-3 
0.4683 .10-3 
5.8824 •10-3 


1.3039.10 

5.1781.10 
1.2850.10' 

5.4812.10 
4.2587.10' 
4.8358.10' 


0.5 


4.3373*10 
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the  range  [-T  ] ,  this  assumption  can  be  viewed  as  a  necessary  condition  for 

c  c 

the  operation  of  phase-coded  poise  DPSK. 

3.2.3  Effects  of  Additive  Noise:  Characterization 

Consider  the  response  of  in-phase  and  quadrature  correlation  receivers 
matched  to  v(t)  (and  assumed  synchronized  in  time)  to  a  s ingle,  transmitted 
data  pulse  signal 

s(t)  =  b0v(t)  ;  0  <  t  <  T  , 


where  E  is  the  energy  per  data  bit.  In  the  absence  of  additive  noise,  the  sum 


of  the  squares  of  the  filter  outputs  is 


16o2E  r 


m  16o2E 


J  g(?)  rJ( UD  d? 


G(ft)  PV(Q)  da  , 


(3.35a) 


(3.35b) 


by  Rayleigh's  theorem,  where  Py(Q)  is  the  inverse  Fourier  transform  of 
Ry(l5l).  The  corresponding  output  due  to  the  additive  white  Gaussian  noise 
with  (two-sided)  spectral  intensity  NQ/2  is  8NqT. 

For  a  nonselective  Rayleigh  fading  channel,  g(£)  =  &(£)  so  that  the 
frequency  correlation  function  G(Q)  is  equal  to  unity  across  the  band  of 
interest.  In  this  case,  the  signal-to-noise  ratio  (snr)  at  the  output  of  the 
matched  filter  becomes 


snr  =  2a  f  (3.36) 
N0 

since,  by  definition,  Ry(0)  *  T2.  Notice  from  (3.35)  that  the  function  Py(Q) 
can  be  viewed  as  a  power-spectrum  of  the  single-pulse,  or  "one-shot"  matched 
filter  output  prior  to  sampling,  in  non-selective  fading  so  that  (3.35b)  and 


i 


the  resalting  signal-to-noise  ratio  is  maximized  when  the  frequency 
correlation  function  G(Q)  is  flat  across  the  band. 

In  contrast,  one  of  the  essential  characteristics  of  frequency-selective 
fading  is  the  fact  that  the  frequency  correlation  function  G(Q)  is  not  flat 
but  resembles  the  transfer  function  of  a  band-limited  data  channel.  We  point 
out  that  this  is  only  a  resemblance  and  not  an  equivalence.  Rather,  in  the 
absence  of  additive  noise,  the  function  G(Q)  is  a  measure  of  the  degree  of 
statistical  correlation  between  two  received  spectral  components  separated  by 
Q  Hz.  Thus,  in  the  case  of  frequency-selective  fading  (3.35b)  indicates,  in  a 
statistical  sense,  what  fraction  of  the  received  power  (due  to  a  single  pulse) 
is  available  at  the  output  of  the  matched  filter.  This  observation  is  made 
apparent  in  the  resulting  expression  for  the  single-pulse  signal-to-noise 
ratio  for  frequency-selective  fading,  given  by 

2or2E  2ff2E 

snr  =  - -  r  g(Q)  B  (Q)  dfl  =  - T  f  gU)  R2(UI)  d$ 

N  T2  v  u  -j*  la, 


2VT2 


(3.37) 


which  implies  that  qy  is  the  key  parameter  in  determining  the  signal-to-noise 
ratio  for  single-pulse  matched  filter  detection.  Notice  that  if  we  neglect 
the  effects  of  intersymbol  interference  on  the  error  probability  in 
frequency-selective  fading  so  that  both  qy  and  q^  (cf.  (3.12))  are  taken  to  be 
zero,  the  average  error  probability  (3.7)  becomes 


4°  E  2qv  ,-l 

N0  ’  1*  J 


(3.38) 


*  •  *  »  '  »  •  ’  •  *  *  •'»*».  * 
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The  results  in  Section  3.1  demonstrate  that  for  small  signal-to-noise 
ratios,  the  average  error  probability  for  frequency-selective  Rayleigh  fading 
is  well  approximated  by  »  (2  +  4<t^E/Nq)-*,  which  is  the  probability  of 
error  for  non-selective  fading.  While  this  approximation  is  valid  independent 
of  data-pulse  shape,  the  above  discussion  implies  that  (3.38)  is  a  much  better 
approximation  of  the  error  probability  for  small  signal-to-noise  ratios  than 
Pns.  In  fact,  for  a  given  pulse  shape,  the  asymptotic  error  probability  for 
small  signal-to-noise  ratios  is  given  by  (3.38).  Moreover,  for  all  examples 
of  data-pulse  shape  and  delay  density  considered  here,  (3.38)  may  be  taken  as 
a  lower  bound  for  the  average  probability  of  error  for  all  values  of  snr. 
Hence,  the  parameter  qy  in  some  sense  determines  the  position  of  the 
"inverse-linear"  lower  bound  on  the  error  probability  as  a  function  of 
signal-to-noise  ratio  S  *  2o2E/N(j.. 

The  differences  between  the  signal-to-noise  ratios  in  (3.36)  and  (3.38) 
are  probably  not  a  major  concern  for  signals  with  relatively  small  time- 
bandwidth  products,  viz.  the  rectangular  and  sine  pulses.  However,  for  large 
time-bandwidth  product  pulses,  such  as  the  phase-coded  pulse  with  moderately 
large  N,  this  factor  becomes  significant  and  plays  a  central  role  in  the 
choice  of  pulse  shape.  For  a  channel  with  Gaussian  delay  power-density 
spectrum,  the  values  of  the  parameter  qy/T^  are  given  in  Table  3.7  for  the 
rectangular  and  sine  pulse  shapes,  and  for  the  phase-coded  pulse  shapes  using 
the  m-sequences  of  length  N  *  7,  15,  31,  and  63  discussed  above.  The  results 
in  Table  3.7  show  that  for  small  to  moderate  values  of  M/T,  a  large  percentage 
of  the  received  signal  power  is  utilized  by  the  matched  filter  detector  for 
both  the  rectangular  and  sine  pulses  for  this  channel.  If  M/T  =  0.1,  for 
example,  the  matched  filter  output  for  the  rectangular  and  the  sine  pulse 


52 


i 

r.y, 

& 


r-' 

:#■ 


N*: 

i 


I'-’v 


fe* 


& 

£ 

>  A 

* 

sg 


represents  about  85%  and  91%  of  the  total  received  signal  power,  respectively. 
In  comparison,  only  about  9%  and  5%  of  the  received  signal  power  is  utilized 
for  matched  filter  detection  of  the  phase-coded  pulse  for  N  -  31  and  N  =  63. 

In  contrast  to  the  analogous  results  for  the  normalized  rms  multipath 
spread,  the  parameters  for  phase-coded  pulses  uniformly  exhibit  a  strong 
dependence  on  the  sequence  length.  It  can  be  shown  that,  for  a  large  class  of 
channel  models  and  typical  sequence  lengths  N,  the  parameter  is  closely 
approximated  by  k/N  where  k  depends  on  the  channel.  For  the  above  example, 
this  constant  is  approximately  2.62,  1.35,  and  0.33  for  M/T  equal  to  0.05, 
0.1,  and  0.5,  respectively.  The  low  signal- ^o-noise  ratio  asymptotic  error 
probability  (3.38)  becomes 

,~1 


(3.39) 


for  the  phase-coded  pulse  with  sequence  length  N. 

In  the  previous  section,  it  was  demonstrated  that  the  phase-coded  pulse 
waveform  is  at  least  as  susceptible  to  the  effects  of  intersymbol  interference 
as  the  rectangular  pulse.  Thus,  phase— coded  pulses  provide  little  potential 
for  improvement  in  the  limiting  error  probability.  Moreover,  when  coupled 
with  the  data  in  Table  3.7,  these  results  indicate  that  the  average  error 
probabilities  for  the  phase-coded  pulse  are  relatively  large  for  all  practical 
values  of  S  *  2o2E/N£>.  For  example,  from  Table  3.7  for  M/T  =  0.05  and  N  =  63, 
we  see  that  for  fixed  error  probability,  the  additional  signal-to-noise  ratio 
required  for  the  phase-coded  pulse  relative  to  the  rectangular  pulse  is  about 
0.4614/0.0424  or  10.37  dB. 

The  relationships  between  irreducible  error  probabilities,  available 
signal  power,  and  system  error  performance  can  be  seen  by  examining  Fig.  3.7. 
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The  average  probability  of  error  is  given  as  a  function  of  S  =  2<j^E/Nq  for  the 
rectangular,  sine,  and  phase-coded  pulse  (N  =  31)  for  the  Gaussian  delay 
power-density  spectrum  with  rms  delay  M  =  T/10.  The  low  signal-to-noise  ratio 
asymptote  given  by  (3.38)  and  the  irreducible  error  probabilities  (which 
appear  as  horizontal  asymptotes)  are  shown  for  each  pulse.  Notice  that  while 
the  limiting  error  probability  for  the  phase-coded  pulse  with  N  =  31  is 
roughly  equivalent  to  that  of  the  rectangular  pulse,  the  actual  error 
probability  for  the  phase-coded  pulse  is  higher  than  for  the  rectangular  pulse 
for  values  of  signalrto-noise  ratio  less  than  35  dfi.  Also,  notice  that  the 
sine  pulse  exhibits  the  best  performance  of  the  four  pulse  shapes  considered 
for  the  entire  range  of  signal-to-noise  ratios.  These  results  are  in 
agreement  with  the  data  presented  in  Tables  3.3  and  3.7  which  show  that  the 
normalized  rms  multipath  spread  ng  for  the  sine  pulse  is  significantly  smaller 
than  the  corresponding  parameter  for  the  rectangular  pulse.  Moreover,  the 
parameter  which  is  related  to  the  utilizable  signal  power,  is  largest  for 
the  sine  pulse  while  the  corresponding  parameter  for  the  phase-coded  pulse 
indicates  severe  performance  degradation  in  the  presence  of  additive  noise. 

Finally,  we  point  out  that  the  performance  of  phase-coded  pulse  DPSK  in 

Rayleigh  fading  is  related  to  previous  results  concerning  the  performance  of 

coherent  spread-spectrum  communications.  In  particular,  the  results  in  [26] 

and  [31]  demonstrate  that  the  effects  of  interference  caused  by  multiple 

propagation  paths  for  relative  delays  in  the  range  [-T+T  ,t-T  ]  can  be 

substantially  reduced  by  the  choice  of  the  signature  sequence.  More 

importantly,  it  is  shown  that  the  spread-spectrum  receiver  (which  is 

essentially  a  coherent  matched  filter)  severely  discriminates  against  signals 

with  propagation  delays  outside  the  range  [-T  ,j  ]. 

c  c 


This  property  is 
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generally  desirable  for  communications  in  a  multipath  environment  when  there 
is  at  least  one  relatively  strong  path.  However,  it  is  precisely  this  property 
that  degrades  the  error  probability  of  spread-spectrum  communications  in  a 
frequency-selective  Rayleigh  fading  channel. 

3.3  Performance  Approximations  for  Other  DPSK  Systems 
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In  the  previous  section  it  was  demonstrated  that  bounds  on  the  error 
probability  of  DPSK  in  a  WSSUS  frequency-selective  Rayleigh  fading  environment 
can  be  obtained  when  the  channel  delay  spectrum  is  normalized  with  respect  to 
the  autocorrelation  function  of  the  specified  data-pulse  shape.  The  resulting 
bounds  show  that  the  limiting  error  probability  is  relatively  insensitive  to 
variations  of  the  channel  delay  power-density  spectrum  (for  fixed  normalized 
rms  multipath  spread). 

The  bounds  on  system  performance  represent  useful  results  since  the 
limits  on  the  effects  of  frequency-selectivity  on  DPSK  communications  have  not 
been  previously  characterized.  Moreover,  the  closeness  of  these  bounds  imply 
that  it  is  not  necessary  to  fully  characterize  the  channel  in  order  to  obtain 
an  accurate  estimate  of  system  performance.  This  is  of  considerable 
importance,  since  in  practice  it  is  difficult  to  precisely  describe  the  fading 
process  [13] . 
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Another  implication  of  the  above  results  is  that  it  may  be  possible  to 
estimate  the  average  probability  of  error  of  a  proposed  system  from  the 
tabulated  error  probabilities  in  Section  3.2.  In  particular.  Figs.  3.3  and 
3.4  indicate  that  the  irreducible  error  probability  for  the  rectangular  pulse 
as  a  function  of  the  normalized  rms  multipath  spread  closely  approximates 
the  irreducible  error  probability  for  the  sine-pulse  as  a  function  of 


Similarly,  Figs.  3.3  and  3.4  indicate  that  the  irreducible  error  probability 
for  the  rectangular  pulse  as  a  function  of  pr  closely  approximates  the 
irreducible  error  probability  for  the  phase-coded  pulse  as  a  function  of  p^. 
In  Fig.  3.8,  the  irreducible  error  probabilities  for  the  rectangular  pulse, 
the  sine  pulse,  and  the  phase-coded  pulse  (N  =  31)  are  shown  as  a  function  of 
their  respective  normalized  multipath  spreads,  pv.  The  results  in  Figs.  3.3- 
3.5  indicate  that  there  is  little  to  be  gained  by  considering  more  than  one 
delay  power-density  spectrum.  Hence,  for  this  example  we  have  represented  the 
fading  channel  by  a  Gaussian  delay  power-density  spectrum. 

The  results  in  Section  3.2  and  Fig.  3.8  indicate  that  there  is 
essentially  a  single  relationship  (represented  by  Fig.  3.8)  between  the 
irreducible  error  probability  Pj.  and  the  normalized  rms  multipath  spread  pv 
for  DPSK  in  frequency-selective  Rayleigh  fading.  We  make  the  following 
observation  concerning  this  result:  The  data  in  Fig.  3.8  in  some  sense 
represents  the  irreducible  error  probability  as  a  function  of  p^  for  a  generic 
(time-limited)  data-pulse  waveform  v(t).  Hence,  the  limiting  error 
probability  for  a  particular  Rayleigh  fading  channel  and  a  particular  time- 
limited  data-pulse  waveform  u(t)  (say)  can  be  closely  approximated  by  means  of 
a  single  rms-type  [13]  channel  parameter, 

\  -  t  iu*J^  gU>R2<5>  d?}1/2  -  l  §  •£  g(5>R2(UI)  dU1/2 

*  t  5  *f  g($)R2(UI)  d{}1/2  (3.40) 

Q  “« 

with  the  autocorrelation  function  R2($)  defined  by  (3.14).  Of  course,  the 
last  equality  in  (3.40)  follows  only  if  the  adjacent-pulse-limited  ISI 
assumption  holds.  (However,  the  data  presented  in  Section  3.2  indicates  that 
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this  mast  be  the  case  if  the  performance  is  to  be  acceptable  at  all.)  The 
resalting  approximation  for  the  irreducible  error  probability  is  found  as  the 
corresponding  valne  of  Pj  in  Fig.  3.8. 

For  small  time-bandwidth  prodnct  pulse  shapes,  a  good  estimate  of  the 
average  error  probability  as  a  function  of  the  signal-to-noise  ratio  S  is 
realized  as  a  curve  which  lies  above  Pns  in  (3.18)  and  Pj  (calculated  from 
(3.40))  and  which  asymptotically  approaches  Pqs  and  Pj  for  extreme  values  of 
the  signal-to-noise  ratio.  If  it  is  also  possible  to  measure  the  quantity 

CD  0D 

2\l^  “  l« I)  d?  -  J  G(Q)  PU(Q)  dQ  ,  (3.41) 

then  a  more  accurate  approximation  (especially  for  large  time-bandwidth  pulse 
shapes)  can  be  found  as  above  by  replacing  P  by  the  low  signal— to— noise 
ratio  asymptote  given  in  (3.38)  for  the  parameter  iiu  calculated  from  (3.41). 

There  may,  of  course,  exist  design  situations  for  which  an  approximation 
of  error  .  performance  is  unacceptable  for  a  system  evaluation.  In  these 


situations,  it 

may  be  necessary  to  construct 

and 

test  the  proposed  system 

or 

obtain  a  more 

comple  te 

characterization 

of 

the  channel 

.  However, 

the 

approximations 

discussed 

above  could  be 

of 

significant 

utility  in 

the 
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consideration  of  DPSE  systems  employing  complicated  pulse  shapes  and  in 
applications  to  fading  channels  which  are  difficult  to  fully  characterize. 
For  example,  if  phase-coded  pulse  waveforms  are  under  consideration  as 
candidates  for  the  pulse  shape,  these  approximations  may  be  useful  in  the 
selection  of  signature  sequences,  as  well  as  in  the  comparative  system 
evaluations  which  would  provide  insight  into  the  potential  benefits  of  hybrid 
frequency— hopping/direct— sequence  systems  versus  simple  frequency— hopped 
systems  employing  DPSK  modulation. 
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In  the  analysis  of  the  preceding  sections,  we  have  assumed  that  the 
response  of  the  fading  channel  to  a  deterministic  input  signal  s(t)  is 
represented  by  a  zero-mean  complex  Gaussian  random  variable  (i.e.,  Rayleigh 
fading).  There  are,  however,  situations  where  the  fading  channel  is  more 
accurately  modeled  as  Rician.  This  channel  model  arises  when  the  received 
signal  contains  a  nonfaded  specular  component  as  well  as  a  Rayleigh  faded 
component. 

The  narrowband  model  for  the  received  signal  in  the  case  of  Rician  fading 
is  given  by  (2.3)  with  the  parameter  a  >  0  representing  the  strength  of  the 
specular  component.  As  in  the  case  of  Rayleigh  fading,  the  decision  statistic 
in  (3.5)  can  be  written  as  a  quadratic  form  of  Gaussian  random  variables  U 
and  V.  The  resulting  conditional  probability  of  error,  given  a  particular 
sequence  of  data  bits  (b,^)  is  transmitted,  is  related  to  the  probability 
distribution  of  the  ratio  of  Chi-square  variables,  each  having  two  degrees  of 
freedom;  the  non-centrality  parameters  depend  on  the  characteristics  of  the 
system.  Expressions  for  the  conditional  error  probabilities  are  derived  in 
(19]  and  [36],  where  in  each  case  the  result  is  obtained  through  the  inversion 
of  the  characteristic  function  of  Hermitian  quadratic  form  of  complex  Gaussian 
random  variables  [37]  by  means  of  a  double  integral  evaluated  by  Price  [38]. 

Results  of  previous  investigations  of  the  performance  of  digital 
communications  over  Rician  fading  channels  indicate  that  the  error  probability 
is  upper  and  lower  bounded  by  the  corresponding  results  for  the  Rayleigh 
channel  (no  specular  component)  and  the  additive  white  Gaussian  noise  channel 
(uo  fading),  respectively.  In  the  remainder  of  this  section,  we  show  that  the 


average  probability  of  error  for  DPSK  in  Rician  fading  is  completely 
determined  by  tbe  parameter  a  in  (2.3)  and  the  parameters  necessary  for  the 
evaluation  of  the  error  probability  for  Rayleigh  fading  considered  in  the 
previous  sections.  Thus,  the  results  for  Rayleigh  fading,  which  are 
essentially  the  results  for  the  "worst-case"  Rician  channel,  are  also  useful 
in  determining  the  error  performance  in  the  more  general  Rician  case. 

It  can  be  shown  that  the  conditional  error  probabilities  for  the  Rician 
channel  depend  on  the  characterization  of  the  Rayleigh  faded  portion  of  the 
received  signal  only  through  the  parameters  rjy,  r\y,  and  Hy  defined  in  (3.12) 
and  discussed  in  Section  3.1.  In  particular,  assuming  without  loss  of 
generality  that  b^_^  *  b^  for  the  transmitted  signal  s(t)  given  by  (3.1),  the 
conditional  error  probability  [19]  can  be  written  as 


Kkt)  -  .VmS -  [  ^(-‘>‘vY(*i>  =  1  • 

2  »  “uu^— i^“W^— i^ 


exp 


[- 


A(bi)+B(fei) 


]  *  V  >'A(&i)B(&i))  ,  (3.42) 


where  Q(.,*)  is  Marcum's  Q  function  [6],  and  IQ(x)  is  the  modified  Bessel 
function  of  the  first  kind  of  order  zero  and  argument  x.  The  parameters  A  (Ik) 
and  B(J^)  in  (3.42),  which  are  related  to  the  non-centrality  parameters  of  the 
corresponding  Chi-square  variates,  are  given  by 


A(&.)  - 

4a  ET 

[®mj<li)+®vv*-i*  -2*bibi-l*  y ®OU^ii^®W^i^  ] 

[  (3.43a) 

1 

■uU<bi)BVy(fei) 

and 

B(fc.)  - 

4a2ET 

[BOU^i)+®yy(bi)  +2(bibi_1)  ^ »uo^i^n,W^i^  ] 

.  (3.43b) 
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The  moments,  “^(ii^  (3.43)  are  the  central  moments  of  the  random  variables 

U  and  V  defined  as 

^(bj)  =  E{  (X  -  Edl^D-fY  "  EtYlb.])*!^)  . 

Using  (2.3)  and  (3.6)  for  the  signal  defined  by  (3.1),  we  see  that  E{U}  = 
o  \Z8ET.b.  and  E{VJ  =  av/lET-b.^  so  that  the  moments  in  (3.43)  are  precisely 
the  same  as  those  defined  in  (3.10)  for  the  analysis  of  Rayleigh  fading 
channels. 

It  is  easy  to  see  that  (3.42)  reduces  to  the  error  probability 
expressions  for  Rayleigh  fading  (see  (3.8)  and  (3.9))  when  a  =  0.  In  fact,  if 
a  =  0,  then  A(b^)  =  B(b^)  =  0,  and  (3.42)  becomes 


p<5>i> 


</ mup^i^ntvy(— i^  *  “uv^-i* 

2  y  "^(^“w^i) 


(3.44) 


since  Q(0,0)  =  Iq(0)  *  1,  which  is  identical  to  (3.8)  for  the  case  b^^  =  b^. 
Alternatively,  if  there  is  no  faded  component  in  the  received  signal,  then 
2a'“  =  0  (cf.  (2.4))  so  that  the  moments  in  (3.43)  become 


and 


®uu  =  “w  =  8N0T 
muv  =  °- 


II  =  b ^ ,  then  A(b^)  =  0,  and  B(]^)  =  2a  E/Nq  and  (3.42)  becomes 


P(b 


No 

(since  Q(0,x)  3  exp(-x/2)  [6]),  which  is  the  average  error  probability  of  DPSK 


fox  the  additive  white  Gaussian  noise  channel. 
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CHAPTER  4 

PERFORMANCE  OF  BINARY  FSK  COMMON I CATIONS  OVER 

FREQUENCY-SELECTIVE  FADING  CHANNELS 

In  this  chapter,  we  consider  the  evaluation  of  the  average  probability  of 
error  for  binary  FSK  communications  over  WSSUS  fading  channels.  As  in  Chapter 
3,  the  analysis  is  first  carried  out  for  the  case  of  Rayleigh  fading  channels 
(i.e.,  with  the  channel  response  given  by  (2.3)  for  a  =  0) .  It  is  then  shown 
that  the  error  probability  for  the  more  general  Rician  channel  is  completely 

specified  in  terms  of  the  key  system  parameters  used  in  the  performance 

evaluation  for  Rayleigh  fading. 

Previous  results  on  the  performance  of  FSK  in  frequency-selective  fading 
include  those  found  in  [2,17,18].  In  [17],  Bello  and  Nelin  consider  the 
performance  of  FSK  where  the  frequency  separation  between  the  two  signals  is 
assumed  to  be  an  integer  multiple  of  1/2T.  The  channel  is  modeled  as  a 
Rayleigh  fading  channel  with  a  Gaussian  delay  power-density  spectrum.  The 
degree  of  frequency-selectivity  is  defined  as  the  distance  between  the  1/e 
points  of  the  frequency  correlation  function  G(Q).  It  is  shown  that  the 

system  error  performance  depends  on  the  relative  phases  as  well  as  the 

frequency  separation  of  the  two  transmitted  signals.  The  authors  claim  that 
FSK  is  generally  less  sensitive  to  the  effects  of  frequency-selective  fading 
than  rectangular  pulse  DPSK.  The  authors  also  state  that  for  large  signal- 
to-noise  ratios,  FSK  produces  lower  error  rates  than  rectangular  pulse  DPSK 
for  channels  which  are  highly  selective.  However,  these  claims  are  not 
supported  by  the  numerical  results  for  the  irreducible  error  probability  in 
[16]  which  differ  from  the  results  in  [17]  by  an  order  of  magnitude. 
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In  [2]  and  [18],  the  performance  of  FH/FSK  for  Rician  fading  channels  is 
considered  for  a  triangular  delay  power-density  spectrum  model.  However,  in 
both  cases  the  authors  assume  that  the  two  FSK  signals  are  orthogonal  at.  the 

receiver.  In  Section  4.2,  we  show  that  this  is  equivalent  to  assuming  that 

the  frequency  separation  between  the  FSK  tones  is  infinite. 

In  this  chapter,  the  effect  of  frequency-selective  fading  on  the 

performance  of  noncoherently-detected  FSK  is  investigated  for  a  number  of 
signaling  formats.  The  parameters  that  characterize  the  effects  of 
intersymbol  interference  and  additive  channel  noise  are  identified.  It  is 
shown  that  the  performance  of  FSK  communications  over  frequency-selective 
Rayleigh  fading  channels  can  be  well  approximated  in  terms  of  these  key 
parameters  which  are  related  to  rms-type  channel  measurements.  The  results 
for  FSK  are  compared  to  the  results  for  DPSK  in  Chapter  3.  It  is  found  that 

the  error  performance  of  these  two  systems  depends  on  a  number  of  common 

factors  and  can  be  quite  similar  in  some  cases.  The  technique  for  obtaining 
bounds  on  the  irreducible  error  probability  for  DPSK  is  applied  to  the 

analysis  of  FSK.  Finally,  we  consider  the  approximation  of  the  performance  of 
other  FSK  systems  in  terms  of  the  key  system  parameters. 

4.1  System  Model 

The  information  source  for  the  FSK  transmitter  is  modeled  as  a  sequence 
of  mutually  independent  random  variables,  each  taking  on  values  in  the 

set  (0,1}  with  equal  probability.  The  binary  data  sequence  (b^  is  given  by 

the  mapping 


bi  =  (-1) b i 


from  {0,1}  to  {-1,1}  for  each  integer  i.  Using  the  narrowband  signal  models 
discussed  in  Chapter  2,  the  transmitted  signal  s(t)  is  defined  by 

s ( t )  =  ^  v(t-iX)e"jtbinht/T  +  0(i,bi)]  (4.1) 

for  t  €  [iT,(i+l)T3  for  each  integer  i.  The  real-valued  waveform  v(t)  is 
assumed  to  be  time  limited  to  the  interval  [0,T3  such  that 


so  that  E  is  the  energy  per  data  bit. 

The  parameter  h  in  (4.1)  is  known  as  the  deviation  ratio  or  modulation 
index  of  the  FSK  signal  set  {39}  and  is  related  to  the  frequency  deviation  f^ 
from  the  carrier  frequency  f^  (which,  of  course  does  not  appear  in  the 
narrowband  signal  representation)  by  h  *  2fdT.  The  two  transmitted  signals 
are  referred  to  as  the  mark  signal  and  the  space  signal .  We  follow  the 
convention  in  {39}  by  denoting  the  signal  transmitted  at  frequency  f  -  f .  as 

C  Q 

the  mark  signal  and  the  signal  transmitted  at  frequency  f£  +  fd  as  the  space 
signal.  From  (4.1)  we  see  that  the  mark  signal  corresponds  to  the  case  when 
=  -1  (i.e.,  bj  *  1),  and  the  space  signal  is  transmitted  when  b^  =  1 
(i.e.,  E^=o).  The  phase  angle  0(i,b^),  which  depends  on  both  the  value  of 
the  i-th  data  bit  and  the  interval  in  which  it  is  transmitted,  represents  the 
phase  of  the  transmitted  signal  s(t)  at  time  t  -  0. 

Depending  on  the  method  used  to  generate  the  transmitted  signal,  the 
phase  angles  0(i,-l)  and  0(i,+l),  corresponding  to  the  mark  and  space  signals, 
respectively,  may  be  related  or  may  be  modeled  as  statistically  independent 
random  variables.  For  example,  if  the  transmitted  signal  is  obtained  by 
switching  between  two  uncoupled  oscillators,  the  two  phase  angles  are  modeled 


65 


as  independent  random  variables  uniformly  distributed  on  [0,2n]. 
Alternatively,  if  the  transmitted  signal  is  a  continuous-phase  FSK  (CPFSK) 
signal,  then  the  phase  angles  are  related  by  0(i,+l)  =  0(i+l,-l)  +  2irhi.  We 
assume  in  any  case  that  no  abrupt  phase  transition  occurs  between  two 
consecutive  data-pulse  intervals  if  the  two  corresponding  data  bits  are 
identical.  That  is,  we  assume  that  0(i,+l)  =  0(i+l,+l)  and  that  0( i ,— 1)  = 

0<i+l,-l). 

By  using  narrowband  signal  models  to  represent  the  signals  in  the  FSK 
system,  we  have  removed  the  dependence  of  the  orthogonality  of  the  two 
transmitted  signals  on  the  carrier  frequency  f  Rather,  the  two  signals  are 
said  to  be  orthogonal  if  and  only  if 

fT  y2(t)e-j[2nht/T  +  0<O,+1)-0(O,-1)]  dt  .  fT  y2(t)e-j[2nht/T]  dt  .  0  # 

Jo  •'o 

Unless  stated  otherwise,  "orthogonality"  is  used  to  describe  the  relationship 

of  the  mark  and  space  signals  at  the  transmitter.  Hence,  for  narrowband 

signal  models,  the  orthogonality  (or  lack  of  orthogonality)  between  the  mark 

and  space  signals  is  completely  determined  by  the  pulse  waveform  v(t)  and  the 

modulation  index  h.  This  is  equivalent  to  the  assumption  that  both  f  and  h 

are  integer  valued,  or  that  (f  T  -  h)  >>  1. 

c 

The  FSK  receiver  is  modeled  as  the  noncoherent  matched  filter  receiver 
shown  in  Fig.  4.1,  which  is  discussed  in  [39].  During  each  time  interval, 
[iT, (i+l)T],  the  receiver  forms  the  decision  statistic  Z represented  by 


i  *  4 1  J(  1)T  r(t)v*(t)eJ,Tt/T  dt  | 2  -  4 1  JU+1)T  r(t)v*(t)e_jnt/T  dt  | 


(i+l)T 


lul2  -  |V|2  . 


(4.2) 


Figure  4.1.  Noncoherent  matched  filter  receiver 
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where 


t  -  2  fU+1,T  r(t)»*(t)eJ'"^1  dt  . 
JiT 

V  =  2  f(l+1)T  r(t)v*(t)e-jnth/T  dt  . 

J  •  T* 


(4.3a) 


(4.3b) 


From  (4.2)  we  see  that  the  decision  statistic  for  the  i-th  data  bit  is  a 
quadratic  form  of  complex  Gaussian  random  variables  D  and  V. 

Under  the  assumption  that  p(t-x.g)  »  0;  |g|  >  T,  (i.e.,  adjacent-pulse- 
limited  ISI) ,  the  output  statistic  depends  on  at  most  three  consecutive 
data  bits  represented  by  b^  =  (bj_y,  b^,  b^+j).  The  probability  of  error  can 
be  written  as  the  average  of  probabilities,  each  conditioned  on  the  event  that 
one  of  the  eight  possible  sequences  of  data  bits  are  transmitted,  i.e.. 


F.  ■  i  i  P<ii>  • 

b . 

—  l 


(4.4) 


where  P(b^)  *  Pr  (error  occurslb^  transmitted}.  In  [11,17],  Bello  and  Nelin 
show  that  the  conditional  error  probabilities  are  given  by 


P(bi)  =  (2  +  y(b.))' 


where  the  "equivalent"  signal-to-noise  ratio  y(»)  may  be  written  as 


(4.5) 


Y(ii) 


2l(nup(bi)  ~  mvv(bi))  I _ 

f (®u|j(}>  ^)  +  ®yy(]l  j) )  —  4  l®jjy(b  ^)  I  —  I  (®u^j(b  ^)  -  myy(b  ^) )  I 


(4.6) 


a  function  of  the  moments  ®Xy(fei)  =  E{XY*lbi)  of  the  random  variables  U  and  V. 


I 


Using  (2.5)  and  (4.3),  the  moments  in  (4.6)  are  found  to  be 


,f(i+l)T  r(i+l)T  -T  „  . 

*Wii>  =  8oZJ  j  f  p(t-i,4)8(t-?)v *(t)eJnth/T 

J  iT  J  iX  iT 

•  s*(x~5) v(x)e“Jnxll/,T  dgdxdt  +  8Nqt,  (4.7a) 
,f(i+l)T  f(i+l)T  fT  .  . 

■W<li>  -  Sc M  I  I  p  ( t-x .  § )  s  ( t-5 )  v  ( t )  e~ J T 

J  iT  J  ij  ^ 

•s*(x-$)v(x)ej"th/T  d£dxdt  +  8Nqt,  (4.7b) 

and 


_  f(  i+l)T 

f(i+l)T  , 

8a  | 

I  f 

JiT 

JiT  - 

•s*(x-5)v(x)ejBxh/T  d£dxdt  +  8N0  J  v(t)v*(t)e-»2nth/T  dt.  (4.7c) 


By  substituting  g(«)  for  p(»,»)  in  (4.7),  we  see  that  the  moments  needed 
to  evaluate  y(bj)  and  hence  Pe  can  be  written  as 


2  T 

■Wii)  "  L  2lRv(h,bi.4)  I2  +  IVh.bi-i.5)!2  +  lRv(h.bi+1.^)|2 


+  2Re{e~jteT(i'^i)_(brbi-l)nh5/TJ.Rv(h,bi.5)-R*(h,bi_1,5) 


+  2Re{eJIeT(l+1'^i)+(bi+l'1)nh]-Rv(h,bi,4)-R*(h.bi+1.4)  ]  d$ 


“w^)  -  l6.—  j  «($)•[  2lRvOi.-bi.§)l2  +lR*(h,-bi_1,$)l2  +lR*(h,-bi+1,5)|: 


+  2Re{e"j[eT(i'^i)'(bi'bi-l),th5/T]-R*(h,-bi,5)-Rv(h.-b 
+  2Re{ejl0T(i+1'ii)+(bi+l+1),rhl.R*(h,-bi,5).Rv(ll,-bi+1,5)  ]  d$ 


+  8n0T,  (4.8b) 


and 


•V^)  =  26_  jT  g($)-[  j2ni?/T)-Rv(b.bi,?)-Rv(h.-bi,5) 


T  "0 


+  Rv(h'bi_1.5)*Rv(b.-bi_1.«)  +  e-j2nb(?_T)/T-Rv(h.bi+1.5)*Rv(b.-bi+1,5) 
+  ej[0T(i'bi)“(brbi-l)nb4/T]-Rv(h.bi_1,§)>Rv(b.-bi.5) 

+  e-j[eT(i+l.l.p+2nh5/T+nh(bi+1-l)].Rv(hbi+i^).iv(h#_bif0 
+  eJIeT(i+1^i)"2,Th5/T+nb(bi-l+1)].Rv(h,bi,5)-Rv(b.-bi+1.5)  ]  d? 


+  8Nq  J  v(t)v*<t)ej2ntb/T  dt  .  (4.8c) 


The  phase  transitions,  ^(i.b^),  which  denote  the  differences  between  the 
phases  of  the  mark  and  space  signals  at  time  t  =  0  is  given  by 


eTu.b.)  = 


e<i,b.)  -  eu-l.bj^)  ;  bi  t  b^ 


(4.9) 


;  bi  =  bi-l  * 


The  correlation  functions  R  (h.b^)  and  R^h.bi.S)  in  (4.8)  are  defined 
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EyOi.bi.S)  =  J  v(t)  v(t-5>  ejnht(1_bi,/T  dt 


(4.10a) 


RyCh.bj,?)  =  f?  v(t)  v(t+T-5)  ejnht(1-bi)/T  dt  (4.10b) 

Jq 

for  0  <.  £  .<  T.  Notice  that  Rv(h,b^,£)  and  Ry(h,b^,£)  are  related  to  the 
aperiodic  autocorrelation  functions  defined  by  (3.10a)  and  (3.10b).  In  fact, 
using  (3.10a)  and  (4.10a),  it  is  easy  to  see  that  Ry(h,l,5)  =  Ry(lj). 

Similarly,  from  (3.10b)  and  (4.10b)  we  have  that  Ry(h,l,£)  =  Ry($).  In 
contrast,  the  functions  Rv(h.-1,£)  and  Ry(h,-1,£)  are,  in  general,  complex¬ 
valued  functions  that  depend  on  the  modulation  index  h  as  well  as  the  waveform 
v(t).  Unfortunately,  it  is  difficult  to  express  the  result  of  these  integrals 
in  compact  form  for  arbitrary  modulation  index.  Finally,  notice  that  the 

function  Rv(h.-1,0)  is  simply  a  measure  of  the  orthogonality  of  the  two  FSK 
signals. 


Using  (4.6),  (4.8),  and  the  fact  that  we  consider  only  binary  symmetric 
signal  sets  (cf.  (4.1)).  it  can  be  shown  that  r^)  =  yi-bp.  This  property, 
along  with  (4.4),  implies  that  P#  in  (4.5)  can  be  reduced  to 

pe  =  4  IP(-1.+1.-1)  +  P(+l,+l,-l)  +  P(-l,+l,+l)  +  P(+1,+1,+1)J,  (4.11) 
so  that  we  need  only  consider  the  case  when  a  space  signal  is  transmitted 
during  the  i-th  data-bit  interval. 


In  the  present  analysis,  we  consider  a  number  of  examples  of  binary  FSK 
signal  formats.  Vithin  the  class  of  CPFSK  signals,  we  consider  modulation 
indices  h  *  1/2,  1,  2,  and  10  for  a  rectangular  waveform  v(t).  The  example 
of  rectangular  pul: s  FSK  for  h  *  10  is  considered  as  an  approximation  to  the 
case  of  "infinite"  frequency  separation  between  the  two  FSK  tones.  Notice  that 
the  CPFSK  signal  with  h  ■  1/2  produces  what  is  known  as  a  minimum-shift-keyed 


Wfsv, 
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or  MSK  signal  [39].  Of  coarse,  it  is  not  common  practice  to  consider  the 
performance  of  MSK  signals  for  noncoherent  detection  since  MSK  signals  do  not 
comprise  an  orthogonal  signal  set.  By  examining  R^fh.-l.O)  in  (4.10a),  it  is 
easy  to  see  that  h  =  1  is  the  smallest  modulation  index  producing  an 
orthogonal  signal  set  for  the  rectangular  pulse.  We  include  the  case  of  h  = 
1/2  as  an  example  in  order  to  gain  a  better  understanding  of  the  role  of  the 
modulation  index  in  determining  the  performance  of  noncoherent  FSE  in  the 
presence  of  frequency-selective  fading.  Since  we  assume  that  knowledge  of  the 
phase  of  the  transmitted  signal  is  unavailable  to  the  receiver,  the  phase 
transitions  for  BISK  signals  are  related  by  ©jd.bj)  =  0  and  Gjd+l.b^  =  n  or 
equivalently,  0j(i,l^)  =  it  and  O-jd+l.b^)  =  0  for  differing  consecutive  data 
bits. 

Other  examples  of  FSE  signals  considered  include  those  generated  by 
phase-coded  waveforms  defined  by  (3.4).  The  signature  sequences  for  this 
pulse  waveform  are  the  m-sequences  discussed  in  Chapter  3.  Notice  that  for 
the  phase-coded  pulse,  the  orthogonality  of  the  transmitted  signals  depends 
only  on  the  modulation  index  h  so  long  as  the  pulse  waveform  is  composed  of  a 
sequence  of  rectangular  pulses.  In  this  case,  Rp(h,-1,0)  =  Rr(h,-1,0),  where 
the  subscripts  "p”  and  "r"  denote  the  phase-coded  pulse  and  the  rectangular 
pulse,  respectively.  Hence,  h  =  1  is  the  smallest  modulation  index  for  which 
the  phase-coded  pulse  signals  are  orthogonal.  However,  for  fixed  modulation 
index,  the  spectral  spreading  effect  of  the  direct-sequence  waveform  causes 
significantly  more  overlap  of  the  spectra  of  the  mark  and  space  signals  than 
for  the  rectangular  pulse.  For  a  given  sequence  length  N,  we  can  specify  a 
set  of  orthogonal  phase-coded  pulse  signals  with  roughly  the  same  degree  of 
spectral  overlap  as  rectangular  pulse  FSE  with  h  =  1  by  let’  *  h  =  N  for 


phase-coded  pulse  signals.  In  what  follows,  we  consider  examples  of  phase- 
coded  pulse  FSK  where  the  modulation  index. is  taken  as  the  sequence  length, 
i.e.,  h  =  N,  as  well  as  examples  with  h  =  1. 

In  Fig.  4.2,  the  results  of  the  numerical  evaluation  of  the  average  error 
probability  of  CPr5K  with  modulation  index  h  =  1  are  shown  as  a  function  of 
the  signal-to— noise  ratio  S  =  2o^E/Nq  for  a  Gaussian  delay  spectrum  with  rms 
multipath  spread  p  =  0.05.  The  probability  of  error  for  orthogonal  FSK  in 
non-select ive  Rayleigh  fading,  given  by 

Pns  “  <2  +  2°2e/n0>-1  '  (4.12) 

is  also  shown.  From  Fig.  4.2  we  see  that  for  small  signal-to-noise  ratios, 
the  average  error  probability  is  largely  determined  by  the  value  of  P  As 


the  signal-to-noise 

ratio 

becomes  large, 

the 

effects  of 

intersymbol 

interference  become 

more 

significant  and 

the 

average  error 

probability 

asymptotically  approaches 

the  irreducible 

error 

probability 

for  FSK  in 

frequency-selective 

fading 

defined  by  Pj  =  Pg 

Is  -> 

w.  The  results 

in  Fig.  4.2 

imply  that  Pj  (shown  as  a  horizontal  asymptote)  is  a  key  performance  parameter 
for  the  design  and  evaluation  of  practical  systems. 


The  irreducible  error  probabilities  Pj  for  several  examples  of  FSK  signal 
sets  and  the  four  examples  of  delay  spectra  given  in  Chapter  2  are  shown  in 
Figs.  4.3  and  4.4  as  a  function  of  the  rms  multipath  spread  p  =  M/T.  In  all 
cases,  the  results  are  obtained  using  standard  numerical  integration 
techniques.  In  Fig.  4.3,  the  results  are  shown  for  rectangular  pulse  CPFSK 
with  modulation  indices  h  «  1/2  and  1,  and  for  rectangular  pulse  FSK  with  h  = 
1  and  G^.(i,— 1)  *  0y(i+l,— 1)  *  0-j>  -  «,  where  0j  denotes  the  common  phase 


transition. 


Notice  that  for  modulation  index  h  =  1,  the  irreducible  error 


probabilities  exhibit  significant  sensitivity  to  the  shape  of  the  delay 


density.  However,  the  variations  caused  by  changes  in  the  shape  of  the  delay 
density  are  small  compared  to  the  variations  between  the  examples  for  0^,  =  o 
and  0^,  =  n.  The  results  in  Fig.  4.3  for  the  two  examples  of  CPFSK  (i.e.,  h  = 
1/2  and  h  =  1)  seem  to  indicate  that  the  irreducible  error  probability  is  also 
sensitive  to  small  changes  in  the  modulation  index.  However,  in  Section  4.2 
it  is  shown  that  the  larger  values  of  Pj  for  CPFSK  with  h  =  1/2  (MSK)  are 
primarily  due  to  orthogonality  considerations  rather  than  the  effects  of 
intersymbol  interference. 

In  Fig.  4.4,  the  irreducible  error  probabilities  for  rectangular  pulse 
CPFSK  with  h  =  10  and  phase-coded  pulse  FSK  with  h  =  1  and  signature  sequence 
length  N  =  31  with  0T  =  0  are  shown.  For  this  example  of  rectangular  pulse 
CPFSK.  we  see  that  the  irreducible  error  probability  is  not  sensitive  to 
changes  in  the  shape  of  the  delay  density  for  large  values  of  rms  multipath 
spread.  Also  notice  that  the  irreducible  error  probability  for  phase-coded 
FSK  exhibits  the  same  erratic  behavior  as  the  results  for  phase-coded  DPSK. 

4.2  Performance  Bounds  and  Approximations 

In  this  section,  the  relationships  between  irreducible  error  probability, 
intersymbol  interference,  and  signal  orthogonality  are  explored  by  first 
characterizing  the  frequency-selective  fading  mechanism  in  the  delay  domain. 
The  results  for  FSK  are  compared  to  the  analogous  results  for  DPSK  in  Chapter 
3.  We  show  that,  in  some  cases,  the  limiting  error  performance  of  FSK  is 
dependent  on  parameters  that  are  closely  related  to  the  key  parameters  for 
DPSK.  We  then  focus  on  the  frequency-domain  characterization  of  the  fading 
channel  and  identify  the  system  parameters  which  provide  an  indication  of 


average  error  probability  with  respect  to  the  degree  of  frequency-selectivity 
and  the  degradation  due  to  additive  noise. 


4.2.1  Effects  of  Intersvmbol  Interference :  Characterization 

Consider  the  effects  of  the  interference  from  the  signal  corresponding  to 
on  the  detection  process  for  the  signal  corresponding  to  bQ  =  +1.  If  we 
consider  the  portions  of  the  received  signal  corresponding  to  positive 
propagation  delays  only,  the  decision  statistic  in  (4.2)  is  given  by  ZQ  =  |u I 2 

A 

-  I V I  with  the  random  variables  D  and  V  given  by 

jv„. 

[v(t-5)e“j,lh(t“5)/T  +  v(  t+T-5)  e~J  ^b-lffb^  t~$)  /T  +  0(-l  .b^)]  jd?dt  (4>13a) 

and 

[v(t-5)e"jnh(t”^/T  +  v(t+T-§)e"JIb-lnh(t~5)/T  +  »b-i>  J  jd^dt.  (4.13b) 

Using  (2.4),  (2.3)  and  (4.7),  the  moments  b^j  and  oyy  are  found  to  be 
2_  X 

“OD  *  ill—  J  *(?)*[  I*v(h,l,?)|2  +  lRy(h,b_1,?)  I2 

+  2Re(e'jleT(0’io)"<1_b-l)n^/T1.Ev(h.1.?)*R*(h.b_1.?)}  ]  d$  (4.14a) 

and 

"W  -  J*  g(?)*[  lR*(h,-l,?)|2  +  lR*(h,-b_1,5)|2 

+  2Re(e"j[eT(0'i0)-(1-b-l),Tb^/T].R*(i1,-l,5).Rv(h,-b_1,5)}  ]  d«-  (4.14b) 
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la  view  of  the  complicated  form  of  (4.14a)  and  (4.14b),  it  is  useful  to 
consider  the  evaluation  of  these  moments  for  specific  examples  of  signaling 
formats.  In  particular,  for  the  case  when  the  modulation  index  h  is  integer- 

valued  and  v(t)  is  the  rectangular  pulse  given  by  (3.2).  the  autocorrelation 
functions  in  (4.10)  are  given  simply  as 


T  -  5 


;  b  =  l 


R(h.b,5)  = 


(4.15a) 


-§eJ  nh^/T, sinc( h£/T)  ;  b  =  -1 


;  b  =  1 


R(h,b,5)  =  J 


(4.15b) 


$*jnh?/T*sinc(W;/T)  ;  b  -  -1  . 


where  smc(x)  sm(nx) /(nx) .  The  bracketed  terms  appearing  in  (4.14)  depend 

only  on  the  transition  phase  9^0.^)  and  the  functions  in  (4.15).  For  this 
example,  the  terms  appearing  in  the  expression  for  are  given  by 


if  b_!  =  1,  (4.16a) 


(T-5)2  +  £2  sine2 ( h^/T) 


+  2<T-5)5cos(9T(o,b0)-nh4/T)sinc(h5/T),  if  b_1  =  _  (4>16b) 

Similarly,  the  bracketed  terms  for  ^  in  (4.14b)  become 

°.  ..  . 

lf  b_!  =  1  (4.17a) 


i  r  ««*Va>V.aA  A/’  r 


$2  +  52sinc2(h$/T)  -  25co8(©T(o,b0)-jrh{/T)§*inc(h?/T),  if  b..^  -1.  (4.17b) 


Under  the  assumption  that  the  delay  power-density  spectrum  g($)  is  symmetric, 
it  is  easy  to  see  that  similar  expressions  result  from  the  consideration  of 
the  signal  components  corresponding  to  negative  delays  and  the  effects  of  the 
signal  corresponding  to  bi+^. 

By  examining  (4.16)  and  the  expression  for  Zi  in  (4.2),  we  can  interpret 
as  the  average  "power"  in  the  desired  signal  component  of  the  received 
signal  when  a  space  signal  is  transmitted.  Similarly,  m^y  represents  the 
average  power  in  the  interfering  (nndesired)  component  of  the  received  signal. 
(The  roles  of  m^  and  myy  are  simply  reversed  if  the  desired  signal  happens  to 
be  a  mark  signal.)  Notice  that  (4.17a)  indicates  that  no  intersymbol 
interference  is  produced  by  the  successive  transmission  of  (b_^,bg)  if 
b_^  “  bQ,  while  if  b_^  -  -bQ,  the  intersymbol  interference  produced  by  the 
data  pulse  preceding  bQ  i*  given  by  the  integral  of  (4.17b). 

Now  consider  the  output  of  the  matched  filters  for  both  positive  and 
negative  delays  for  the  transmitted  sequence  bj  »  (+l,+l,+l)  in  the  absence  of 
additive  noise.  For  the  rectangular  pulse  waveform  and  integer-valued 
modulation  index,  we  see  from  (4.7),  (4.16),  and  (4.17)  that  the  moments  m^, 
■yy,  and  m^jy  are  given  by 

■n^+l.+l'+l)  -  16q2g  fT  g(?M  2T2 ]  d?  -  320^,  (4.18a) 

T  J0 

*yy(+l.+l,+l)  -  0.  (4.18b) 

and 


■Dy(+1.+1,+1)  -  0  . 


(4.18c) 
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Hence,  the  total  average  power  at  the  output  of  the  space  signal  filter  in  the 
absence  of  data  modulation  is  given  by  (4.18a)  as  32o2ET.  Using  (4.18a)  and 
(4.17b),  we  see  that  the  portion  of  detector  output  due  to  interfering  signals 
relative  to  the  total  output  is 

°VV(-1.+1,-1)  =  T~2 *  fT  g($)  [  S2  +  42 sine2 (hJj/T) 

mpjj(+lf+l,+l)  0 

-?[cos(eT(0,b0)-jih5/T)+cos(eT(l,b1)-jrh§/T)]5sinc(h5/T)  ]  d$  .  (4.19) 

From  the  discussion  above,  we  conclude  that  the  effects  of  intersymbol 
interference  are  well  characterized  by 

»»v(h.e)  =  /— myy.(:1,t1,~1.)  ,  (4.20) 

V  ■UD(+i.+i.+i) 

which  we  take  as  the  definition  of  the  normalized  rms  multipath  spread  for  FSK 
communications  in  frequency-selective  fading.  From  (4.19)  it  is  clear  that 
Hv(h,8)  in  (4.20)  depends  on  the  phase  transitions  9j(i,l>£)  as  well  as  the 
modulation  index  h. 

If  it  is  assumed  that  the  delay  power-density  spectrum  g(£)  is  a  unimodal 
function,  (4.19)  provides  a  means  to  predict  the  role  of  the  modulation  index 
and  the  phase  transitions  in  determining  the  effects  of  intersymbol 
interference.  A  positive,  real-valued  function  f(x)  is  said  to  be  unimodal  if 
f  is  non-increasing  for  increasing  I x I .  For  example,  the  four  models  of  delay 
spectra  discussed  in  Chapter  2  are  unimodal  functions.  If  g(£)  is  unimodal 
and  h  is  integer  valued,  it  is  easy  to  show,  using  (4.19),  that  the 
sensitivity  to  the  effects  of  intersymbol  interference  is  minimized  if  the 
common  phase  transition  0T  «  0  (i.e. ,  CPFSK) .  Similarly,  (4.19)  indicates 
that  the  effects  of  intersymbol  interference  are  maximized  when  9^,  =  n>  which 
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is  the  Halting  case  of  phase-discontinnous  FSK.  These  results  are  in 
agreement  with  the  irreducible  error  probabilities  shown  in  Fig.  4.3  for  the 
rectangular  pulse  with  h  *  1.  For  this  example,  the  assumptions  0^.  =  0  and 
0j  *  n  produce  limiting  error  probabilities  that  can  differ  by  more  than  an 
order  of  magnitude.  From  this  discussion,  it  is  also  clear  that,  for  unimodal 
delay  densities  and  integer- valued  modulation  indices,  the  irreducible  error 
probability  for  random  phases  must  lie  between  the  values  of  Pj.  for  0^=0  and 
*  n. 

As  the  modulation  index  increases,  (4.19)  indicates  that  the  normalized 
rms  multipath  spread  becomes  less  dependent  on  the  phase  transition  0^..  la 
fact,  it  is  easy  to  show  that  for  increasing  h,  the  normalized  rms  multipath 
spread  given  by  (4.20)  approaches 

r  -lfT  i  i1/2 

Pv(».0)  *  [  2T  *J  g({)  d?  ]  -  M/T,  (4.21) 

where  M  is  the  rms  delay  defined  in  (2.3).  We  conclude  that  for  large 
modulation  index  and  fixed  rms  delay,  the  irreducible  error  probability  as  a 
function  of  N/T  becomes  independent  of  the  shape  of  the  delay  densitj  as  well 
as  the  phase  transition  0^.. 

Now  suppose  that  the  modulation  index  is  some  fixed  integer.  If  g(()  is 
unimodal,  then  for  large  values  of  rms  delay,  the  value  of  |iy(h,0)  is 
dominated  by  2IH^  *  2T^pv(«»,0) ,  which  corresponds  to  the  integral  of  the  ^ 
term  in  (4.19).  This  may  be  seen  by  noting  that  the  integral  of  the 
sinusoidal  terms  in  (4.19)  becomes  small  for  increasingly  wide  delay  spectra. 
Alternatively,  as  the  rms  delay  (and  hence  the  first  term  in  (4.19))  becomes 
small,  the  sinusoidal  terms  become  significant  so  that  (iy(h,0)  is  influenced 
by  the  shape  of  the  delay  density  as  well  as  the  phase  transition  0_.  The 


irreducible  error  probabilities  in  Fig.  4.4  for  h  =  10  provide  a  good  example 
of  the  behavior  of  systems  with  large  modulation  indices.  Notice  that  Pj  for 
this  example  depends  on  the  shape  of  the  delay  spectrum  only  for  rms  delays 
less  than  0.01. 

We  next  examine  the  characteristics  of  the  parameters  discussed  above  for 
an  example  of  non-orthogonal  signals.  In  particular,  for  rectangular  pulse 
CPFSK  with  h  =  1/2  (i.e.,  MSK) ,  m^p(+l ,+i ,+i)  =  32o2ET  as  before.  However, 
the  normalized  rms  multipath  spread  is  given  by 

H2(h.*)  =  T"2J^  g(5)  (  s2+[cos.(n|/2T)  ]2-2$cos (n/2+7t$/2T) [" ]  )  d^.  (4.22) 

By  comparing  (4.22)  and  (4.19),  we  see  that  if  the  delay  density  is  unimodal, 
MSK  is  considerably  more  sensitive  to  the  effects  of  intersymbol  interference 
than  signals  with  integer-valued  modulation  index.  It  can  be  shown  that  this 
increased  sensitivity  is  partly  due  to  the  lack  of  orthogonality  between  the 
two  MSK  signals.  Notice  that,  in  the  absence  of  additive  noise,  (4.8b) 
becomes 

«“W  =  1—.E  fT  g(?)  (  [~Sn/2T2Tj  f*  52sinc2U/2T))  d£  =  32o2ET.[— 2]  (4.23) 

TO  jt 

^or  llj  =  (+1,+1,+1).  Equation  (4.23)  implies  that  if  three  consecutive  space 

signals  are  transmitted,  the  output  of  the  mark  signal  filter  is  more  than 

half  that  of  the  space  filter,  a  property  which  holds  even  for  nonfading 

channels.  In  comparison,  recall  from  (4.18b)  that  for  integer-valued 

modulation  index,  n»yy(+l ,+l ,+l)  =  0.  This  result  and  the  data  presented  in 

Fig.  3.3  suggest  that  MSK  signals  are  of  little  interest  for  noncoherent 

detection,  especially  for  applications  to  frequency-selective  fading  channels. 
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In  Table  4.1,  numerical  values  of  the  normalized  rms  multipath  spread 
defined  by  (4.20)  are  shown  as  a  function  of  the  rms  delay  M  for  the  Gaussian 
delay  power-density  spectrum.  The  results  are  listed  for  signals  with 
modulation  indices  h  =1,2.  and  10,  for  0T  =  0  (CPFSK)  and  0T  =  n.  Notice 
that  in  the  case  of  CPFSK,  there  is  substantial  variation  between  the 
normalized  rm.  multipath  spreads  for  different  modulation  indices.  In  fact, 
for  very  small  values  of  rms  delay,  this  variation  is  nearly  equal  to  h.  Also 
notice  that  in  all  cases,  the  difference  between  the  values  of  normalized  rms 
multipath  spread  for  0^  =  0  and  0^.  =  n  are  very  large.  It  is  also  interesting 
to  note  that  for  the  case  when  *  n,  there  is  almost  no  variation  for 
different  modulation  indices.  Numerical  evaluations  of  the  normalized  rms 
multipath  spread  for  MSK  signals  revealed  that  this  parameter  ranges  between 
0.64  for  M/T  =  0.001  and  0.82  for  M/T  =  0.5,  indicating  that  the  parameter 

^(h.0^)  is  not  a  useful  measure  of  the  effects  of  intersymbol  interference 
for  this  case. 

For  the  consideration  of  phase-coded  pulses,  it  is  much  more  difficult  to 
derive  analytical  expressions  which  provide  insight  to  the  character  of  the 
normalized  rms  multipath  spread  ny(h,0).  Alteratively,  consider  the  moment 
®UU  in  (4.7)  for  the  phase-coded  pulse  with  modulation  index  h  =  N,  where  N  is 
the  sequence  length.  If  we  neglect  the  terms  corresponding  to  the  additive 

■jj|j(+l,+l ,+l)  is  given  by 
2  T 

“ou  -  ifLf  Jo  g(?)*2[Rp(h.l.5)+Rj(h,1.5)+2ip(h,1.5)Rp(h,1,?)]  d£.  (4.24) 

Recall  that  Rp(h,l,15)  and  Rp(h,l,$)  are  given  by  Ip($)  and  Rp({). 
respectively,  where  these  latter  functions  are  the  autocorrelation  functions 
given  by  (3.10)  for  the  analysis  of  0PSK.  Hence,  in  the  absence  of  data 


Table  4.1.  Normalized  rms  multipath  spread  nv(h,0^.)  as  a  function  of  the  rms 
delay  M  for  rectangular  pulse  FSK  with  9™  =  0  and  (0T  =  n) 


7.4081.10"6 

1.4814.10*? 

7.4049*10  3 

(2.3915.10  3) 

(2.3914.10  3) 

(2.3896*10  3) 

1.3602.10”4 

2.7192.10"4 
(9.9 938- 10“3) 

1.3419.10**3 

(9.9984.10-3) 

(9.8481. 10**3) 

5.4384.10-4 

1 .0859*1 0-3 

5.1151.10-3 

(1.9988-10-2) 

(1.9950.10-2) 

(1.8839.10**2) 

1.0001.10"2 

2.5776.10-2 

5.4470.10-2 

(9.8481 .10-2) 

(9.4197.1 0  2) 

(5.5126.10-2) 

5.0002 .10~2 
(1.8839-10-1) 

8 ,6507 • 10-2 
(1.6123.10-1) 
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(1.0250.10-1) 

5.0000.10-1 

4.5710-10-1 

4 . 4052  * 10-1 

(4.8907.10-1) 

(4.5099*10  1) 

(4.4026 .10_1) 
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modulation,  the  average  power  at  the  output  of  the  space  signal  filter  is 
given  by  16<x^EHp/T,  where  the  parameter  Hp  is  defined  by  the  relation  (3.22). 
The  values  of  H^/J^  for  an  m-sequence  of  length  N  =  31  are  listed  in  Table  3.1 
for  a  Gaussian  delay  density  as  a  function  of  the  rms  delay  M. 

In  Table  4.2,  the  values  of  the  pp(h,0)  with  h  =  N  are  shown  for  phase- 
coded  pulses  with  sequence  lengths  N  ■  7  and  31.  Also  listed  in  Table  4.2  are 
the  values  of  pp(h,0)  for  modulation  index  h  *  1  (shown  in  parentheses).  In 
each  case  the  results  are  obtained  for  ©T  -  0.  Surprisingly,  these  results 
indicate  that  the  normalized  rms  multipath  spreads  are  generally  smaller  for 
h  «  1  than  for  h  -  N.  From  the  similarity  of  the  results  for  h  =  N  and  h  =  1, 
we  conclude  that  the  sensitivity  to  the  effects  of  intersymbol  interference  is 
la-rgely  determined  by  the  autocorrelation  properties  of  the  phase-coded  pulse 
waveform  rather  than  the  value  of  the  modulation  index.  From  this  we  might 
also  assume  that  the  normalized  rms  multipath  spread  is  relatively  independent 
of  the  phase  transition  0T.  The  evaluation  of  p  (h,^)  with  ©T  *  n  for  the 
cases  listed  in  Table  4.2  shows  that  pp(h,n)  differs  from  the  values  listed  in 
Table  4.2  by  at  most  a  factor  of  two.  Also  notice  that  there  is  no  clear 
dependence  of  pv(h,0j)  on  the  sequence  length,  N.  Finally,  we  point  out  that 
in  comparison  to  the  corresponding  results  for  phase-coded  DPSK  listed  in 
Table  3.3,  it  is  clear  that  phase-coded  pulse  FSK  is,  in  general,  much  more 
sensitive  to  the  effects  of  intersymbol  interference. 

The  discussion  leading  to  the  expression  for  m^  in  (4.18c)  also  provides 
an  interesting  interpretation  of  the  role  of  this  parameter  in  determining  the 
performance  of  FSK.  Recall  from  Chapter  2  that  the  rms  delay  M  is  related  to 
the  bandwidth  of  the  frequency  correlation  function  G(Q) .  For  small  rms 
delays  (large  frequency-correlation  bandwidths) ,  there  is  significant 
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Table  4.2.  Normalized  rms  multipath  spread  pv(h,0)  as  a  function  of  the  rms 

delay  11  for  phase-coded  pulse  FSK  for  N  =  7.31  with  h  =  N  and 
(h  =  1) 
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correlation  between  the  fading  processes  acting  on  the  two  signals.  The 

parameter  that  characterizes  the  relationship  between  the  effect  of  frequency 
correlation  and  the  orthogonality  of  the  transmitted  signals  is  m^(+i ,+i ,+i) , 
(or  more  precisely,  a»^y(+l  ,+l  ,+l)  /mjjjj(+l  ,+l  ,+l) ) .  For  the  rectangular  pulse 
waveform  and  integer— valued  modulation  index,  we  see  from  (4.18c)  that 

s  0.  In  contrast,  it  is  easy  to  see  from  (4.7c),  that  for  MSK 

signals. 


muv/muu 


-2  T  2jT^ 

«DV(+l,+l,+l)/mUTJ(+l.+l.+l)  -  T_  j  g(?)*[  “ ~ 

2  0 


(4.25) 


which  not  only  shows  that  the  received  signals  are  highly  correlated  but  also 
shows  that  the  value  of  the  correlation  is  independent  of  both  g  and  M. 
Notice  that  the  square  of  the  magnitude  of  (4.25)  corresponds  to  the  result  in 
(4.24),  confirming  that  the  lack  of  orthogonality  of  MSK  signals  is  a 
significant  factor  in  determining  resulting  error  probability  for  both 
selective  and  nonselective  fading  channels.  Finally,  we  point  out  that  the 
close  relationship  between  the  values  of  aipy(+l ,+l ,+l)  and  myy(+l ,+l ,+l)  for 
the  rectangular  pulse  waveform  does  not  hold  for  phase-coded  pulse  waveforms. 

In  fact,  for  N  *  31  and  Gaussian  delay  density  with  M/T  =  0.001,  numerical 
evaluations  of  these  parameters  reveal  that  m^  *  m^y/io  and  m^y  =  myy/100  for  h 
=  31  and  h  *  1,  respectively. 

The  discussion  leading  to  the  definition  of  the  normalized  rms  multipath 
spread  in  (4.20)  and  the  results  in  Tables  4.1  and  4.2  suggest  that  the 
normalized  rms  multipath  spread  is  a  good  indicator  of  the  effects  of 
intersymbol  interference  for  both  rectangular  pulse  and  phase-coded  pulse  with 
integer-valued  modulation  index.  Since  the  irreducible  error  probability  is 
the  performance  parameter  most  dependent  on  the  effects  of  intersymbol 


1  .aaas.iii  m..*  >n^.aiu 
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interference,  we  expect  that  if  the  delay  spread  is  normalized  with  respect  to 

Pv(h,0j),  Pj  should  exhibit  minimal  sensitivity  to  either  the  signal  format  or 
the  shape  of  the  delay  power-density  spectrum  g(£).  Numerical  evaluations  of 

Pj  as  a  function  of  the  normalized  rms  multipath  spread  show  that  of  all 
examples  of  pulse  waveform,  modulation  index,  and  delay  density  considered, 
the  cases  exhibiting  the  most  variation  are  rectangular  pulse  waveforms  with  h 
=  1  for  0^,  =  o  and  0^.  =  it.  The  irreducible  error  probabilities  for  these  two 
cases  are  shown  in  Fig.  4.5  as  a  function  of  Pr(h,©T)  for  the  Gaussian, 
exponential,  and  triangular  delay  densities.  Notice  that  even  for  these 
extreme  cases  the  maximum  variation  for  fixed  normalized  rms  multipath  spread 
is  less  than  a  factor  of  1.4.  These  results  show  that  the  limiting  error 
probability  is  almost  completely  determined  by  the  normalized  rms  multipath 
spread  Pv(h.0T)  defined  in  (4.20). 

Before  concluding  the  discussion  on  the  characterization  of  the  effects 
of  intersymbol  interference,  it  is  interesting  to  draw  comparisons  between  the 
normalized  rms  multipath  spread  parameters  for  DPSK  and  FSK.  In  particular, 
notice  from  (3.23)  and  (4.21),  that  for  large  h,  the  parameter  Hv(h,0)  =  pr 
which  is  the  normalized  rms  multipath  spread  for  rectangular  pulse  DPSK. 
Thus,  for  fixed  rms  delay  M  and  sufficiently  large  modulation  index,  the 
irreducible  error  probability  for  rectangular  pulse  DPSK  is  a  good 
approximation  for  Pj  for  FSK. 


Next,  notice  that  the  parameter  for  sine  pulse  DPSK  is  given  by 


g(5)  [52cos2(n?/T)  +  252sinc2U/T)  -  2£2cos(nS/T)  sinc^/T)  ]  d£  .  (4.26) 


For  the  case  when  the  phase  transition  0T  =  o,  (4.20)  becomes  (except  for  a 


constant  factor) , 
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rT 

J  g($)  [«2  +  2?2$inc2(5/T)  -  2?2cos(«5/T)*inc(5/T)]  d*.  (4.27) 

J0  L 

Hence,  ve  see  that  sine-pulse  DSPE  is  considerably  less  sensitive  to 
intersymbol  interference  effects  than  FSE  for  any  integer— valued  modulation 
index,  even  for  the  minimizing  case  of  phase  transition,  0^,  =  o. 


4.2.2  Effects  of 


l :  Bounds 


As  in  the  consideration  of  DPSE  in  Chapter  3,  the  insensitivity  of  the 
irreducible  error  probability  for  FSE  as  a  function  of  the  normalized  rms 
multipath  spread  suggests  the  possibility  of  obtaining  useful  bounds  on  the 
probability  of  error  with  respect  to  the  parameters  of  the  fading  channel.  In 
this  section,  the  method  of  obtaining  bounds  on  the  irreducible  error 
probability  for  DPSE  communications  in  a  frequency-selective  Rayleigh  fading 
environment  is  applied  to  the  evaluation  of  Pj  for  FSE. 

It  is  clear  from  the  definition  of  ny(h,0T)  in  (4.20)  and  the  discussion 
above  that  the  normalized  rms  multipath  spread  is  dependent  on  the  phases 
®j(i,-l)  and  0j(i+l,-l).  However,  the  results  in  the  previous  section 
indicate  that  0T(i,-l)  -  0T(i+l,-l)  *  0T  =  0  and  0T(i,-l)  =  0T(i+l,-l)  =  @T  = 
n  in  some  sense  represent  limiting  cases  in  terms  of  the  normalized  rms 
multipath  spread.  In  what  follows,  the  bounds  for  the  irreducible  error 
probability  are  considered  for  these  cases  only.  The  resulting  bounds  for 
these  cases  are  representative  of,  if  not  limits  for,  the  results  for  other 
choices  of  phases. 

4. 2. 2.1  Rectangular  pulse 

Under  the  assumption  of  adjacent-pulse-limited  ISI,  the  moments  defined 
in  (4.7)  can  be  written  (for  integer-valued  modulation  index)  as 
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■^(ti)  -  16o2E/T.| 

[  2T2  +  (Itb^j— 1}  +  I(bi+1— l})f'(g)]  . 

(4.28a) 

■yy^)  -  16o2E/T.| 

[  (1{bi-lfc11  +  ICbi+1=-l})*T2»iJ(h,0T)]  . 

(4.28b) 

and 

■gyi]^)  -  16o2E/T. 

[  Ilbi-l““ 1{bi+1“- l)fw(g)]  * 

(4.28c) 

where  the  normalized  rat  multipath  spread  jiy(h,0T)  is  given  by  (4.20),  and 
I{.}  is  the  indicator  function.  The  functionals  f'(g)  and  f"(g)  in  (4.28)  are 
given  by 

f'(g)  -  f  g($)  ft2  -  2«T  +  $2sinc2($/T) 

J0 

+  2(T-5)?2cos(ifC/T+0T)*inc(5/T)]  d{,  (4.29) 

and 

T 

f"(g)  -  J  g(?)  [2(?-T)«eJ(,,,l«/T,.sinc(5/T) 

+  cos(eT)(T-«)52*»inc(«/T)  -  ?2cos(8T) ,inc2(«/T) ]  d{  .  (4.30) 

Thus  we  may  proceed  as  in  Chapter  3  by  characterizing  the  functions  f’(g)  and 
f*(g)  for  delay  speotra  g(g)  satisfying  the  constraint,  given  by 
T 

T2  g($)  [$2  +  {2sinc2(«/T)  -  2?2cos(9rn?/T)sinc(5/T)] 

“  (4.31) 

Notice  that  t"  is,  in  general,  a  complex-valued  function  of  g.  From  (4.S)  and 


(4.6)  we  see  that  the  average 
■jjyQ^)  only  through  its  magnitude 


probability  of  error  depends  on  the  moment 
Hence,  in  order  to  obtain  bounds  for  the 


V-V-~«0V-' 


error  probability,  it  is  necessary  to  characterize  the  range  of  values  of 
lf"(g)l  for  g  satisfying  (4.31).  Recall  from  Chapter  3  that  the  minimum  and 
maximum  values  of  a  continuous  linear  functional  are  found  by  searching  over 
the  set  of  extremal  g.  Since  lf"(g) !  represents  a  non-linear  functional  of  g, 
this  technique  cannot  be  directly  applied  to  find  the  minimum  and  maximum  of 
lf”l.  However,  it  can  be  shown,  through  an  application  of  the  triangle 
inequality,  that  the  maximum  of  lfM(g)l  for  g  satisfying  (4.31)  is  obtained 
for  some  extreme  g.  The  characterization  of  the  possible  values  of  If "I  is 
completed  by  assuming  that  the  minimum  of  |f"(g)|  is  zero.  Also,  notice  that 
if  *  ^i+i  “  ”1»  depends  only  on  f"  +  f"  “  2Re{f")  so  that  the 

problem  of  a  non-linear  functional  does  not  arise  in  this  case.  The  resulting 
bounds  for  the  limiting  error  probability  are  found  by  searching  over  the 
volume 

5  =  <  ^in’^a*1  X  tO.|f"lBax]  X  [2Re{f"}Bin.2Re{f»}Bax]  }. 

Numerical  evaluations  of  the  resulting  bounds  for  the  irreducible  error 
probability  reveal  that  the  maximum  variation  between  the  upper  and  lower 
bounds  for  fixed  Pr(h,0j)  occurs  for  the  case  h  =  1  and  =  0.  The  results 
for  this  case  are  shown  in  Fig.  4.6  as  a  function  of  Pr(l,0).  Notice  that  the 
maximum  variation  for  this  exaaq>le  is  still  less  than  a  factor  of  3.  We  point 
out  that  the  bounds  in  Fig.  4.6  not  only  represent  the  case  of  the  maximum 
variation,  but  also  represent  upper  and  lower  bounds  for  all  examples  of 
rectangular  pulse  waveforms  considered.  That  is,  the  bounds  for  rectangular 
pulse  FSK,  with  modulation  indices  h  ■  1,  2,  and  10  for  the  cases  0^,  =  o  and 
“  n,  all  lie  between  the  bounds  shown  in  Fig.  4.6. 
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4. 2. 2. 2  Phase-coded  pulse 

The  task  of  obtaining  bounds  for  the  irreducible  error  probability  for 
phase-coded  pulse  FSE  is  sore  difficult  than  for  DPSE  since  Pj.  for  FSE  depends 
on  essentially  twice  as  aany  parameters  (a  number  of  which  are  complex  valued) 
as  the  limiting  error  probability  for  DPSE.  Under  the  assumption  of 
adjacent-pul se-liaited  ISZ.  the  function  myy(-i,+l,-l)  Can  be  written  as 

■yyC-l.+l.-l)  -  |»J(h,eT)Hp  .  (4.32) 

As  in  the  case  for  phase-coded  pulse  DPSE,  there  is  no  obvious  way  to 
characterixe  the  moments  m^J^) ,  myyffcj) ,  and  m^y^)  in  terms  of  the 
normalized  rms  multipath  spread  iip(h,*)  and  a  few  linear  functionals  of  g. 
Also,  we  have  seen  that  m^+l  ,+l  ,+l) ,  which  is  closely  related  to  the 
parameter  defined  in  (3.22)  is  dependent  on  the  shape  of  the  delay  power- 
density  spectrum.  Moreover,  unless  the  rms  delay  is  zero,  both  m^yt+l ,+l ,+l) 
and  myy(+i,+i,+l)  can  be  significant  in  the  determination  of  the  irreducible 
error  probability. 

However,  representative  bounds  on  the  irreducible  error  probability  can 
be  obtained  by  assuming  a  specular  multipath  model  for  the  fading  channel,  as 
in  Chapter  3,  with  a  sampling  epoch  chosen  so  that  a  path  with  average  power 
of  at  least  1/2N  of  the  total  power  of  the  delay  spectrum  corresponds  to  a 
zero  relative  delay.  Ve  proceed  by  using  the  model  for  the  discrete  density 
spectrum  g(kTc)  end  rewriting  the  equivalent  signal-to-noise  ratio  Y(ii)  in 
(4.6)  as 

. (4.33) 
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The  bounds  for  the  irreducible  error  probability  are  found  by  characterizing 
the  minimum  and  maximum  values  of  the  parameters  (mp^(ji^)-«yy(^i) )  and 
(mp^(£i)ayy(j^)-|atjjy(&£)  !  )  for  each  combination  of  the  interfering  data  bits 
(1*4-1  *nd  bi+1)  under  the  constraint 


®yy("“l  #~1)  m 


Using  the  assumed  symmetry  properties  of  the  delay  spectrum  g(kT  ),  the  search 

c 

for  the  minimum  and  maximum  irreducible  error  probability  can  be  reduced  to 


the  six-dimensional  search  over  the  limiting  values  of  the  two  functions 

described  above  for  the  cases  kH«bi+1*l,  b i-l”b i+l“_1 '  and  bi-l*“bi+l*  Th® 
resulting  bounds  for  the  irreducible  error  probability  of  phase-coded  pulse 
waveforms  are  given  in  Tables  4.3  and  4.4,  respectively,  for  m-sequences  of 
length  N  -  7  and  N  ■  31  discussed  above.  Notice  that  the  variation  between 
upper  and  lower  bounds  is  not  large,  and  that  the  bounds  for  the  phase-coded 
pulse  are  very  similar  to  those  obtained  for  the  rectangular  pulse  above. 

4.2.3  £ff eg.t.1  si  AddlUY-9  Noise:  Characterization 

Consider  the  response  of  the  space  and  mark  filters  (producing  the 
sampled  outputs  U  and  Y)  to  a  sinale  transmitted  data  pulse  signal, 

s(t)  -  yy  v(t)e_^bOnllt/T  ;  0  ±  t  ±  T  , 

where  E  is  the  energy  per  data  bit.  In  the  absence  of  additive  noise  and 
assuming  a  space  signal  is  sent,  the  second  moments  of  the  outputs  of  the 
space  and  mark  filters  are  given  by 


16«2E 


£  g(?)  |*v(h,l, |{|) I2  d$ 


.  16o_E  f  0(Q)  ||v(0)|^»|V(a)|2]  d0  , 
T  a* 


(4.34a) 


Table  4.3.  Irreducible  error  probability  bounds  for  N  =  7  phase-coded  pulse 
FSK  with  respect  to  the  nornalized  rms  multipath  spread  |i  (h,0), 
for  h  *  7,  and  Oj  *  o 


H_(h,0) 


Lower  Bound 


Upper  Bound 


0.0005 

1.2499 .10”' 

1.2504.10  ' 

0.001 

4.9999.10"7 

5.0060.10"7 

0.005 

1.2499.10-5 

1 .2887 .10-5 

0.01 

4.9998.10-5 

5.6829.10"5 

0.05 

1.2488.10"3 

4.9551.10-3 

Table  4.4.  Irreducible 
with  respect 
h  -  31.  and 

error  probability  bounds  for  N  *  31  phase-coded  pulse 
to  the  normalized  rms  multipath  spread  u  (h,0),  for 
0T  =  0 

t*(h,0) 

Lower  Bound 

Upper  Bound 

0.0005 

0.001 

0.005 


1.2499.10- 7 
4.9999.10"7 

1.2499.10- 5 


1 .4131 •10-7 


9.2896.10' 


4.4880.10' 


97 


and 


g<?)  lKT(k.-l.UI)|2 

2 

160  E  ("  G(Q)  [V(Q-2nh/T)V*(0)®V*(0-2nh/T)V(Q)]  dfl  ,  (4.34b) 

T 


where  V(0)  is  the  inverse  Fourier  transform  of  the  basic  poise  shape  v(t),  and 
9  denotes  convolution.  Notice  that  (4.34a)  represents  a  measure  of  the  common 
spectra  between  the  data  poise  waveform  and  the  frequency  correlation  function 
G(Q).  Similarly,  (4.34b)  represents  a  measure  of  the  common  spectra  of  v(t), 
a  frequency  shifted  version  of  v(t),  and  the  correlation  function  G(Q).  The 
corresponding  output  of  each  filter  due  to  the  (real)  additive  white  Gaussian 
noise  with  (two-sided)  spectral  intensity  N^/2  is  8NqT. 


For  a  nonaeleotive  Rayleigh  fading  channel,  the  signal-to-noise  ratio  at 
the  output  of  the  space  filter  becomes 


snr  • 


0 


(4.35) 


since,  from  (4.10a),  ^(h,l,0)  »  T2.  Thus,  it  is  easy  to  see  that  the  one- 
shot  signal-to-noise  ratio  at  the  output  of  the  space  signal  filter  is 
maximized  when  the  rms  delay  spread  is  zero.  In  the  case  of  frequency- 
seleotive  fading,  the  signal-to-noise  ratio  at  the  output  of  the  desired 
signal  filter  is 


snr  * 


2  c2E  « 

N</  ^ 


g(t)  R2(h,l, I(|)  d( 


(4.36) 


Recall  from  (4.10a)  and  (3.14a),  that  Rj(h,l,$)  -  R^U) .  so  that  tjv  defined  in 
(3.12a)  for  the  consideration  of  DPSK  communications  is  also  the  hey  parameter 
in  determining  the  signal-to-noise  ratio  for  the  space  signal  filter  for 
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single— pulse  matched  filter  detection.  If  we  neglect  the  effects  of 
intersymbol  interference  and  assume  that  the  modulation  index  is  sufficiently 
large  that  the  signal  components  of  both  and  aTJV  are  approximately  zero, 
the  average  error  probability  (4.5)  becomes 


d  f  „  2g  E  2^v  -T1 
6  L  N0  T2  J 


(4.37) 


which  is  identical  to  (3.38)  except  for  a  factor  of  two  appearing  in  the 
signal-to-noise  ratio. 

The  results  in  Section  4.1  (see  Fig.  4.2)  demonstrate  that  for 
rectangular  pulse  FSK  with  integer— valued  modulation  index,  the  low  signal- 
to— noise  ratio  average  error  probability  in  frequency-selective  Rayleigh 


fading  is  well  approximated  by  Pn#  -  (2  +  2«2E/N0)_1  which  is  the  probability 


of  error  for  orthogonal  FSK  in  non-selective  fading.  While  it  is  not 
surprising  that  PQS  represents  a  lower  bound  for  the  average  error 
probability,  it  might  not  have  been  anticipated  that  this  performance 
parameter  for  orthogonal  signals  would  provide  an  adequate  low  signal-to-noise 
ratio  approximation  of  P^  for  selective  fading.  In  particular,  notice  from 
(4.35)  that  even  in  the  case  when  a  single  data  pulse  is  transmitted,  the 
relative  delays  introduced  by  the  frequency-selective  channel  produce  non-zero 
outputs  in  both  filters.  Thus,  in  general,  the  transmission  of  a  space  signal 
also  produces  an  output  due  to  this  signal  at  the  mark  signal  filter. 
However,  it  can  be  shown  that  PQj  is  an  accurate  low  signal-to-noise  ratio 
approximation  of  the  average  probability  of  error  in  selective  fading  for 
rectangular  pulse  FSK  for  any  integer-valued  modulation  index  h. 

The  discussion  above  indicates  that  the  low  signal-to-noise  ratio  error 
probability  is  largely  determined  by  the  output  due  to  the  desired  signal  at 


J 
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the  filter  Batched  to  that  signal,  which  in  torn,  depends  only  on  the 
paraaeter  tj  Thus,  we  expect  that  for  a  given  pulse  shape,  the  asymptotic 
error  probability  for  small  signal-to-noise  ratios  is  well  approximated  by 
(4.37)  which  is  independent  of  the  modulation  index  h.  Numerical  evaluations 
of  the  average  error  probability  in  (4.4)  and  P  defined  in  (4.37)  show  that 
this  is  the  case  for  both  rectangular  pulse  and  phase-coded  pulse  FSK,  for  all 
coabinations  of  delay  spectra,  aodulation  indices  and  phase  transitions 
considered.  In  fact,  the  average  probability  of  error  for  phase-coded  . %e 
FSK  is  both  well  approxiaated  by  (4.37)  for  saall  values  of  signal-ti  ..>ise 
ratio  and  lower  bounded  by  (4.37)  for  all  values  of  snr.  Hence,  just  ,  'n 
OPSK,  the  paraaeter  deteraines  the  position  of  the  "inverse-linear”  lower 
bound  on  FSK  error  probability  as  a  function  of  signal-to-noise  ratio  S  * 
2«*2E/N0.  Of  course,  the  factor  of  two  (or  3  dB)  difference  between  (4.37)  for 
FSK  and  (3.38)  for  DPSK  reflects  the  disadvantage  of  noncoherent  versus 
differentially  coherent  detection  in  Rayleigh  fading  channels. 

In  Table  3.7  of  Section  3.2,  the  values  of  if  /T2  are  listed  for  the 

v 

rectangular  pulse,  and  for  the  phase-coded  pulse  using  the  a-sequences  of 
length  N  »  7,  13,  31,  and  63  for  a  Gaussian  delay  power-density  spectrum.  It 
can  be  shown  that  for  phase-coded  pulse  FSK, 


(4.38) 


independent  of  the  aodulation  index  h,  where  N  is  the  length  of  the  signature 
sequence  and  the  constant  k  depends  on  the  channel.  Moreover,  the  discussion 
in  Seotioa  3.2  concerning  the  performance  degradation  of  the  phase-coded  pulse 
relative  to  the  rectangular  pulse  for  DPSK  applies  to  the  consideration  of 
phase— coded  pulse  versus  rectangular  pulse  FSK.  In  particular.  Table  3.7 


W'nV.V. 
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shows  that  for  a  Gaussian  delay  density  with  M/T  =  0.1,  the  additional 
signal-to-noise  ratio  requirement  for  phase-coded  pulse  FSK  relative  to 
rectangular  pulse  FSK  for  fixed  error  probability  is  about  S  dB  and  9  dB  for 
sequence  lengths  N  =  7  and  N  =  31,  respectively. 

In  Figure  4.7,  the  average  error  probability  for  several  examples  of  FSK 
signalling  formats  is  shown  as  a  function  of  the  signal-to-noise  ratio,  S. 
For  comparison  with  the  results  for  DPSK  in  Figure  3.7,  the  results  are  shown 
for  rectangular  pulse  FSK  with  modulation  indices  h  =  1,  and  10  (6^.  =  o)  and 
phase— coded  pulse  FSK  (N  3  h  —  31,  8^  —  0)  for  the  Gaussian  delay  power- 
density  spectrum  with  rms  delay  M  =  T/10.  Also  shown  are  the  low  signal-to- 
noise  ratio  asymptotes  given  by  (4.37),  and  the  irreducible  error 
probabilities  which  appear  as  horizontal  asymptotes.  (The  low  signal-to-noise 
ratio  error  probability  given  by  (4.37)  is  identical  for  both  examples  of 
rectangular  pulse  FSK).  Notice  that  in  contrast  to  the  results  for  DPSK  in 
Section  3.2,  the  limiting  error  probability  for  the  phase-coded  pulse  is 
considerably  larger  than  that  of  either  example  of  .rectangular  pulse  FSK. 
Also  notice  that  the  rectangular  pulse  with  h  =  1  exhibits  the  best 
performance  of  the  three  FSK  signals  considered  for  the  entire  range  of 
signal-to-noise  ratios.  These  results  are  in  agreement  with  the  discussion 
leading  to  (4.37)  which  shows  that  the  phase-coded  pulse  experiences  severe 
performance  degradation  in  the  presence  of  additive  noise;  and  the  data 
presented  in  Tables  4.1  and  4.2  which  show  that  the  normalized  rms  multipath 
spread  py(h,8j)  for  rectangular  pulse  FSK  (with  modulation  index  h  =  1,  and 
phase  transition  8^  *  0)  is  significantly  smaller  than  the  corresponding 
parameter  for  phase-coded  pulse  FSK. 


4.3  Performance  Approximations  for  Other  FSK  Systems 

In  Section  4.2  it  was  shown  that  meaningful  bounds  for  the  irreducible 
error  probability  of  FSK  in  a  WSSUS  frequency-selective  Rayleigh  fading 
channel  can  be  obtained  under  assumptions  that  correspond  to  the 
characteristics  of  many  practical  systems  [16.17].  These  bounds  demonstrate 
that  the  normalized  rms  multipath  spread  py(h,eT)  defined  in  (4.20)  is  the  key 
parameter  in  determining  the  limiting  error  probability  for  FSK  systems. 
Moreover,  the  closeness  of  these  bounds  (for  fixed  normalized  rms  multipath 
spread)  indicates  that  it  is  sufficient  to  specify  only  this  parameter  in 
order  to  obtain  an  accurate  estimate  of  the  irredncible  error  probability  for 
a  particular  fading  channel  and  FSK  signalling  format. 

Since  in  many  practical  systems,  the  limiting  error  probability  is  of 
primary  concern  in  the  evaluation  of  system  performance,  the  bounds  in  Section 
4.2  represent  a  characterization  of  the  limits  of  the  effects  of  frequency- 
selectivity  on  FSK  communications  for  a  wide  variety  of  signalling  formats  and 
frequency-selective  fading  channel  models.  As  in  the  case  of  DPSK  in  Section 
3.3,  these  results  also  imply  that  if  it  is  possible  to  measure  the  normalized 
rms  multipath  spread  parameter  for  a  proposed  FSK  system,  then  we  can  obtain 
an  estimate  of  the  average  probability  of  error  for  this  system  from  the  error 
probabilities  in  listed  in  Section  4.2. 

In  particular,  it  was  shown  in  Section  4.2  that  the  irreducible  error 
probability  as  a  function  of  the  normalized  rms  multipath  spread  for 
rectangular  pulse  FSK  for  any  integer-valued  modulation  index  and  any  phase 
transition  is  bounded  by  the  results  presented  in  Fig.  4.5.  Also,  comparing 
the  results  in  Tables  4.3  and  4.4  with  Fig.  4.5,  we  see  that  the  irreducible 
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error  probability  for  phaae-ooded  poise  FSK  as  a  function  of  |xp(b>Of)  closely 
approximates  the  irreducible  error  probability  for  rectangular  pulse  FSK  as  a 
function  of  |ir(h,8y).  Hence,  the  results  presented  in  Fig.  4.5  represent  what 
is  essentially  a  single  relationship  between  the  irreducible  error  probability 

Pj  and  the  normalized  rms  multipath  spread  ^(h.Oj)  for  FSK  in  a  frequency- 
selective  Rayleigh  fading  environment.  Thus,  we  can  take  the  data  in  Fig.  4.5 
as  representing  the  irreducible  error  probability  as  a  function  of  pv(h,0j) 
for  a  generic  FSK  signalling  format  with  time-limited  data-pulse  waveform 
v(t).  The  limiting  error  probability  for  a  particular  FSK  system  may  be 
approximated  by  means  of  a  single  channel  measurement. 


Hv(h,e) 


,(-l,+l,-l) 


(4.39) 


Using  this  measurement,  the  resulting  approximation  fox  the  irreducible  error 
probability  is  found  as  the  corresponding  value  of  Pj  in  Fig.  4.5. 

Notice  that  this  performance  parameter  is  more  complicated  than  the 
corresponding  result  in  (3.40)  for  DPSK.  Unfortunately,  it  is  not  clear  that 
a  less  complicated  form  of  channel  measurement  would  provide  information  that 
is  adequate  for  determining  the  limiting  error  probability.  However,  for  the 
special  case  of  rectangular  pulse  FSK,  the  results  in  the  previous  section 
indicate  that  an  accurate  estimate  can  be  obtained  from  a  less  complicated 
channel  measurement.  In  particular,  for  rectangular  pulse  waveforms, 
Bgij(+1,+1,+1)  is  independent  of  the  modulation  index  h  and  the  phase 
transition  6^  (of.  (4.18a)).  If,  in  addition,  the  modulation  index  is  large, 
the  channel  suy  be  characterized  by  measuring  the  value  of 


fT  1/2 

Hv(-,0)  -  [  T  g(«)*«  J 


(4.40) 


which  it  equal  to  both  the  normalized  rms  multipath  spread  for  rectangular 
pulse  DPSK  and  the  multipath  spread  p  =  M/T.  For  rectangular  pulse  FSK  with 
smaller  modulation  index  where  the  phase  transitions  are  modeled  as  uniformly 
distributed  random  variables.  (4.19)  and  the  results  in  Table  4.1  indicate 
that  (4.40)  also  gives  a  reasonable  estimate  of  the  limiting  error  probab  ility 
for  this  case. 

A  good  estimate  of  the  average  error  probability  for  small  time-bandwidth 
product  pulse  shapes  is  realized  as  a  curve  lying  above  Pq<  in  (4.12)  and  Pj 
calculated  from  (4.40)  and  the  data  in  Fig.  4.6  which  asymptotically 
approaches  Pn#  and  Pj  for  extreme  values  of  the  signal-to-noise  ratio.  For 
large  time-bandwidth  pulse  shapes,  a  more  accurate  approximation  may  be 
obtained  by  additionally  measuring  the  quantity 

tiD  9 

2V1*  *  JL  •  <«>  -  |  g({)  iJ(UI)  dt  .  (4.41) 

and  replacing  PQs  by  the  low  signal-to-noise  ratio  asymptote  given  in  (4.37). 

If,  in  the  absence  of  a  complete  channel  characterization,  it  is  possible 
to  obtain  close  estimates  of  the  parameters  discussed  above,  these 
approximations  could  be  used  in  the  preliminary  evaluation  of  FSK  systems 
employing  complicated  signals.  Alternatively,  if  a  fairly  accurate  model  for 
the  delay  power-density  spectrum  can  be  obtained,  the  definition  of  the 
normalized  rms  multipath  spread  (4.20),  the  data  in  Fig.  4.5,  and  the 
expression  for  the  low  signal-to-noise  ratio  error  probability  in  (4.37)  may 
be  used  for  comparative  evaluations  of  FSK  signalling  formats  and  data-pulse 
shapes. 

Finally,  we  point  out  the  several  basic  similarities  between  the  above 
results  for  the  approximations  of  FSK  and  the  approximations  for  DPSK 


discussed  in  Section  3.2.  In  particular,  the  low  signal-to-noise  ratio 
asymptotic  error  probabilities  for  any  particular  data-pulse  shape  are 
identical  for  DPSE  (3.38)  and  FSE  (4.37)  except  for  a  factor  of  two  increase 
in  effective  noise  power  for  FSE.  Also,  the  normalized  ns  multipath  spread 
for  rectangular  pulse  DPSE  is  similar  to  the  normalized  rms  multipath  spread 
for  rectangular  pulse  FSE,  especially  for  large  modulation  indices.  The 
relationship  between  the  low  signal-to-noise  ratio  performance  of  the  two 
systems  might  have  been  predicted  by  noting  that  DPSE  can  be  viewed  as  a  form 
of  binary  orthogonal  signalling  with  signal  duration  2T.  In  contrast,  the 
relative  effects  of  intersymbol  interference  on  the  two  binary  orthogonal 
signalling  systems  do  not  appear  to  be  related  to  the  basic  similarities 
between  system  models.  In  particular,  it  is  not  immediately  apparent  that 
rectangular  pulse  FSE  exhibits  the  same  sensitivity  to  intersymbol 
interference  as  rectangular  pulse  DPSE  only  for  very  large  modulation  indices. 
Nor  was  it  expected,  prior  to  the  development  in  Section  4.2,  that  it  is 
possible  to  choose  a  time-limited  pulse  waveform  for  DPSE  (sine-pulse)  that 
produces  lower  error  probabilities  than  any  standard  (i.e.,  rectangular  pulse) 
form  of  FSE  communications  over  the  same  fading  channel. 

4.4  Applications  j&  Ftff IMMrSt jOL  lifiilfi  Ch&afl.e  1.1 

In  this  section  it  is  shown,  in  a  development  that  closely  parallels 
Section  3.4,  that  when  the  received  signal  contains  a  nonfade d  specular 
component  as  well  as  a  Rayleigh  faded  component,  the  average  error  probability 
depends  on  the  same  hey  parameters  identified  in  the  preceding  sections  of 
this  chapter.  With  the  parameter  a  >  0  in  (2.3)  representing  the  strength  of 
the  specular  component,  the  decision  statistic  Z.  in  (4.2)  is  a  quadratic  form 
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of  (nonzero  mean)  Gaussian  random  variables  U  and  V.  As  in  tbe  analogous 
situation  for  DPSK  communications,  the  probability  of  error,  given  a 
particular  sequence  of  data  bits  Q^)  i*  transmitted,  is  given  by  a  value  of 
the  probability  distribution  of  the  ratio  of  Chi-square  variables  with  non¬ 
centrality  parameters  dependent  on  the  characteristics  of  the  FSK  system.  In 
the  remainder  of  this  section,  we  show  that  the  average  probability  of  error 
for  FSK  in  Rician  fading  is  completely  determined  by  the  parameter  a  in  (2.3) 
and  the  moments  “gptb^) ,  ayyCj)^) ,  and  defined  in  (4.7)  for  the 
evaluation  of  the  error  probability  for  Rayleigh  fading  channels. 


In  particular  if  we  assume,  without  loss  of  generality,  that  b. 


1  in 


(4.1),  so  that  a  space  signal  is  transmitted  during  the  i-th  signalling 


interval,  the  conditional  error  probability  [19]  can  be  written  as 


Q(  v  A(bi)  .  Vstfep) 


r  \/^mPP^i^'Hl>W^fei^2  ~  +  ^aUP^i)~lllW^— P  I 

2  \/(,BOT<fei)+®w(£i>>2  ”  4^mUV*-i^2 


A^+JKbj)  _ _ 

‘«P[ - 2 - J  *  ty  '/A(bi)B(ii))  , 


(4.42) 


where  Q(*,»)  is  Karcum's  Q  function,  and  I_(x)  is  the  modified  Bessel 


function.  The  parameters  AQ^)  and  Bd^)  in  (4.42)  are  given  by 


A^)-  8o2ET.[  app(fej)  *  aW(fej)  ~  \/(Buu^i>-t-">w<fei>>2  ~4*mpy(fei>  i2"  j  (4.43,) 

\A“jjjj(ki7+myy(fcj7)  ^  l®XJv(fei)  I2 


B(ki)«  8o2ET.[  aQU<fej)  *  ■w^j)  ♦  yCBu(ti)^n(ti»2  ~4*Bl}V(fei)  *2  J.  (4.43b) 

^*UB^i^+aw^i^^  ~  ^“ov^iiyi2 
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The  noaents,  in  (4.43)  ere  the  centre!  noaents  of  the  rendoa  veriebles 
V  end  V  defined  in  Chepter  3.  Using  (2.3)  end  (4.3)  for  the  signel  defined  by 

e 

(4.1)  end  essnaing  e  epece  signel  is  sent,  we  see  thet  E{U)  ■  a  \/  8ET  end  E(V} 
•  0  so  thet  the  noaents  in  (4.43)  ere  precisely  the  seae  es  those  defined  in 
(4.7)  for  the  enelysis  of  Reyleigh  feding  channels. 

It  is  eesy  to  see  thet  (4.42)  reduces  to  the  error  probebility 
expressions  for  Keyleigh  feding  (see  (4.6)  end  (4.8))  when  a  *  0.  In  this 
cese,  A(iA)  -  B(ii)  *  0,  end  (4.42)  becoaes 


2  y*“0U^i*+aW*ki*)2  “ 


r(l»i)  I 


],  (4.44) 


which  is  identical  to  (4.6)  for  the  ceee  b£  «  i.  Alternatively,  if  there  is 
no  feded  conponent  in  the  received  signal  so  thet  the  channel  is  siaply  the 
additive  white  Gaussian  noise  channel,  then  2c2  *  0  in  (2.4)  end  the  aoaents 
in  (4.43)  becoae 

“uu  “  “w  "  8N0T 

end 
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CHAPTER  5 

t 

ADAPTIVE  EQUALIZATION  TECHNIQUES  FOR  DIGITAL  COMMUNICATIONS 

OVER  FREQUENCY-SELECTIVE  FADING  CHANNELS 


In  Chapters  3  and  4,  the  performance  of  differentially  coherent  (DPSK) 
and  noncoherent  (FSK)  communications  via  VSSUS  frequency-selective  fading 
channels  was  considered.  It  was  shown  that  the  average  error  probability  for 
these  systems  is  degraded  by  two  effects  that  are  each  closely  related  to  the 
multipath  spread  of  the  fading  channel.  The  first  effect  is  that  of  loss  of 
available  signal  energy  due  to  time-dispersion  of  the  desired  signal  beyond 
the  sampling  window  of  the  matched  filter  receiver.  While  this  effect  causes 
minor  performance  degradation  for  small  time-bandwidth  product  signals  and  for 
adjacent-pul ae-limited  ISI,  it  can  become  a  major  consideration  for  large 
time-bandwidth  product  signals,  and  for  channels  which  are  more  severely 
frequency-selective  (dispersive  in  time). 


The  second  adverse  effect  is  due  to  the  intersymbol  interference 
introduced  by  the  superposition  of  delayed  versions  of  the  transmitted  signal 
at  the  channel  output.  The  results  in  previous  sections  demonstrated  that  it 
is  this  effect  which  most  severely  limits  the  performance  of  digital 
communications  over  frequency-selective  channels.  In  fact,  we  see  from 
Figures  3.7  and  4.7  that  for  rms  multipath  spread  p  as  small  as  0.05, 
intersymbol  interference  can  result  in  irreducible  error  probabilities  that 
are  unacceptably  large  for  many  applications.  Notice  that  the  parameter  p  ” 
M/T  is  a  linear  funotion  of  both  the  rms  delay  M  defined  in  (2.3)  and  the 
data-bit  rate  1/T.  Within  the  past  decade,  there  has  been  increased  interest 
in  applications  for  high-frequency  (HF)  and  troposcatter  communications 


[7,14.23,40]  for  which  the  rms  multipath  spread  p  -  0.5.  For  all  examples  of 
delay  spectra  and  modulation  formats  considered  in  previous  chapters,  the 
irreducible  error  probabilities  resulting  from  this  degree  of  frequency- 
selective  are  substantially  larger  than  10~*.  Thus,  in  order  to  achieve 
reliable  communications  over  dispersive  fading  channels,  it  is  necessary  to 
employ  more  complicated  receiver  structures  which  are  less  sensitive  to  the 
effects  of  intersymbol  interference. 

For  channels,  such  as  telephone  lines  and  line-of-sight  (LOS)  microwave 
links,  both  linear  and  nonlinear  equalizers  (typically  in  the  configuration  of 
a  tapped-de lay- line  (TDL))  have  been  effectively  used  to  reduce  the  effects  of 
intersymbol  interference  [20].  More  recently,  there  has  been  interest  in 
applications  of  similar  equalization  techniques  to  improve  the  performance  of 
digital  communications  over  classes  of  randomly  time-variant  linear  channels 
[8,23].  The  majority  of  channels  for  which  equalization  techniques  have  been 
considered  are  accurately  modeled  as  VSSUS  frequency-selective  Rayleigh  fading 
channels.  While  it  is  clear  that  the  random  character  of  the  fading  channel 
presents  a  number  of  additional  difficulties  (e.g.,  zero-mean  unequalized 
signal  amplitude)  which  must  be  overcome  through  the  use  of  diversity  or  some 
form  of  redundant  transmission,  both  simulation  studies  [8,24,41-43]  and 
experimental  evidence  [14]  indicate  that  for  a  given  data  rate,  reliable 
communications  can  be  achieved  for  fading  channels  with  rms  delays  which  would 
produce  unacceptable  limiting  error  probabilities  if  the  system  were  not 
equalized.  Alternatively,  for  a  fading  channel  with  a  given  value  of  rms 
mltipath  spread,  equalization  allows  an  increase  of  an  order  of  magnitude  in 
data  rates  for  a  specified  bit-error  probability  [7]. 


i 
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Studies  of  equalization  techniques  for  non-fading  dispersive  channels 
(e.g.,  high-speed  digital  telephone  links)  have  been  concentrated  on  linear 
equalizers,  in  general,  and  on  linear  TDL  equalizers,  in  particular  [21]. 
However,  one  of  the  characteristics  of  VSSUS  frequency-selective  channels  is 
that  the  tine-variant  transfer  function  frequently  exhibits  nulls  at  various 
portions  of  the  spectrun  which  result  in  "deep"  fades  of  the  spectrum  of  the 
received  signal  at  these  frequencies  [9,13].  The  primary  goal  of  the 
equalizer  is  the  elimination  (or  reduction)  of  the  effects  of  intersymbol 
interference.  More  generally,  the  goal  of  equalization  is  to  realize  an 
equalizer  filter  such  that  the  combination  of  the  channel  and  equalizer  yields 
an  equivalent  "channel"  with  flat  frequency  response  [44].  Since  practical 
linear  TDL  equalizers  can  only  realize  equalizer  filters  with  Z-transform 
transfer  functions  consisting  entirely  of  zeroes  [21],  linear  equalizers  are 
not  well  suited  for  applications  to  VSSUS  frequency-selective  fading  channels 
[8]. 

There  are  several  nonlinear  equalizer  structures  that  can  be  used  to 
compensate  for  channels  with  nulls  in  the  frequency  response  [21].  The 
nonlinear  equalizer  structure  that  has  received  the  most  attention,  because  of 
its  relatively  uncomplicated  implementation  and  demonstrated  effectiveness,  is 
the  class  of  TDL  decision-feedback  equalizers  (DFE)  [8,23],  The  form  of  DFE 
receiver  considered  for  applications  in  fading  channels  consists  of  a  linear 
TDL  equalizer  operating  on  the  received  signal  at  some  intermediate  frequency 
(or  passband  [45]),  a  filter  matched  to  the  transmitted  data  pulse  shape,  and 
an  additional  TDL  operating  on  the  results  of  the  (nonlinear)  decisions  of 
previous  data  symbols.  The  performance  of  the  TDL  equalizer  strongly  depends 
on  the  effectiveness  of  the  linear  TDL  portion  of  the  equalizer  [24,44].  In 


what  follows,  the  characteristics  of  DFE  equalizers  are  considered  by  first 
discussing  the  characteristics  of  linear  TDL  equalizers. 

5.1  Linear  Equalizer  Characteristics 

A  model  for  a  linear  tapped-de lay- line  equalizer  with  +  k2  +1  taps  is 
shown  in  Fig.  5.1.  Using  narrowband  signal  representation,  the  input  to  the 
equalizer  is  given  by 


r(t)  -  f  h(t,{)  s(t-$)  d?  +  n(t)  .  (5.1) 

LflD 

The  time-varying  channel  response  function  h(t,£)  is  assumed  to  be  that  of  a 
WSSUS  Rayleigh  fading  channel  as  discussed  in  Chapter  2  and  n(t)  is  the 
equivalent  low-pass  Gaussian  noise  with  (one-sided)  spectral  intensity  N^.- 

The  tap-weight  vectors  c  and  4  in  Fig.  5.1  represent  the  real  and 
imaginary  parts  of  the  equivalent  complex-valued  tap-weight  vector  x  =  c  -  jd. 
We  have,  for  simplicity,  dropped  the  dependence  on  the  time  reference  t.  The 
remainder  of  the  receiver  is  composed  of  in-phase  and  quadrature  correlation 
receivers  such  as  those  shown  in  Fig.  3.1,  which  operate  on  the  in-phase  and 
quadrature  outputs  I(t)  and  Q(t)  of  the  TDL  equalizer.  The  sampled  output  of 
the  composite  receiver  is  given  by 

T 

®  “  2  f  v(t)  ej0  r.(t)  dt  ,  (5.2) 

J0  ° 

where  r^(t)  “  I(t)  -  jQ(t)  is  the  output  of  the  TDL  equalizer  and  the  slowly 
time-varying  phase  0  represents  the  combined  phases  of  the  transmitter,  the 
channel,  and  the  receiver.  It  is  convenient  to  represent  the  output  U  in 


terms  of  the  vector  sum  of  the  equalizer  output  as 


'v  VS  V 


where  w'  represents  the  conjngete  transpose  of  the  column  vector  w,  T  is  the 

duration  of  the  data  pulse  waveform  v(t),  and  E  is  the  transmitted  energy  per 

bit.  The  tap-weight  vector  w'  in  (5.3)  is  given  by  w*  =  [w*  ,  . ,.,w* . 

~K1  0 

]  ,  where  (Kj+^+DT  is  the  width  of  the  time  window  of  the  received  signal 
r(t)  represented  by  the  equalizer  output  rQ(t).  The  vector  1  represents  the 
sequence  of  data  bits  which  have  influence  on  the  receiver  output  D.  Notice 
that  the  length  of  &  depends  on  the  number  of  equalizer  taps  as  well  as  the 
time  duration  of  the  channel  impulse  response  function  h(t,£). 

Simulation  studies  [24]  of  the  performance  of  linear  TDL  equalizers  are 
often  based  on  the  assumption  that  the  channel  response  to  an  impulse  at  t  =  0 
is  non-zero  only  for  relative  delays  5  in  the  range  [-T.T],  Notice  that  this 
is  .similar  to  the  adjacent-pulse-limited  ISI  assumption  used  in  Chapters  3  and 
4.  In  what  follows  we  assume,  for  purposes  of  discussion,  that  the  response 
function  h(t,{)  is  identically  zero  for  |t-$l  >  T.  Also,  notice  that  there  is 
some  ambiguity  as  to  the  time  reference  of  the  equalizer  receiver  output 
relative  to  the  transmitted  signal.  In  fact,  the  determination  of  time 
reference  (timing  phase)  is  a  major  concern  for  the  implementation  of  both 
linear  and  decision-feedback  equalizers  [23,46].  We  adopt  the  common  [24,44] 
assumption  that  the  time  reference  is  chosen  so  that  the  output  of  the  zero-th 
t*P«  vq,  corresponds  to  the  maximum  channel  output  due  to  the  desired  data 
bit,  bQ.  Under  the  above  assumptions,  the  data  sequence  &  affecting  the 
receiver  output  U  is  given  by 


IbK2+l'  bI2 


where  K^K^+Ij+l  is  the  total  number  of  equalizer  taps. 


The  matrix  R  in  (5.3)  contains  elements  which  represent  the  combined 
impulse  response  of  the  transmitted  waveform  v(t),  the  fading  channel  response 
function  h(t,£),  and  the  matched  filter  correlation  receiver.  Under  the 
adjacent-pul se-limited  ISI  assumption  R  is  given  by 


R  - 


R(-N,l),  R(-N.O),  R(-N.-l),  0  .  0. 

0  .  R(-N+l ,1) .  R(-N+l ,0) ,  R(-N+l,l) ,  0, 


0 

0 


.  (5.4) 


0 


R(N,1) ,  R(N.O),  R(N.-l) 


where 


R(n.k)  =  f  R(UI)  h(nT,  kT+£)  d£  .  (5.5) 

iT 

The  function  R(.)  in  (5.5)  is  the  aperiodic  autocorrelation  function  of  the 
time-limited  pulse  waveform  v(t)  defined  by  (3.14a).  It  is  common  to 
represent  the  combined  response  of  the  transmitter,  fading  channel,  and 
receiver  in  this  way.  In  fact,  the  rows  of  R  are  closely  related  to  the 
discrete-time  channel  response  functions  discussed  in  [9]  and  the  equivalent 
tapped— delay— line  model  for  doubly  selective  fading  channels  considered  in 
[10],  Note  that  the  first  argument  of  R(n,k)  denotes  the  value  of  this 
integral  for  the  channel  response  at  time  t  =  nT.  The  variation  of  the  value 
of  (5.5)  as  a  function  of  n  is  related  to  the  time-selectivity  of  the  fading 
channel.  For  the  present  discussion  we  assume  that  the  value  of  R(n,k)  is 
independent  of  n  for  n  €  [-Kj-l,  Kj+ll • 


One  of  the  key  characteristics  of  both  ELF  and  troposcatter  channels  is 
that  the  response  function  h(t.£)  is  time  varying.  When  the  changes  of  the 
characteristics  of  the  response  function  are  sufficiently  rapid  that  the  above 


assumption  does  not  apply,  the  channel  is  modeled  as  a  time-selective  fading 
channel  [9].  Typically,  the  degree  of  time-selectivity  is  measured  in  terms 
of  the  width  of  the  Doppler  power-density  spectrum  which  is  the  temporal 
Fourier  transform  of  the  function  p(t-x,*)  defined  in  (2.4).  Definitions  of 
this  type  which  are  commonly  used  in  the  published  literature  are  the  distance 
between  the  1/e  points  or  the  half-power  bandwidth  of  the  Doppler  spectrum 
(e.g.,  see  [3,9]).  Alternatively,  the  time-selectivity  of  the  VSSUS  fading 
channel  may  be  grossly  characterized  by  the  smallest  value  of  t  for  which 
p(T,0)  «  0.  This  value,  which  is  inversely  proportional  to  the  width  of  the 
Doppler  spectrum,  is  referred  to  as  the  channel  "decorrelation  time"  [7,24]. 

Since  both  HF  and  troposcatter  fading  channels  are  time-selective  to  some 
degree,  a  time-varying  equalizer  must  be  employed.  For  applications  to  WSSUS 
fading  channels,  the  equalizer  must  adapt  to  the  time-varying  characteristics 
of  the  channel.  One  of  the  earliest  applications  of  adaptive  equalizers  was 
to  high-speed  data  communications  over  switched  telephone-line  networks  [22]. 
In  this  situation,  the  receiver  is  designed  to  adapt  to  one  of  a  number  of 
links  before  the  data  is  transmitted.  Equalization  of  fading  channels 
requires  continuous  adaptation  in  step  with  the  changing  characteristics  of 
the  fading  channel.  Moreover,  it  is  crucial  that  the  equalizer  adapt  in  a 
time  period  which  is  much  shorter  than  the  channel  decorrelation  time.  For 
many  applications,  this  implies  that  the  decorrelation  time  must  be  much 
larger  than  the  data  pulse  duration  T.  Fortunately,  this  constraint,  which  is 
referred  to  as  a  learning  constraint  [40],  does  not  preclude  the  use  of 
equalization  techniques  for  most  fading  channel  applications.  Typical  values 
of  the  decorrelation  time  are  on  the  order  of  10^ *T  for  HF  systems  and  10^'T 
for  troposcatter  links  [7]. 


In  order  to  reduce  the  effects  of  intersymbol  interference,  the  tap- 
weights  w  are  adjusted  (adaptively)  to  minimize  some  cost  function  which, 
hopefully,  is  related  to  the  resulting  system  performance.  When  the  receiver 
bases  the  binary  decision  for  each  successive  data  bit  on  a  single  sampled 
output,  the  cost-function  that  is  most  often  employed  is  the  mean-square-error 
(USE)  [8,23.44]. 

The  consideration  of  the  performance  of  adaptive  equalization  techniques 
for  slowly  time-varying  fading  channels  must,  it  seems,  be  based  on  the 
following  basic  assumption.  At  each  time  t  and  for  each  sample  function 
h(t,()  of  the  channel  impulse  response,  the  TDL  is  adjusted  in  a  manner  such 
that  the  equalizer  is  able  to  "track"  the  slow  changes  in  the  channel.  In 
this  case  it  is  not  clear  that  knowledge  of  the  channel  statistics  (zero-mean 
Gaussian  for  a  Rayleigh  fading  channel)  can  be  utilized  to  obtain  a  meaningful 
interim  measure  of  system  performance.  On  the  other  hand,  if  we  assume  that 
the  equalizer  is  unable  to  track  the  channel,  then  the  TDL  would  only  serve  to 
add  additional  intersymbol  interference  as  well  as  correlated  Gaussian  noise. 
In  this  case  it  is  reasonable  to  assume  that  it  would  be  better  not  to  attempt 
channel  equalization. 

Consider  the  receiver  output  U  for  a  fixed  channel  impulse  response 
function  h(>,(),  which  is  a  function  only  of  the  relative  delay  £.  The  mean- 
square-error  (MSE) ,  as  defined  for  the  consideration  of  fading  channels 
[8,24,45]  and  for  the  consideration  of  deterministic  channels  [45,47],  is 
given  by 

MSE  -  E^{|u  -  kV  8ET  bQ|2}  ,  (5.6) 
where  E  is  the  transmitted  energy  per  bit  and  E^f.}  denotes  the  conditional 


expectation  given  R.  The  expectation  in  (5.6)  is  taken  with  respect  to  the 
data  sequence  for  i^O  and  the  additive  noise  samples  a.  The  scale  factor  k 
in  (5.6)  is  a  positive  real  constant  vhich  does  not  affect  the  form  of  the 
minimum  MSE  filter.  Using  the  vector  representation  in  (5.3)  for  U,  we  have 

MSE  -  Er(|  /Sp  l'  1  i  +  2*'  a  -  k /8ET  bQ|2}  .  (5.7) 


If  we  assume  that  the  data  sequence  is  a  sequence  of  independent,  identically 
distributed  random  variables  taking  on  the  values  -1,  +1  with  equal 
probability,  (5.7)  becomes 

MSE  =  f  I’  I  I  I'  S  +  8N0T-w'w  +  k2 •  8ET  -  16kE.Re [w'Rg]  ,  (5.8) 

where 

*  [R  (*»“’ ^)  ,  • . .  ,R  <».0>,  ...,  R  (‘tlj)]  . 

It  can  be  shown  [23,44,45]  that  the  solution  for  the  MSE  optimum  tap— weight 
vector  w  is  given  by 

w'  A  w  “  k.  Re(w*  Rq}  (5.9) 

or 


where  the  Hermitian  matrix  A  i*  written  as 


A 


T”1  S  I' 


(5.10) 


Using  (5.9)  and  (5.8),  the  minimum  MSE  (MMSE)  is  given  by 
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Alternatively,  the  MMSE  can  be  written  aa 


WISE  -  kz8ET  -  k8E.Re{w'  R.) 


(5.12) 


Notice  that,  for  a  particular  channel  response  matrix  R,  the  second 


sioment  of  the  filter  output  D  (for  b,  =  i)  ia  (5.3)  is  related  to  the  mean- 


square-error  by 


EjdDl2)  *  MSE  -  k28ET  +  16kE.Re{w'  R^ 


(5.13) 


Using  the  expression  for  the  MMSE  in  (5.12)  we  have 


ER(|U|2)  -  MMSE  -  k28ET  +  16kE.Re{w'  R^ 


(5.14) 


ES(|UI2)  -  k8E.Re{w'  !„}  >  0  . 


(5.15) 


We  have  seen  that  the  MSE  or,  in  particular,  the  MMSE  is  related  to  the 


second  moment  of  the  receiver  output  U  (cf.  (5.14)).  It  is  tempting  to 


consider  the  possibility  of  characterizing  an  "average"  MMSE  by  taking  an 


additional  expectation  with  respect  to  the  channel  statistics.  Using  the 


expression  for  MSE  in  (5.8)  and  averaging  with  respect  to  the  channel  ensemble 


(WSSUS  zero-mean  Gaussian)  we  have 


MSE  -  w'  C  w  +  8NqTw'w  +  k28ET  , 


where  the  matrix  C  denotes  the  expected  value  of  the  quadratic  form  RR' .  If 


there  is  no  additional  restriction  on  w,  then  this  implies  that  w  =  0  produces 


the  minimum  MSE.  Alternatively,  if  we  insist  that  s'w  =  1 ,  then  it  is  easy  to 


see,  using  an  example  where  the  fading  channel  is  nonselective,  that  the  above 


minimization  produces  the  wrong  result. 


w.’.  .‘voviv  :v 
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It  is  interesting  to  note  that  the  last  tern  in  (5.8),  which  is  also  the 
right-hand  side  of  (5.15),  is  jnst  a  constant  times  the  real  part  of  the 
output  D  due  to  the  desired  signal  (corresponding  to  bg»l).  Then  we  see  that 

(5.15)  implies  that  the  optimum  tap  vector  w  produces  a  "one-shot"  receiver 
response  whose  real  part  has  non-zero  mean-value.  Hence,  we  see  that  a  linear 
TDL  equalizer  designed  to  minimize  the  MSE,  not  only  reduces  the  average 
contribution  of  interfering  adjacent  data  signals,  but  also  provides  phase- 
coherence  between  the  transmitter  and  receiver.  Moreover,  using  (5.8)  and 

(5.15) ,  it  can  be  seen  that  for  small  values  of  E/N^,  the  linear  TDL  equalizer 
receiver  is  closely  related  to  the  maximal  ratio  diversity  combiner  for  a 
one-shot  receiver  [48]  where  the  "diversity"  links  are  produced  by  the  time- 
dispersive  properties  of  the  channel.  In  a  frequency-selective  fading 
environment,  the  received  signal  can  be  characterized  as  the  superposition  of 
delayed  versions  of  the  transmitted  signal.  Thus,  the  efficient  combining  of 
these  signal  components  at  the  receiver  results  in  a  potential  for  diversity 
gain.  The  diversity  that  is  due  to  the  time-dispersive  nature  of  the 
frequency-selective  fading  channel  is  commonly  referred  to  as  "implicit" 
diversity  [25]  as  opposed  to  explicit  forms  of  diversity  which  can  be  realized 
by  redundant  transmissions  over  channels  separated  by  space,  time,  frequency, 
or  antenna  polarization  [6] .  (Recall  that  maximal  ratio  combining  is  a 
linear,  coherent  diversity  combining  technique  which  achieves  the  maximum 
possible  output  signal-to-noise  ratio  for  a  single  data  pulse  at  each  sampling 
instant.)  If  the  matrix  &  in  (5.10)  is  dominated  by  the  term  due  to  the 
additive  noise,  and  if  the  contributions  from  interfering  data  signals  are 
ignored,  the  statistic  U  in  (5.3)  becomes 


u  *  w'  So  *0  +  2*'  fl 
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(5.16) 


With  U  defined  in  this  way,  the  signel^to-noise  ratio  is  given  by 


SNR  ”  — —  d»  D. 
H.I*  ®°*° 


which  is  precisely  the  definition  of  the  instantaneous  signal-to— noise  ratio 


for  aaxiaal  ratio  coabining  considered  in  [6]  and  [48]  where  the  order  of  the 


iaplicit  diversity  is  related  to  the  nuaber  of  significant  eleaenta  of  ]L.  in 


light  of  the  relationship  between  a  MMSE  linear  IDL  equalizer  and  coherent 


aaxiaal  ratio  coabining,  it  is  not  reasonable  to  assuae  that  it  is 


advantageous  to  consider  signalling  techniques  other  than  coherent  PSK.  In 


fact,  it  is  well  known  that  coherent  binary  PSK  produces  smaller  error 


probabilities  than  other  coaaon  fora  of  binary  signalling  (i.e.,  DPSK,  and 


coherent  and  noncoherent  FSK)  for  aaxiaal  ratio  combining  in  a  Rayleigh  fading 


environment  [6]. 


Ve  have  seen  that  for  a  given  channel  impulse  response,  the  MMSE-TDL 


equalizer  allows  the  use  of  coherent  PSK  detection.  Using  the  solution  for 


the  MSE-optimal  w  in  (5.9)  for  k  *  1,  the  random  variable  U  in  (5.3)  can  be 


written  as 


D  A"1  at  +  2*6  A-1  a  • 


(5.17) 


where  £,  and  A  are  all  determined  by  the  channel  response  function  and  the 


statistics  of  the  additive  noise.  Conceptually,  the  conditional  bit-error 


probability  for  a  particular  h(»,$)  and  a  particular  data  sequence  &  can  be 


obtained  by  finding  the  signal-to-noise  ratio  and  coaputing  P@ (j^) *Q( SNR1^) 


where  Q(x)*l-F(x) :  the  standard  Gaussian  cumulative  distribution  function. 


The  average  probability  of  error  is  then  determined  by  averaging  Pe(£) 


over 


✓  / 
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the  possible  date  sequence  combinations.  Unfortunately,  the  evaluation  of 
both  the  conditional  and  average  error  probabilities  is  quite  cumbersome  even 
for  non-fading  channels  and  the  results  concerning  the  performance  of  MUSE 
equalizers  are  almost  always  in  the  form  of  simulation  studies  [21,44].  An 
exact  evaluation  of  the  average  error  probability  for  fading  channels  appears 
hopeless,  since,  in  this  case,  it  is  necessary  to  perform  additional  averaging 
with  respect  to  the  channel  response  function  h. 

5.2  Decision-Feedback  Equalizer  Characteristics 

The  receiver  structure  of  a  DFE  is  similar  to  that  of  a  linear  TDL 
equalizer  receiver  with  the  addition  of  a  feedback  loop  which  forms  a  linear 
combination  of  the  values  of  previously  detected  data  symbols  [8,23,42.49,50]. 
(The  nonlinear  nature  of  the  equalizer  is  due  to  the  effect  of  the  hard 
decisions  used  to  form  the  output  of  the  feedback  section.)  The  device  used 
to  form  the  feedback  output  D  can  be  (and  usually  is  [8])  represented  as  a 
linear  TDL  with  tap  spacing  necessarily  equal  to  the  data  symbol  duration  T. 
If  we  assume  the  current  decision  concerns  the  zero-th  data  bit  b^,  the  output 
of  the  feedback  section  of  the  DFE  is  given  by 

D  -  T  d i  b_ ,  ,  (5.18) 

S-l  1  1 

where  N  is  the  number  of  feedback  taps  and  is  the  value  of  the  decision  on 

the  i-th  data  bit.  Using  vector  representation  for  the  parameter  D,  the 

expression  for  the  decision  statistic  U  for  a  DFE  can  be  written  as 

U  "  w'  £f  b  +  4'^  +  2*»  a  .  (5.19) 

where  4'  ”  . . .  ,d^]  and  ■  [b_^, . . .  ,b_^] .  Notice  that  there  are  now  Kj  + 
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1  taps  for  the  linear  TDL  portion  of  the  equalizer  so  that  1^=0 .  The  matrix 

£f  in  (5.19)  is  obtained  from  (5.4)  by  deleting  the  last  Kj  +1  columns  and  the 

l**t  ^  rovs  of  the  original  channel  response  matrix. 

The  USE  criterion  for  the  DFE  is  defined  similarly  as  for  the  linear  TDL 
equalizer  and  is  based  on  the  assumption  [8,24,44]  that  the  previously 
detected  symbols  in  the  feedback  filter  are  correct,  i.e.,  6^  =  for  i<0. 

Using  this  assuaiption  and  the  definition  in  (5.6),  the  NSE  for  the  DFE  is 
given  by 

MSE  “  Ej^{|  \]8&.  £f  bf  +  +  2w'  n  -  k  >/8ET  *b0|2}  *  <5.20) 

Minimization  of  (5.20)  with  respect  to  w  shows  that  the  linear  TDL  that 

minimizes  the  MSE  is  given  by  the  relation  [44] 

S'  A f  S  *  k*  Me  (s'  Af©] 
or 

S  "  k'A"1  Rf0  ,  (5.21) 

where  the  matrix  is 

N  T 

Af  =  t"1  Sf  ££  +  jT“  •!  •  ‘  (5.22) 

The  coefficients  of  the  DFE  are  then  chosen  so  that  the  resulting  output  D 
exactly  cancels  the  intersymbol  interference  due  to  previously  detected 
symbols  or  "precursors"  [23].  In  this  case,  the  coefficients  d^  can  be 
written  in  terms  of  the  coefficients  of  the  linear  TDL  section  by 

4  -  S'  Sj  .  (5.23) 

where  the  matrix  is  obtained  from  (5.4)  by  deleting  the  first  Kj  coluams 
and  the  first  I2+i  rows  of  g. 
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For  any  practical  implementation  of  DEE  receivers  (especially  for 
Rayleigh  fading  channels)  it  is  not  reasonable  to  assume  that  all  previous 
data  decisions  are  correct.  Moreover  in  the  case  when  a  decision  error  is 
auuie,  the  output  from  the  feedback  portion  of  the  equalizer  adversely  affects 
subsequent  decisions  causing  error  propagation.  For  nonfading  channels,  the 
increase  in  the  average  error  probability  due  to  this  effect  is  considered  in 
[20]  using  a  Markov  chain  model.  The  analysis  is  carried  out  under  the 
assuaption  that  the  tap  weights  for  the  linear  TDL  forward  filter  are  chosen 
according  to  a  zero-forcing  (ZF)  criterion  [22,44].  The  ZF  cost  function  is 
closely  related  to  the  definition  of  the  MSE  in  (5.7)  by  neglecting  the  term 
corresponding  to  additive  channel  noise.  The  zero— forcing  criterion  is  simply 
the  requirement  that  the  intersymbol  interference  contribution  at  the  output 
of  the  forward  filter  is  forced  to  zero,  at  the  cost  of  an  enhancement  of  the 
additive  noise.  It  should  be  pointed  out  that  the  zero-forcing  criterion  is 
not  generally  applicable  to  fading  channels  since  ZF  can  be  achieved  only  for 
channels  with  limited  intersymbol  interference  [44].  Assuming  ideal  DFE 
characteristics  and  a  particular  nonfading  channel  response  function,  it  is 
shown  in  [20]  that  error  propagation  is  not  a  major  source  of  degradation  for 
typical  values  of  signal-to-noise  ratio  and  channel  characteristics.  A 


similar  conclusion  is  obtained  in  [40]  where  a  simulation  study  of  the  error 
propagation  effect  is  carried  out  for  a  single-echo  and  an  exponentially 
decaying  multiple  echo  channel  model. 

From  the  discussion  of  linear  TDL  equalizer  characteristics  in  Section 
5.1,  we  see  that  the  forward  filter  portion  of  the  DFE  has  the  potential  of 
providing  a  degree  of  implicit  diversity  gain  as  well  as  some  reduction  of  the 
effects  of  intersymbol  interference,  at  the  expense  of  an  increase  in  the 
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effective  noise  power.  In  contrast,  the  feedback  portion  of  the  DFE  is 

designed  to  eliminate  the  distortion  dne  to  previously  transaitted  data 
• 

signals  at  the  expense  of  a  loss  of  diversity  gain  and  soae  degradation  due  to 
error  propagation.  Hence,  the  DFE  receiver  can  be  viewed  as  a  compromise 
between  these  two  configurations. 

Finally,  we  point  out  that  for  both  linear  TDL  and  DFE  receivers,  the 
ability  to  track  the  channel  characteristics  strongly  depends  of  the  available 
signal-to-noise  ratio  [44].  Clearly,  if  the  received  signal  is  dominated  by 
the  additive  noise,  there  is  little  hope  of  maintaining  equalization.  In  a 
Rayleigh  fading  environment,  the  instantaneous  SNR  can  fall  below  the  level 
necessary  to  maintain  equalization  so  often  that  the  equalizer  performance  is 
severly  degraded.  In  order  to  overcome  this  difficulty,  practical  DFE  modems 
are  designed  to  process  signals  received  from  L  explicit  independent  diversity 
links  [23].  This  is  accomplished  by  employing  a  separate  forward  TDL  filter 
for  each  of  the  L  explicit  links,  appropriately  combining  the  forward  filter 
outputs  by  means  of  automatic  gain  control  devices,  and  using  a  common  data 
pulse  filter  and  feedback  loop  TDL  to  process  the  combined  output  of  the 
forward  filters. 

5.3  DFE  Performance  Considerations 

The  published  analytical  work  on  the  evaluation  of  the  average  bit-error 
probability  for  DFE  receivers  in  frequency-selective  fading  is  almost  entirely 
due  to  Honsen  [24,40].  The  analysis  in  [24]  is  carried  out  for  a  system 
without  explicit  diversity  and  is  based  on  the  fundamental  assumption  that  the 
feedback  portion  of  the  equalizer  precisely  cancels  the  intersymbol 
interference  due  to  pulses  transmitted  prior  to  bQ  (the  desired  bit).  Under 


this  assumption,  it  is  easy  to  see  that  the  resulting  analysis  is  independent 
of  the  parameters  of  the  feedback  portion  of  the  equalizer.  Hence,  the 
development  in  [24]  can  be  viewed  as  an  analytical  technique  that  can  be 
applied  to  linear  TDL  equalizers. 

Notice  that  the  receiver  output  U  in  (5.3)  can  be  written  as 

0  -  \/?w'  lf0b0  +  \ m*'  !f£f  +  2w'n  ,  (5.24) 

where  and  are  obtained  by  deleting  the  column  of  Rf  and  the  entry  of  &£ 
corresponding  to  the  desired  bit  bQ.  The  first  term  of  (5.24)  then  denotes 
the  desired  signal  component,  the  second  term  is  the  output  due  to  interfering 
signals,  and  the  last  term  is  due  to  the  additive  channel  noise.  In  [24], 
Nonsen  approximates  the  interfering  data  bits  b^  for  i  >  0  as  independent 
zero-mean  complex  Gaussian  random  variables  with  variance  y2.  (Recall  that 
the  intersymbol  interference  due  to  previous  data  pulses  is  assumed  to  be 
eliminated  by  the  feedback  portion  of  the  equalizer.)  Vith  the  intersymbol 
interference  from  future  data  pulses  modeled  as  additive  Gaussian  noise,  the 
minimization  of  the  MSE  with  respect  to  the  coefficients  of  the  forward  filter 
yields 


2'  I  *  -  k.Re{w'  If0}  , 

where  the  matrix  g  is  given  by 

B  -  f  V  If  +  ^.1  .  (5.25) 

Substituting  the  solution  for  the  MMSE  forward  filter  tap  weights  into  (5.24), 
the  statistic  U  becomes 

\/^  ifoB  1®f0b0  +  ft  -fOB  1^f^  +  2-f()fi  • 


U 


(5.26) 
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Under  the  Gaussian  noise  characterization  for  the  interfering  data 
pulses,  we  see  that  the  last  two  terms  in  (5.26)  represent  the  contributions 
due  to  additive  noise.  Monsen  defines  an  “equivalent  signal-to-noise  ratio" 

as 


SNR  = 


™  ( R*  R-^  R  )  ^ 

x  t&fo  —  — fp< 

8N0T  1*0  fi"1  If0 


-1 n  ,  R_1  R 

NqT2  Zto  fi  fifO  ' 


(5.27) 


which  is  the  ratio  of  the  output  due  to  the  desired  data  pulse  to  the  average 


of  the  output  due  to  both  interfering  data  pulses  and  additive  noise.  Hence, 
we  see  that  the  intersymbol  interference  has  been  replaced  with  interference 


due  to  additional  additive  Gaussian  noise. 


The  second  approximation  employed  in  [24]  is  to  replace  the  matrix  B  by 

its  mean  value,  E(B).  However,  even  with  these  simplifying  assumptions,  it  is 

still  necessary  to  find  the  eigenvalues  of  the  correlation  matrix  ECR^r^}  of 

zero-mean  correlated  complex  Gaussian  random  variables  with  means  and 

covariances  that  depend  on  the  channel  response  function  h(.,£).  When  the 

eigenvalue  distribution  is  obtained  for  a  particular  channel  realization,  it 

is  then  necessary  to  find  the  coefficients  of  a  partial  fraction  expansion  of 

y 

a  polynomial  in  the  eigenvalues.  The  variance  y  of  the  interfering  data  bits 
is  chosen  in  [24]  to  be  0.5  to  allow  the  best  agreement  between  calculated  and 
simulated  results. 


Since,  in  reality,  the  average  bit-error  probability  is  dominated  by  the 
probability  of  a  wrong  decision  for  the  worst  case  intersymbol  interference, 
it  would  be  useful  to  supplement  the  analysis  in  [24]  by  characterizing  the 


receiver  output  D  for  the  worst-case  data  sequence.  Using  (5.24)  we  have 


D  "  SfO  b0  +  If  £  +  2E'fi  •  (5.28) 

so  that  w'(&£Q  +  R^.  £)  represents  the  sampled  output  of  the  receiver  for  a 
particular  data  sequence.  If  we  assume  h(.,()  =  0  for  1(1  >  T  and  that  the 
forward  TDL  filter  consists  of  only  2  taps,  (R^  +  r  £)  is  given  by 


R( • ,-1)1  r b2R( • ,1)  +  b1(R(.,0) 

R( • ,0)  b1R(.,l) 


(5.29) 


Unfortunately,  even  for  this  simple  example  it  is  not  clear  that  the  worst 
case  interference  and  the  resulting  effect  on  the  error  probability  can  be 
adequately  characterized. 


5.3.1 


Fading  Channel  Approximation  of  DFE  Performance 


Since  the  forward  TDL  equalizer  in  each  explicit  diversity  branch  is 
assumed  to  achieve  coherence  with  the  channel,  an  interesting  approximation  of 
DFE  system  performance  (which  proves  to  be  a  reasonable  estimate  of  the 
performance  of  practical  DFE  receivers)  is  the  error  probability  of  coherent 
PSK  for  maximal  ratio  combining  of  independent  nonselective  Rayleigh  fading 
channels.  If  L  denotes  the  order  of  the  explicit  diversity,  the  probability 
of  error  of  maximal  ratio  PSK  is  given  by  [6], 


pc(L) 


Q(  V  2s)  *f|  (s)ds  , 
0  L 


(5.30) 


where  s  is  the  instantaneous  combined  signal-to-noise  ratio  and  f^(s)  is  the 
probability  density  function  of  s  for  diversity  L.  Assuming  identical 
diversity  branch  statistics  (i.e.,  average  branch  signal-to-noise  ratios),  the 
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probability  density  function  f^(s)  is  [6] , 


*»<•) 


.L-l 


<L-U!  ‘  SL  ‘ 


exp(-s/S)  , 


(5.31) 


where  the  common  per  branch  average  signal-to-noise  ratio  S®2<t^E/Nq.  The 


results  of  the  computation  of  PC(L)  as  a  function  of  S  are  shown  in  Fig.  5.2 
for  diversity  L  -  1,  2,  and  4.  As  a  comparison  for  the  results  in  Chapters  3 
and  4  for  the  case  of  diversity  reception,  the  error  probability  for  square 
law  diversity  combining  of  DPSK  in  nonselect ive  Rayleigh  fading  is  also  shown 
as  a  dashed  line  for  the  same  orders  of  diversity.  When  the  decision 
statistic  for  DPSK  signalling  is  of  the  form  of  Z  in  (3.5),  "square  law" 
combining  should  be  viewed  as  quadratic  form  combining  (see  [6])  of  the 
received  signals  given  by  the  sum  of  L  independent  variables  for  n=l,L 
[17].  Since  DPSK  signalling  can  be  represented  as  a  particular  form  of 
orthogonal  signalling  on  the  interval  [0,2T],  this  quadratic-form  combining, 
which  is  the  optimal  combining  technique  for  DPSK  in  Rayleigh  fading  if  no 
channel  estimate  is  available,  is  equivalent  to  square  law  combining  of  the 
orthogonal  signals.  It  can  be  shown,  using  the  results  in  [6]  or  [48],  that 
the  average  error  probability  for  square  law  DPSK  in  Rayleigh  fading  is  given 


by 


Pd(L)  *  PeL  *  f  1  (1-Pe)n  C(L-l+n, n)  , 
n*0 


(5.32) 


.-1 


where  Pe=(2+2S)  i*  the  average  error  probability  for  DPSK  in  nonselective 
Rayleigh  fading  (cf.  (3.18)),  and  C(m,n)  is  the  binomial  coefficient.  The 
corresponding  result  for  orthogonal  FSK  can  be  obtained  from  Fig.  5.2  by 
finding  the  value  of  the  error  probability  corresponding  to  a  3  dfi  reduction 
of  signal-to-noise  ratio  S.  The  effect  of  the  limiting  error  probability  on 
the  square  law  diversity  performance  of  DPSK  and  FSK  can  be  approximately 


and  square  law  DPSK  for 
ratio,  S 


determined  by  finding  tbe  irreducible  error  probability  for  a  particular  value 
of  rms  multipath  spread  |i  from  Fig.  3.3  or  Fig.  4.3,  which  is  the  exact  result 
for  L*l.  The  approximation  of  the  limiting  error  probability  for  L=2,4  is 
then  obtained  by  finding  the  corresponding  value  of  S  for  the  square  law,  L=1 
curve  and  determining  Pd(L)  with  L=2  and  L=4  for  this  value  of  S. 

5.4  Simulation  and  Experimental  Results  on  DFE  Performance 

5.4.1  Simulation  Results 

One  of  the  earliest  published  records  of  computer  simulation  results  for 
the  performance  of  DFE  for  fading  channels  appears  in  [23].  For  this 
simulation  study,  the  Rayleigh  fading  channel  is  modeled  with  a  rectangular 
delay  power-density  spectrum  (cf.  (2.16))  with  the  parameter  TQ  ranging 
between  0  and  0.5/T  where  T  is  the  data  pulse  duration.  Notice  that  for  all 
values  of  TQ  considered  in  this  simulation,  the  frequency-selective  channel 
satisfies  the  adjacent-pulse-limited  ISI  assumption  for  the  received  signal 
prior  to  equalization.  From  the  definition  of  the  rms  multipath  spread  in 
Chapter  2,  we  see  that  these  values  of  T^/t  correspond  to  rms  multipath 
spreads  ranging  from  zero  (i.e,,  nonselective  fading),  to  a  maximum  value  of 
0.3  for  Tq  =  0.5/T.  The  simulation  is  carried  out  for  a  DFE  consisting  of  a 
single  three-tap  forward  filter  with  data  pulse  interval  tap  spacing,  and  a 
three-tap  feedback  filter  (i.e.,  K^+i  =  N  =  3).  The  results  of  this  study 
suggest  that  the  implicit  diversity  gain  resulting  from  the  correct  resolution 
of  the  multipath  signal  components  can  be  realized  in  a  simulated  fading 
environment.  This  implicit  diversity  gain  produced  error  probabilities  less 
than  the  probability  of  error  for  coherent  detection  in  nonselective  fading 
(i.e.,  Pc)  by  *n  order  of  magnitude  at  S  =>  10  dB  (representing  a  gain  of  12 
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dB) ,  and  more  than  two  orders  of  magnitude  at  S  ■  15  dB  (for  a  gain  of  more 
than  20  dB)  for  an  rms  multipath  spread  of  0.3.  The  results  in  [23]  also 
indicate  that  significant  reduction  in  error  probability  (an  order  of 
magnitude  for  S  =15  dB)  is  achieved  for  values  of  p  as  small  as  0.18.  For 
smaller  rms  multipath  spreads  (i.e.,  p  <  0.05),  the  reported  error  probability 
is  closely  approximated  by  Pg  and  is  upper  bounded  by  Pc  for  all  values  of  p 
considered.  It  is  interesting  to  note  that  the  results  in  Chapter  3  indicate 
that  for  this  channel  model  the  irreduc ible  error  probability  of  rectangular 


pulse  DPSK 

for 

p  =  0.3  is  approximately 

6xl0“2 

and 

that 

the 

average 

probability 

of 

error  for  DPSK  is  about 

10“2  at 

S 

“  15 

dB, 

even  for 

nonselective  Rayleigh  fading.  In  contrast,  no  limiting  probability  of  error 
was  exhibited  in  the  DFE  simulation  fox  error  probabilities  as  low  as  10~^ . 
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In  [24]  simulation  results  are  reported  for  an  explicit  diversity 
application  of  a  high  speed  (•»  10  Mbits/second)  troposcatter  system  described 
in  [51].  The  OFE  modem  was  designed  to  employ  rectangular  pulse  quadriphase 
shift  keying  (QPSK)  and  accommodate  up  to  four  independent  diversity  channels, 
each  with  a  forward  TDL  filter  of  three  taps  spaced  at  one-half  the  data 
symbol  interval  T.  The  outputs  of  the  forward  filters  are  coherently  combined 
and  used  as  the  input  to  in-phase  and  quadrature  filters,  each  of  which  are 
followed  by  a  three-tap  feedback  TDL  with  symbol-interval  tap  spacing.  The 
simulation  consisted  of  measuring  the  average  bit-error  probability 
performance  for  a  prototype  modem  using  a  troposcatter  •; nnel  simulator.  The 
channel  simulator  is  realized  as  a  TDL  filter  with  tap-weights  driven  by  the 
outputs  of  Gaussian  noise  generators  so  as  to  realize  a  random  channel  with  a 
fixed  rms  multipath  spread  for  each  simulation.  Using  this  configuration, 
simulation  results  are  obtained  for  values  of  rms  multipath  spreads  as  large 
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as  0.6  for  diversity  L=2  and  0.95  for  diversity  L=4 .  While  the  delay  power- 
density  spectrum  model  for  the  channel  simulator  is  not  explicitly  discussed, 
it  is  reasonable  to  assume  that  a  channel  with  a  rectangular  delay  spectrum  is 
employed  since  this  model  is  used  in  [23] ,  and  the  TDL  channel  simulator  is 
most  easily  implemented  for  this  model. 

The  results  indicate  that  for  dual  diversity  (L=2) ,  and  signal— to—noise 
ratios  of  less  than  20  dB,  the  lowest  error  probabilities  are  obtained  for  rms 
multipath  spreads  between  0.25  and  0.4.  For  these  values,  it  is  shown  that 
the  DFE  modem  produces  error  rates  which  are  lower  than  Pc(2)  by  a  factor  of 
ten  at  S  =  15  dB  (for  a  gain  of  5  dB) ,  and  an  improvement  of  nearly  two  orders 
of  magnitude  at  S  =  20  dB  (representing  a  gain  of  10  dB)  for  S  3  20  dB.  In 
contrast  to  the  results  cited  above,  this  simulation  study  establishes  the 
existence  of  an  irreducible  error  probability  of  approximately  10-6  for  p  = 
0.6  and  10  ®  for  p  *  0.45.  In  comparison,  the  results  in  Fig.  3.3  and  Fig. 
5.2  for  L  *  2,  indicate  that  for  p  *  0.45,  the  irreducible  error  probabilities 
for  both  rectangular  pulse  DPSK  and  sine  pulse  DPSK  are  about  3xl0-^ 
representing  an  equivalent  per  branch  average  signal-to-noise  ratio  of 
approximately  8  dB. 

For  fourth-order  diversity,  simulation  results  are  obtained  for  values  of 

p  as  large  as  0.95.  As  in  the  dual  diversity  simulation,  the  smallest  average 

error  probability  is  obtained  for  an  rms  multipath  spread  of  about  0.45.  In 

this  case  the  resulting  error  probability  is  shown  to  be  less  than  P  (4)  by 

c 

more  than  an  order  of  magnitude  for  S  =  10  dB  (a  gain  of  3  dB)  and  about  two 
orders  of  magnitude  for  S  ■  15  dB  (a  gain  of  5  dB) .  For  one  simulation  of 
fourth-order  diversity,  an  irreducible  error  probability  of  about  10“^  was 
obtained  for  an  rms  multipath  spread  of  0.45.  In  comparison,  the 


approximation  of  the  irreducible  error  probability  for  DPSK  with  L=4  is  about 
—3 

2x10  ,  so  that  the  use  of  DFE  and  coherent  detection  provides  a  reduction  of 
more  than  four  orders  of  magnitude  for  the  limiting  error  probability  in 
similar  fading  environments. 

5.4.2  Results  of  Experimental  Studies 

In  [14] ,  live  link  test  results  are  reported  for  the  DFE  modem  described 
in  [51]  and  used  in  the  channel  simulator  studies  cited  above.  In  this  test, 
the  DFE  modem  was  operated  in  both  dual  (L«*2)  and  quadruple  (L=4)  diversity 
configurations  and  for  two  test  "channels"  with  distinct  fading 
characteristics.  The  two  "channels"  each  consisted  of  a  radio-frequency  (RF) 
transmitter,  fading  channel,  and  an  RF  receiver  producing  an  intermediate- 
frequency  output  (forward-filter  input)  for  each  diversity  branch.  The  two 
channels  considered  in  the  test  were  the  900  MHz  and  4.5  GHz  troposcatter 
radio  systems  at  the  Rome  Air  Development  Center  test  facility.  The  multipath 
profiles  for  these  channels  were  measured  with  a  channel  probe  using  an  m- 
sequence  coded  sounding  signal.  It  was  found  that  the  delay  density  followed 
a  basic  triangular  shape  (cf.  (2.14))  and  that  the  shape  of  the  delay  density 
did  not  vary  greatly  from  test  to  test  (with  the  exception  of  the  900  MHz 
channel  where  a  secondary  triangular  "hump"  was  occasionally  observed) .  Data 
are  presented  in  [14]  for  the  900  MHz  channel  with  a  measured  rms  multipath 
spread  of  11*0.45,  and  for  the  4.5  GHz  channel  with  |t  measured  as  approximately 
0.35.  The  probability  of  error  for  the  DFE  modem  was  measured  by  the  bit¬ 
error  rate  averaged  over  twenty-minute  periods  as  a  function  of  the  median 
received  signal-to-noise  ratio.  In  what  follows,  the  performance  of  the  DFE 
modem  is  discussed  in  terms  of  the  mean  value  of  the  measured  error 
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probabilities  as  a  function  of  the  per  channel  signal-to-noise  ratio.  It 
should  be  pointed  out  that  the  results  in  [14]  exhibit  sample  points  that  vary 
from  this  mean  value  by  as  much  as  a  factor  of  ten  for  a  particular  measured 
signal-to-noise  ratio. 

For  all  cases  reported  in  [14] ,  the  test  results  show  major  discrepancies 

between  measured  DFE  modem  performance  and  the  results  of  simulation  studies 

discussed  above.  In  particular,  for  the  900  MHz  channel  and  dual  diversity, 

test  results  are  reported  for  measured  signal-to-noise  ratios  ranging  from  5 

to  17  dB.  These  results  indicate  that  the  DFE  modem  error  probability  is  on 

the  average  slightly  larger  than  Pc(2)  for  signal-t  >-noise  ratios  less  than 

about  15  dB  (differing  by  at  most  a  factor  of  three  at  S  =  5  dB,  indicating  a 

SNR  degradation  of  about  5  dB) .  For  S  =  17  dB,  the  average  error  probability 

for  the  DFE  modem  was  approximately  equal  to  P„(2).  Similar  results  were 

obtained  with  dual  diversity  for  the  4.5  GHz  channel  with  measured  rms 

multipath  spread  0.35.  For  this  channel,  test  results  are  reported  for 

signal-to-noise  ratios  ranging  from  15  to  22  dB.  The  measured  DFE  error 

probabilities  are  slightly  larger  than  P  (2)  for  S  <  20  dB,  and  differ  from 

c 

Pc(2)  by  less  than  a  factor  of  two  for  the  entire  range  of  signal-to-noise 
ratios  reported. 

The  test  results  for  fourth-order  diversity  show  that  P  (4)  is  an 

c 

excellent  approximation  of  DFE  performance  for  both  channels  considered.  In 
particular,  for  the  900  MHz  channel  (p=0.45),  the  results  show  that  the 
average  measured  DFE  error  probability  is  between  P  (4)  and  2xP  (4)  for  all 
reported  values  of  signal-to-noise  ratio  (0-15  dB)  which  corresponds  to  a 
degradation  of  less  than  3  dB.  Similarly,  test  results  for  the  4.5  GHz 
channel  (|i*0.35)  show  that  the  measured  average  error  probabilities  lie 
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between  1.2xPc(4)  and  3xPc(4)  for  ell  signal-to— noise  ratios  which  range 
between  5  and  17  dB  for  this  channel. 

The  basic  conclusion  that  can  be  drawn  from  the  above  experimental 
results  is  that  the  average  error  probability  for  maximal  ratio  PSK  is  a 
reasonable  estimate  of  the  performance  of  decision-feedback  receivers  for 
practical  fading  channel  communication  systems.  This  is  in  stark  contrast  to 
the  implications  of  the  simulation  results  which  show  that  in  certain  cases, 
the  DFE  can  realize  error  probabilities  two  orders  of  magnitude  lower  than 
PC(L),  even  in  the  case  of  no  explicit  diversity  (L“l) .  It  is  reasonable  to 
assume  that  some  of  the  discrepancy  between  predicted  and  measured  performance 
is  due  to  the  less  than  ideal  characteristics  of  the  RF  transmitter  and 
receiver  employed  for  the  experimental  studies.  However,  by  examining  some  of 
the  very  basic  limitations  of  equalizer  performance,  it  can  be  seen  that  the 
results  of  the  simulation  studies  can  be  characterized  as  "optimistic. " 

Recall,  from  the  discussion  in  Section  5.2,  that  the  feedback  portion  of 
the  DFE  receiver  is  designed  solely  to  eliminate  the  intersymbol  interference 
due  to  previously  transmitted  data  signals.  Hence,  the  implicit  diversity 
gain  can  only  be  realized  by  the  coherent  combining  properties  of  the  forward 
TDL  filter.  Thus,  the  potential  for  diversity  gain  with  respect  to  the 
decision  on  the  desired  data  bit,  b^  i*  characterized  by  the  vector  Rjsq  (cf. 
(5.24)).  For  the  example  of  a  rectangular  delay  delay  power-density  spectrum 
(2.1$)  with  parameter  Tq»T/2,  which  corresponds  to  the  channel  considered  in 
[23],  consists  of  only  two  entries,  R(»,-l)  and  R(*,0),  which  are  Gaussian 
random  variables,  linearly  related  to  the  response  function  h(»,§)  through 
(5.5).  Using  (2.16)  and  the  definition  of  the  rectangular  pulse  correlation 
function  in  (3.14a),  it  is  easy  to  see  that  the  second  moments  of  the 
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components  of  R^q  axe  related  by 

14.E{R2(.,-1)}  -  E{R2(.,0)}  . 

In  general,  the  random  variables  R(«,-l)  and  R(,,0)  are  correlated  Gaussian 
random  variables  by  nature  of  the  overlapping  portions  of  the  rectangular 
pulse  autocorrelation  function.  However,  in  order  to  characterize  a  lower 
bound  for  DFE  receiver  performance,  we  proceed  under  the  assumption  that 
R(*,-l)  and  R(»,0)  are  uncorrelated  and  hence  independent.  (It  is  well  known 
16]  that  correlations  among  diversity  links  can  only  degrade  the  performance 
of  maximal  ratio  diversity  systems.)  In  this  case,  the  error  probability  for 
maximal  ratio  PSK  is  given  by  (5.30)  by  neglecting  the  parameter  L,  where  the 
probability  density  function  of  the  instantaneous  signal-to-noise  ratio  s  is 
given  by 

f(s)  =  i3s  »[exp(-s/S)  -  exp(-14s/S)]  , 

and  where  S  is  the  average  signal-to-noise  ratio  corresponding  to  R(.,0). 

Calculation  of  the  integral  in  (5.30)  for  this  density  function  of  s  shows 

that  for  S  between  10  and  30  dB,  the  maximal  ratio  error  probability  lies 

approximately  mid-way  between  the  curves  corresponding  to  P  (i)  and  P  (2) 

c  c 

shown  in  Fig.  5.2,  and  obtains  a  value  of  5.4xl0~^  for  S  =  30  dB.  Hence,  we 
see  that  even  under  ideal  conditions,  the  effective  order  of  implicit 
diversity  is  between  one  and  two.  Moreover,  it  is  not  likely  that  the  actual 
DFE  performance  would  reflect  much  of  this  potential  diversity  gain  since  i) 
in  reality,  the  implicit  diversity  branches  are  correlated,  ii)  the  one-shot 
MUSE  forward  filter  approaches  the  maximal  ratio  combining  filter  only 
asymptotically  for  small  signal-to-noise  ratios,  and  iii)  the  MMSE  forward 
filter  tap  weights  are  largely  determined  by  the  character  of  the  intersymbol 


interference,  especially  for  high  signal- to-noise  ratios.  Thus,  second-order 
diversity  may  be  considered  to  be  a  gross  upper  bound  for  the  effective  order 
of  implicit  diversity  for  this  example.  In  contrast,  it  is  claimed  in  [23], 
that  results  of  the  simulation  study  for  this  channel  model  show  that  the  DFE 
receiver  realizes  an  implicit  diversity  order  of  approximately  four.  The 
discussion  above  and  the  results  of  experimental  studies  in  [14]  show  that 
this  is  not  the  case. 

We  conclude  that  while  some  potential  for  implicit  diversity  gain  does 
exist,  it  is  not  a  significant  factor  for  MMSE  decision-feedback  equalization. 
Rather,  the  main  attributes  of  decision-feedback  equalizers  are  the  ability  to 
establish  a  coherent  communication  environment,  and  the  substantial  reduction 
of  the  effects  of  intersymbol  interference  for  practical  systems  within  the 
limits  of  the  capabilities  of  the  equalizer.  Thus,  for  a  given  order  of 
explicit  diversity  L,  the  performance  of  DFE  receivers  in  frequency-selective 
Rayleigh  fading  channels  is  both  easily  and  accurately  approximated  by  the 
average  error  probability  for  maximal  ratio  PSK  defined  in  (5.31). 

5.5  Additional  Considerations 


There  are. 

of  course,  a  number  of  design 

problems 

associated  with 

the 

implementation 

of 

adaptive 

equalizers 

which  largely  determine 

the 

effectiveness  of 

the 

equalizer. 

The  aspect  of 

equal izer 

implementation  which 

moat  profoundly  affects  system  performance  is  the  realization  of  the  algorithm 
used  to  adaptively  adjust  the  coefficients  of  the  TDL  filters.  In  practice, 
the  equalizer,  and,  hence,  the  tap-adjustment  hardware,  is  designed  to  operate 
in  two  modes.  At  the  start  of  operation,  a  known  data  sequence  or  training 
sequence  is  transmitted  to  allow  the  receiver  to  adjust  to  the  initial  state 
of  the  channel.  The  required  length  of  the  training  sequence  depends  on  the 


« 


total  number  of  equalizer  taps,  the  characteristics  of  the  channel,  and  the 
speed  of  convergence  of  the  adjustment  algorithm  [23,44,52].  At  the  end  of 
the  start-up  mode,  the  receiver  uses  the  results  of  previous  data  symbol 
decisions  to  continuously  update  the  tap  coefficients  [40,44].  Thus,  the 
primary  performance  criteria  for  tap-adjustment  algorithms  are  the  speed  of 
convergence,  the  sensitivity  to  occasional  decision  errors,  and  the  complexity 
or  cost  of  implementation. 

The  most  commonly  implemented  adjustment  algorithm  for  both  linear  1DL 
and  decisionr-feedback  HSE  equalizers  is  a  form  of  stochastic  gradient  or 
steepest  descent  algorithm  [23,49,50].  In  this  implementation,  the  tap-vector 
w  is  recursively  adjusted  according  to  the  relation 

-k+1  =  *k  +  AefcEk  ' 

where  0^  is  the  vector  of  received  samples  U  given  by  (5.3),  ek  is  the 
difference  between  the  k-th  output  and  the  decision  on  the  k-th  output,  and  A 
is  a  scale  factor  which  must  be  chosen  small  enough  to  ensure  equalizer 
convergence  [44].  While  the  steepest  descent  algorithm  has  the  advantage  of 
ease  of  implementation,  it  is  characterized  by  a  slow  rate  of  convergence 
[47].  It  has  been  shown  [53]  that  steepest  descent  algorithms  require 
training  sequences  of  length  lOxN,  where  N  is  the  total  number  of  equalizer 
taps.  A  primary  constraint  for  the  equalization  of  fading  channels  is  the 
need  to  track  the  sometimes  rapidly  varying  channel.  For  applications  to 
Rayleigh  fading  channels,  it  is  necessary  to  insert  training  sequences  at 
periodic  intervals  to  ensure  stable  equalizer  operation.  For  example,  in  [8] 
the  results  of  a  simulation  study  of  a  DFE  consisting  of  ten  forward  and  ten 
feedback  taps  with  no  explicit  diversity  are  presented.  The  study  was 


implemented  using  a  channel  simulator  to  realize  a  channel  with  average 
decorrelation  time  on  the  order  of  10^  data  pulse  durations.  In  order  to 
maintain  channel  equalization,  it  was  necessary  to  insert  a  100-bit  training 
sequence  at  2000-bit  intervals  in  the  data  sequence. 

The  slow  rate  of  convergence -of  the  steepest  descent  algorithm  is  a 
limiting  factor  in  determining  the  applicability  of  adaptive  equalizers 
implementing  this  algorithm.  Vithin  the  last  decade,  a  number  of  studies  have 
been  directed  toward  finding  rapidly  converging  algorithms  that  can  be 
implemented  with  minimum  additional  complexity  [47,53-56].  Of  these,  the 
algorithms  that  have  received  the  most  attention  are  those  derived  from  the 
work  of  Godard  [53]  who  first  applied  the  Kalman  estimation  algorithm  to  the 
problem  of  equalizer  adjustment.  These  investigations  show  that  it  is 
possible  to  obtain  near  MSE-optimal  equalizer  adjustment  using  training 
sequences  of  length  N,  where  N  is  the  total  number  of  equalizer  taps  [54]. 
Moreover,  it  is  suggested  that  equalizer  adjustment  algorithms  of  this  type 
are  capable  of  tracking  the  occasionally  rapid  changes  experienced  in  a  fading 
environment.  The  two  main  disadvantages  of  these  algorithms  are  additional 
sensitivity  to  additive  noise  [44]  and  increased  complexity  [44].  The 
additional  computational  burden  for  even  the  most  efficient  of  these 
algorithms  is  roughly  2N  tines  the  complexity  of  steepest  descent  algorithms 
for  an  N-tap  equalizer. 


CHAPTER  6 


SUMMARY  AND  CONCLUSIONS 

In  this  thesis,  we  have  investigated  the  effects  of  frequency-selective 
fading  on  the  average  probability  of  error  for  DPSK  and  FSK  communications. 
By  considering  the  effects  of  frequency-selective  fading  in  both  the  frequency 
domain  and  the  delay  domain,  we  have  examined  the  interaction  between  the 
characteristics  of  the  fading  channel,  the  transmitted  signal,  and  the 
quadratic  nature  of  the  detection  process  employed  for  DPSK  and  noncoherent 
FSK  communications.  Ve  have  identified  the  system  parameters  that  are  of 
fundamental  importance  in  determining  the  average  error  probability.  It  has 
been  shown  that  the  performances  of  DPSK  and  FSK  can  be  approximated  in  terms 
of  one  or  two  parameters  that  can  be  obtained  from  rms-type  channel 
measurements.  Lastly,  we  have  examined  the  basic  limitations  of  equalizer 
performance  and  have  developed  a  method  of  obtaining  estimates  of  adaptive 
equalizer  performance  for  practical  systems. 

The  analysis  of  DPSK  and  FSK  for  several  examples  of  multipath  models 
(delay  power-density  spectra)  leads  to  several  conclusions.  For  DPSK  systems, 
the  average  error  probability  as  well  as  the  irreducible  error  probability 
strongly  depend  on  the  shape  of  the  data-pulse  waveform.  These  examples  also 
indicate  that  the  lowest  error  probabilities  are  obtained  when  small  time- 
bandwidth  product  pulses  are  employed.  In  particular,  the  sine-pulse,  which 
is  the  smallest  time-bandwidth  product  pulse  waveform  considered  [39],  also 
produces  the  lowest  error  probability  of  any  pulse  shape  considered.  Large 
time-bandwidth  product  pulse  shapes  (e.g.,  phase-coded  pulses  generated  by  nr- 
sequences)  exhibit  increased  sensitivity  to  the  effects  of  intersymbol 


interference  as  well  as  larger  error  probabilities  for  practical  values  of 
signal-to-noise  ratio. 

Similar  conclusions  apply  to  the  relationship  between  the  data-pulse 
shape  and  the  probability  of  error  for  FSK  systems.  In  addition,  the 
modulation  index  as  well  as  the  relative  phases  between  successive  transmitted 
signals  play  an  important  role  in  determining  the  system  error  probability. 
The  lowest  error  rates  for  FSK  are  obtained  for  continuous-phase  signals  with 
small  modulation  indices. 

Comparing  the  results  for  DPSK  and  FSK,  we  find  that  the  error 
probabilities  for  these  systems  depend  on  a  number  of  common  factors,  and  in 
certain  cases,  the  effects  of  frequency-selectivity  on  both  DPSK  and  FSK 
systems  can  be  characterized  by  the  same  parameter.  A  technique  for  obtaining 
bounds  on  system  performance  in  terms  of  the  hey  system  and  channel  parameters 
applies  to  the  evaluation  of  both  DPSK  and  FSK  systems  for  the  several 
examples  of  channel  models  and  signaling  formats.  Finally,  the  error 
probabilities  for  systems  employing  complicated  pulse  shapes  and  for  channels 
which  are  difficult  to  fully  characterize  can  be  predicted  from  the  hey  system 
parameters  and  the  performance  data  contained  in  this  thesis. 

For  channels  with  rms  multipath  spreads  less  than  about  0.1,  it  i. 
possible  to  substantially  reduce  the  irreducible  error  probabilities  for  both 
DPSK  and  FSK  systems  by  choosing  system  parameters  which  minimize  the 
sensitivity  to  intersymbol  interference.  However,  for  many  practical  systems 
the  effects  of  1SI  can  severely  limit  the  performance  of  these  forms  of 
digital  communications.  In  Chapter  5,  the  characteristics  of  adaptive 
equalizers  coauaonly  employed  for  fading-channel  communications  are  described 
and  the  results  of  several  simulation  and  experimental  studies  are  compared. 


It  is  revealed,  that  while  claims  of  "implicit  diversity  gain"  abonnd  in  the 
literature,  the  significant  attributes  of  adaptive  equalization  techniques  art 
the  ability  to  establish  a  coherent  communication  link  and  the  substantial 
reduction  of  the  effects  of  intersymbol  interference.  Comparing  the  results 
of  simulation  studies  and  experimental  investigations  of  equalizer  performance 
with  the  analytical  results  for  the  error  probabilities  of  DPSK  and  FSK 
systems  reveals  that  reliable  communications  can  be  achieved  through  channel 
equalization  even  in  fading  environments  that  produce  unacceptably  high  error 
probabilities  for  unequalized  systems. 

It  should  be  pointed  out  that  one  of  the  factors  limiting  the  implicit 
diversity  gain  for  signals  with  small  time-bandwidth  products  is  the  inability 
to  resolve  the  multipath  components  in  the  received  signal.  For  rectangular 
data  pulse  waveforms  with  duration  T,  it  is  clear  that  signals  arriving  at  the 
receiver  with  relative  delays  of  less  than  T  seconds  cannot  be  resolved  into 
independently  faded  components.  In  fact,  for  YSSDS  fading  channels,  the 
statistical  correlation  between  signals  received  with  relative  delay  £  is 
completely  determined  by  the  product  of  the  delay  power-density  spectrum  and 
the  autocorrelation  function  of  the  data  pulse  waveform  evaluated  at  £. 

Within  the  last  five  years,  there  has  been  considerable  interest  in  the 
application  of  equalization  techniques  to  both  coherent  and  noncoherent 
spread-spectrum  communication  systems  employing  phase-coded  waveforms.  The 
analysis  [57]  and  test  results  [58,59]  indicate  that  the  use  of  this  type  of 
signal  with  large  time-bandwidth  product  allows  resolution  and  coherent 
(maximal  ratio)  combining  for  both  HF  and  troposcatter  channels.  The 
development  of  useful  bounds  and  approximations  of  the  performance  of  adaptive 
equalizers  for  both  spread-spectrum  and  conventional  communications  in 


frequency-selective  fading  channels  is  an  area  where  much  further  work  is 
needed. 

The  current  interest  in  applications  of  adaptive  equalization  techniques 
to  digital  communications  over  frequency-selective  fading  channels  is  in  large 
part  due  to  the  well-documented  effectiveness  of  these  techniques  in  nonfading 
systems.  While  it  is  true  that  reliable  communications  can  be  achieved 
through  the  use  of  adaptive  equalizers,  it  is  not  clear  that  channel 
equalization  is  the  most  efficient  way  to  reduce  the  effects  of  intersymbol 
interference.  Unfortunately,  the  published  literature  is  void  of  the 
consideration  of  alternative  techniques  that  could  be  used  in  practical 
systems  to  obtain  acceptable  performance. 

One  example  of  a  possible  alternative  which  deserves  consideration  is  the 
use  of  St-ary  orthogonal  signalling  in  conjunction  with  error  correcting  codes. 
In  particular,  suppose  that  the  source  data  rate  R  -  10^  bits/ sec  and  the 
channel  rms  delay  M  =  3x10  ^  sec  are  fixed  design  parameters.  Hence,  for 
binary  signaling  the  rms  multipath  spread  p^  =  M/T  =  0.3.  If  8-ary  orthogonal 
FSK  is  used  to  transmit  the  binary  data,  the  effective  rms  multipath  spread  is 
reduced  to  pg  *  0.1.  Unfortunately,  the  evaluation  of  the  probability  of 
error  for  K-ary  orthogonal  FSK  in  frequency-selective  Rayleigh  fading  is  at 
present  an  unsolved  problem.  While  it  is  true  that  the  introduction  of  an 
error-correcting  code  has  the  effect  of  increasing  the  channel  data  rate  (and, 
hence,  increasing  the  effective  rms  multipath  spread),  it  may  be  possible  to 
realize  error  rates  significantly  lower  than  those  achieved  through  adaptive 
channel  equalization.  The  investigation  of  alternative  methods  of  achieving 
reliable  communications  over  frequency-selective  fading  channels  is  an  area 
which  appears  to  be  both  promising  and  challenging. 
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